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Abstract 

I  present  an  in-depth  theoretical  analysis  of  the  optical  impact  of  transverse  index 
confinement  on  the  lasing  modes  of  microcavity  surface  emitting  lasers.  Using  two  differ¬ 
ent  variational  techniques — the  weighted  index  method  and  the  finite  element  method  I 
solve  Maxwell’s  equations  for  several  azimuthally  symmetric  laser  designs  to  compute  the 
relevant  laser  parameters:  modal  resonance,  field  profile,  confinement  factor,  and  threshold 
gain.  Through  my  weighted  index  analysis,  I  discovered  two  new  effects  of  transverse  field 
confinement:  a  polarization  dependent  change  in  mirror  reflectance  with  aperture  radius, 
previously  only  noted  in  edge-emitting  lasers,  and  a  mode  dependent  blueshift.  Comparing 
my  blueshift  predictions  against  measured  results  for  devices  operating  at  nominally  850  nm 
and  780  nm,  I  found  excellent  agreement  between  theory  and  experiment.  For  a  sufficiently 
sharp  gain  spectrum,  this  blueshift  may  be  used  to  achieve  low  threshold,  single  mode  lasing 
by  exploiting  the  dispersion  in  blueshift  rates  amongst  the  various  transverse  eigenmodes  to 
spectrally  select  the  single  desired  lasing  mode.  Through  my  finite  element  analysis,  I  ascer¬ 
tained  the  physical  mechanisms  responsible  for  diffraction  in  microcavity  lasers.  These  are: 
the  transverse  confinement  of  the  optical  mode,  the  degree  of  tilt  of  the  mode  propagation 
vector  away  from  normal  to  the  mirror  surfaces,  and  the  density  of  parasitic  modes  in  the 
spectral  vicinity  of  the  lasing  mode.  This  knowledge  is  of  vital  importance  for  low  threshold 
design  analyses,  since  diffraction  is  the  dominant  loss  source  in  this  regime.  Based  on  the 
success  of  the  parasitic  mode  loss  picture  for  describing  diffraction  in  the  stratified  surface 
emitting  laser  geometry,  I  augmented  my  weighted  index  model,  via  a  rigorous  electromag¬ 
netic  calculation  of  lasing-to-parasitic  modes  coupling,  to  compute  diffractive  loss  and  total 
threshold  gain.  Incorporating  this  new  weighted  index/parasitic  mode  calculation  with  a 
semiconductor  bandstructure  and  gain  model,  I  computed  spontaneous  emission  rates  for  a 
low  threshold  laser  using  a  new  optical  mode  density  model  specifically  designed  to  address 
the  partial  field  confinement  afforded  by,  for  example,  oxide  apertures.  Finally,  the  optical, 
gain,  and  emission  results  were  combined  to  compute  light  verses  current,  threshold  current, 
and  spontaneous  emission  factors  for  the  same  low  threshold  device. 
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Optical  Physics  of  Microcavity  Surface  Emitting  Lasers 


L  Introduction 

1.1  Foreword 

The  advances  in  semiconductor  and  optoelectronic  technology  over  the  last  quarter 
century  have  opened  the  door  for  the  theoretical  researcher  to  probe  concepts  not  only  of 
purely  scientific  interest  but  of  great  practical  application  as  well.  This  is  perhaps  best  exem¬ 
plified  by  the  archetypical  semiconductor  laser,  whose  understanding  requires  the  mastery  of 
several  branches  of  physics  and  yet  boasts  wide  sweeping  application  ranging  from  compact 
disc  players,  to  environmental  sensors,  to  aircraft  avionics,  etc.  These  lasers  have  established 
a  growing  strong-hold  within  the  optical  source  genre  due  in  a  large  part  to  their  cost,  size, 
and  exceptional  operating  characteristics,  many  of  which  are  unobtainable  by  their  gas  and 
solid  state  counterparts — at  any  cost.  Indeed,  as  anyone  who  has  worked  with  them  can 
attest,  the  idea  of  a  ND^+-YAG,  ruby,  or  CO2  laser  operating  with  an  efficiency  of  50% 
off  batteries  and  costing  pennies  on  the  dollar  belongs  in  an  Arthur  C.  Clarke  novel.  Yet 
semiconductor  lasers  have  achieved  all  these  feats  and  more.  A  prime  example  is  amplitude 
modulation  by  direct  modulation  of  the  drive  current,  a  feature  ideally  suiting  these  lasers 
for  fiber-optic  sources,  the  backbone  of  the  communications  industry.  The  prominence  of 
these  lasers  is  largely  attributable  to  the  advances  in  semiconductor  growth  and  microfab¬ 
rication  technology,  permitting  the  researcher  to  control  device  construction  literally  down 
to  the  atomic  level  and  to  realize  truly  quantum  mechanical  structures — as  Steven  Wright 
quipped — “the  dreams  stuff  are  made  of.”  Most  notably,  this  technology  facilitated  the 
creation  of  hetero junctions  and  quantum  wells,  the  optical  sources  within  most  lasers,  and 
until  then  a  purely  theoretical  entity  of  quantum  mechanics.  Despite  the  relative  maturity 
of  the  semiconductor  laser,  the  research  activity  has  not  slowed,  maintaining  a  frenzied  pace 
driven  by  the  continual  emergence  of  new  applications  and  the  demand  for  increased  per¬ 
formance.  A  large  percentage  of  the  work  is  focused  on  conquering  new  material  systems 
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to  achieve  reliable  lasers  at  wavelengths  outside  of  the  well  developed  red-to-near-infrared 
range.  The  elusive  blue  laser,  the  last  leg  of  the  triad  required  for  high-brightness,  flat-panel 
displays,  and  the  mid-infrared  intraband  quantum  cascade  laser,  which  uses  the  well  estab¬ 
lished  GaAs/AlGaAs  material  system,  are  excellent  examples  of  this  type  of  work.  Although 
this  research  is  exceptionally  significant  in  the  material  science  and  quantum  mechanics  dis¬ 
ciplines,  optically  these  lasers  differ  very  little  from  their  predecessors.  This  is  not  the  case 
though  with  vertical  cavity  surface  emitting  lasers,  or  VCSELs,  another  highly  active  area 
of  laser  research,  and  the  foundation  of  this  thesis. 

The  VCSEL  is  a  relatively  recent  spin-off  of  the  more  traditional  edge-emitting  laser, 
whose  primary  differences — cavity  orientation,  size,  and  design — are  optical  in  nature;  a  typ¬ 
ical  VCSEL  and  edge  emitter  are  illustrated  in  Figure  1.1.  As  the  name  implies,  the  VCSEL 
cavity  is  normal  to  the  growth  direction.  VCSEL  mirrors  are  either  monolithically  grown 
semiconductor  or  post- growth  deposited  dielectric  distributed  Bragg  reflectors  (DBRs):  al¬ 
ternating  quarter-wave  layers  of  high  and  low  index  of  refraction  materials.  A  drawback  of 
the  vertical  cavity  is  a  very  short  gain  region,  often  necessitating  extremely  thick  DBRs  of 
80  layers  or  more  to  achieve  the  greater  than  99%  reflection  required  for  lasing.  On  the 
other  hand,  however,  the  typical  VCSEL  cavity  is  only  a  few  wavelengths  long,  roughly  two 
orders  of  magnitude  shorter  than  most  in-plane  lasers.  An  immediate  result  being  the  lasing 
threshold,  which  scales  roughly  with  active  volume,  is  also  orders  of  magnitude  smaller,  un¬ 
fortunately  at  the  expense  of  the  ability  to  generate  high  output  powers.  Additionally  due  to 
the  short  cavity  length,  VCSELs  only  support  a  single  longitudinal  mode.  The  large  spacing 
between  spectrally  sharp  longitudinal  modes — due  to  the  high  reflectance  DBRs — makes  the 
VCSEL  cavity  the  optical  analog  of  the  quantum  well,  a  quasi-two  dimensional  structure. 
This  feature  has  been  thoroughly  investigated,  focusing  on  the  potential  for  cavity  design  to 
effect  the  quantum  nature  of  light,  enhancing  or  suppressing  spontaneous  emission  to  lower 
the  lasing  threshold  beyond  the  active  volume  limit.  From  the  results  of  analyses  on  the 
simple  planar  VCSEL  came  predictions  that  with  further  optical  confinement,  lower  order 
quantum  systems  could  be  imitated  and  lower  thresholds  obtained. 
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Figure  1.1 


Illustration  of  a  simple  etched  post  vertical  cavity  surface  emitting  laser  and 
an  edge-emitting  laser. 
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The  first  embodiment  of  a  three-dimensionally  index  confined  vertical  cavity  laser 
came  in  the  form  of  an  etched  post,  with  DBRs  providing  one  dimension  of  confinement 
and  the  air-to-semiconductor  interface  providing  the  other  two.  Although  it  would  seem 
this  device  would  provide  the  closest  replica  of  the  classical  “particle  in  a  box”,  it  was 
not  practical.  It  suffered  terribly  from  surface  recombination  electrical  losses,  as  well  as 
optical  scattering  losses,  and  was  therefore  quickly  abandoned.  Nonetheless,  in  its  short  life¬ 
span  a  great  deal  of  analysis  was  performed  on  it,  with  the  most  noteworthy  result  being  a 
prediction  that  in  a  true  zero  dimensional  system  the  threshold  could  be  made  to  completely 
disappear.  The  potential  of  this  device  is  enormous;  a  zero  threshold  device  would  require 
zero  bias  power.  This  untapped  potential  served  to  motivate  researchers  both  interested 
in  the  practical  application  of  such  a  device  as  well  as  those  intrigued  by  the  exciting  new 
physics  embodied  in  them.  The  breakthrough  came  when  the  microelectronics  group  at 
the  University  of  Texas  (U.T.)  found  that  by  replacing  proton  implantation  with  native 
oxidation  for  current  confinement,  they  achieved  the  additional  benefit  of  index  confinement 
of  the  optical  field  without  the  immense  electrical  losses  associated  with  etching  [31].  The 
idea  caught  on  like  wildfire  and  over  the  few  short  years  since  its  inception,  oxidation  has 
essentially  replaced  implantation  as  the  method  of  choice  for  carrier  and  field  confinement. 

However,  the  oxide  is  not  perfect.  It  clearly  provides  a  transverse  index  step,  but 
the  total  effect  of  the  step  is  not  well  understood,  as  evidenced  by  the  large  variety  of 
oxide  designs  throughout  the  research  community.  For  example,  two  leading  universities  in 
VCSEL  research,  the  University  of  California  at  Santa  Barbara  (U.C.S.B.)  and  the  University 
of  Texas,  report  optimized  designs  based  on  wholly  different  concepts  of  the  oxide  optical 
effect.  Santa  Barbara  reports  that  thinner  oxides  placed  further  from  the  gain  region  at 
nodes  of  the  field  standing  wave  are  superior,  minimizing  the  diffraction  loss  of  the  lasing 
mode  [30].  Texas  reports  just  the  opposite;  thicker  oxides  placed  near  the  gain  region  at 
antinodes  are  better,  providing  more  transverse  confinement/waveguiding  to  the  mode  [19]. 
Both  groups  have  enjoyed  success  with  their  respective  approaches,  achieving  extremely  low 
thresholds  [33] .  Their  results  suggest  that  the  primary  effects  of  the  oxide  differ  depending 
on  the  cavity  and  oxide  aperture  design.  These  results  may  also  be  interpreted  to  mean 
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the  use  of  oxides  has  not  yet  been  optimized  and  the  ultimate  threshold  limit  has  not  yet 
been  achieved.  A  further  reduction  in  threshold  must  be  preceded  by  a  better  understanding 
of  the  oxide  effects.  A  complete  model  must  explain  these  two  seemingly  opposing  results. 
Once  in  hand,  this  “Grand  Unified  VCSEL  Optical  Theory”  may  be  used  to  predict  the  ideal 
oxide  aperture  and  VCSEL  cavity  design  for  achieving  the  ultimate  limit  in  low  threshold 
lasing.  It  is  precisely  this  model  which  is  the  goal  of  this  dissertation. 

l.S  Optical  Problem 

The  fact  that  the  optical  problem  has  persisted  so  long  reflects  the  extreme  difficulty 
in  solving  the  differential  equations  governing  the  fields  in  VCSELs.  The  equations  are 
themselves  simple  enough.  They  are,  of  course.  Maxwell’s  equations,  given  in  differential 
form  (in  MKS  units)  as 

VxE  =  (1.1) 

-  dD  t 

Vxff  =  ^  +  J. 

V-.D  =  /?„, 

VB  =  0. 

Here  the  two  source  terms  are  the  current  density  and  the  volume  charge  density  (/>,,). 
Assuming  isotropic  media,  the  electric  magnetic  fluxes  are  related  to  the 

electric  and  magnetic  fields  by  the  material  parameters  {p,  e)  and  the  constitutive 
relations 

D  =  eoSrE,  (1.2) 

B  =  (XQiiril-, 

Pq,  Eo,  and  are  the  free-space  and  relative  permeability,  and  the  free-space  and  relative 
permittivity,  respectively.  The  relative  permittivity  is  related  to  the  more  commonly  used 
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index  of  refraction  by 


6r  —  T)  . 


(1.3) 


Absorptive  loss  and  material  gain  result  from  a  negative  or  positive  imaginary  part  of  or 
rj.  Since  lasing  fields  are  spectrally  sharp— nearly  delta  functions  in  frequency— the  fields 
may  be  approximated  as  harmonic,  depending  on  time  as  exp  (iut),  where  i  =  y/-^  and  u)  is 
the  radial  frequency.  Furthermore,  assuming  non-magnetic  material  (//,.  =  1),  which  is  the 
case  for  all  semiconductor  laser  materials,  the  Maxwell  curl  equations  may  be  combined  to 
derive  the  vector  Helmholtz  equation, 

V  X  V  X  ^  —  (J^HQeQCrE  =  —iujfioJ .  (1-4) 

The  vector  Helmholtz  equation  and  the  MKS  Maxwell’s  equations  serve  as  the  foundation 
for  all  my  optical  laser  analysis.  The  primary  challenge  in  solving  (1.4)  results  from  the 
apertured  VCSEL  geometry.  The  geometry  introduces  two  fundamental  difficulties:  (1)  the 
relative  permittivity,  and  therefore  the  vector  Helmholtz  equation,  is  not  separable  due  to 
the  aperture  or  etched  post;  (2)  the  planar  Fabry-Perot  cavity  has  very  high  Q.  The  non¬ 
separability  is  intrinsic  to  the  apertured  VCSEL  design;  without  it  there  could  be  no  index 
confinement^  Similarly,  the  high  Q  is  vital  for  lasing  to  occur  in  the  short  VCSEL  cavities. 
These  two  intrinsic  VCSEL  attributes  conspire  to  greatly  complicate  optical  modeling  efforts. 


1.3  Semianalytic  Calculations 

Since  analytic  calculations  generally  yield  not  only  a  specific  solution  to  a  given  prob¬ 
lem,  but  insight  into  the  physical  underpinnings  of  the  problem,  they  are  often  preferred 
over  numerical  techniques.  Moreover,  analytic  calculations  are  as  a  rule  much  faster  than 
numerical  approaches,  favoring  them  for  iterative  design  work.  Hence,  the  majority  of  the 
optical  modeling  efforts  have  involved  an  analytic  or  semianalytic  representation  of  the  fields. 
The  most  obvious  approach  is  to  actually  solve  the  full  vector  Helmholtz  equation  in  each 
region  of  constant  index  of  refraction,  coupling  the  solutions  via  the  tangential  field  continu- 

^This  is  not  strictly  true,  but  does  hold  for  all  current  VCSEL  designs.  This  issue  will  be  revisited  later 
in  Chapter  V. 
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ity  conditions  [12].  In  this  approach  the  field  is  represented  by  an  infinite  weighted  sum  of 
eigenfunctions  (for  example,  trigonometric  or  Bessel  functions)  for  each  region,  the  weight¬ 
ing  determined  by  mode  matching  to  adjacent  regions.  But  this  technique  works  well  only 
when  the  refractive  index  profile  separates  in  some  preferred  coordinate  system,  reducing  the 
infinite  expansions  to  a  single  term.  If,  as  in  realistic  VCSEL  structures,  the  refractive  index 
profile  does  not  separate  in  any  coordinate  system,  then  simple  single-term  special-function 
solutions  to  the  governing  partial  differential  equations  do  not  represent  exact  solutions  for 
the  modes  and  analytic  methods  become  quite  cumbersome.  Despite  this  difficulty,  most 
previous  analytic  calculations  have  introduced,  at  some  point  in  the  treatment,  a  single 
product  term  to  describe  a  particular  electromagnetic  field  component  [21,26].  It  is  vital  to 
realize  that  this  is  equivalent  to  assuming  the  underlying  differential  equation  separates. 

Since  separable  descriptions  facilitate  closed-form  expressions,  rapid  calculation,  and 
comparison  with  well-understood  “textbook”  problems,  there  is  considerable  motivation  to 
improve  and  justify  them.  Therefore,  in  Chapter  II,  I  generalize  the  weighted  index  method 
(WIM) — a  separable  approximation — to  compute  cavity  modes  in  cylindrically-symmetric 
dielectric  VCSEL  structures.  I  also  show  (in  Appendix  A),  using  the  calculus  of  variations, 
that  this  technique  provides  the  best  separable  solution  to  the  scalar  Helmholtz  equation.  In 
Chapter  III,  I  apply  the  WIM  to  both  oxide-apertured  and  etched-post  VCSEL  designs  to 
calculate  the  spatial  profile,  optical  confinement  factor,  resonant  frequency,  and  mirror  and 
absorption  losses  of  the  fundamental  and  higher-order  lasing  modes.  Furthermore,  I  compare 
the  WIM  lasing  mode  resonance  predictions  to  measurement  for  two  different  apertured 
device  designs,  validating  the  WIM  results.  Ignoring  diffraction  effects,  I  show  in  Appendix 
B  that  the  threshold  conditions  based  on  the  WIM  optical  losses  are  a  generalization  of  the 
classical  Fabry- Perot  round-trip  amplitude  and  phase  conditions  [65]. 

1-4  Numerical  Solutions 

For  aperture  sizes  approaching  the  field  wavelength,  diffraction  becomes  the  dominant 
loss  mechanism,  rendering  the  original  WIM  incomplete.  In  this  realm,  any  threshold  pre¬ 
diction  based  solely  on  the  WIM  losses  will  be  inaccurate — the  smaller  the  aperture,  the 
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less  accurate  the  results.  I  would  prefer  to  augment  the  WIM  to  incorporate  diffraction  and 
make  it  a  complete  optical  model.  But,  before  this  can  be  accomplished  it  is  crucial  to  un¬ 
derstand  the  true  physical  mechanisms  which  determine  the  diffraction  in  VCSELs.  Ideally, 
this  could  be  done  experimentally.  Several  devices  of  different  oxide  aperture  radii  could  be 
fabricated  and  the  threshold  current  measured.  The  threshold  current,  however,  strongly 
depends  on  the  electrical  parasitics — surface  recombination  and  leakage  current — both  of 
which  vary  significantly  for  small  aperture  VCSELs.  Hence,  it  is  extremely  difficult  to  ferret 
out  optical  loss  information  from  threshold  current  measurements,  necessitating  the  use  of 
modeling  to  understand  the  purely  optical  effects  of  oxide  confinement. 

Volumetric  numerical  methods,  although  computationally  intensive,  are  ideally  suited 
for  rigorous  analysis  of  non-separable  partial  differential  equations  (PDEs)  over  a  closed  do¬ 
main.  These  techniques  approximate  the  true  answer  by  a  numerical  one,  formed  by  discretiz¬ 
ing  either  the  PDE — ^by  replacing  the  derivatives  with  finite  differences — or  the  solution — by 
expanding  the  function  over  a  weighted  basis  set.  These  methods  are  quasi-exact,  in  that 
their  accuracy  increases  with  the  number  of  unknowns  and  becomes  exact  in  the  infinite 
limit.  To  date,  only  one  numerical  method  has  been  applied  to  the  VCSEL  optical  prob¬ 
lem.  G.R.  Hadley  [27]  estimated  the  diffractive  losses  due  to  oxide  apertures  using  a  finite 
difference  eigensolver.  But  his  solution  was  limited  in  two  primary  respects:  (1)  the  finite 
difference  approach  only  works  on  a  Cartesian  grid,  and  (2)  it  may  only  be  used  to  solve  for 
scalar  fields.  His  model  can  not  handle  tapered  oxides,  touted  by  Santa  Barbara  and  Sandia 
National  Laboratories  as  the  optimal  oxide  design.  Furthermore,  for  VCSEL  designs  with 
oxides  in  or  near  the  cavity — such  as  the  University  of  Texas’s  design — the  vector  nature  of 
the  fields  may  be  important.  To  extend  his  results,  in  Chapter  IV,  I  derive  a  new  vector  finite 
element  method  (FEM).  The  FEM  explicitly  solves  for  all  three  electric  field  vector  compo¬ 
nents  in  azimuthally  symmetric  VCSELs^.  The  method  is  based  on  a  variational  solution 
of  the  full  vector  Helmholtz  equation.  Unlike  finite  difference  approaches,  the  FEM  is  not 
limited  to  a  Cartesian  grid  and  can  easily  incorporate  any  azimuthally  symmetric  dielectric 
geometry,  including  complex  loss.  In  this  method,  the  Helmholtz  equation  boundary  value 

present  the  vector  finite  element  functions  in  Appendix  C 
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problem  is  transformed  to  a  large,  sparse,  non-Hermitian,  generalized  eigenvalue  problem. 
For  each  cavity  mode  the  eigenvalues  correspond  to  the  complex  radial  frequency  (a;)  and 
the  eigenvectors  are  the  basis  expansion  coefficients  determining  the  fields.  Since  the  method 
may  only  be  applied  on  a  closed  domain,  in  Appendix  D,  I  derive  a  new  artificial  absorbing 
layer  designed  to  mimic  the  unbounded  nature  of  the  fields.  I  conclude  Chapter  IV  by  ap¬ 
plying  the  FEM  to  a  novel  oxide  apertured,  oxide  DBR  VCSEL,  reported  to  be  the  smallest 
working  (photopumped)  device  ever  built  [67].  From  the  eigenvalue  and  eigenvector,  I  derive 
all  the  pertinent  mode  data:  resonant  wavelength,  confinement  factor,  field  distribution, 
and  total  threshold  gain  including  diffractive  loss.  This  analysis  should  be  the  most  accurate 
field  calculation  performed  to  date.  From  these  results,  I  infer  the  fundamental  mechanisms 
governing  VCSEL  diffraction. 

1.5  Parasitic  Modes 

Perhaps  the  most  referenced  calculation  addressing  diffraction  in  VCSELs  is  a  model 
based  on  an  “unfolded  cavity”  [29].  In  this  model,  a  Gaussian  mode,  whose  shape  is  deter¬ 
mined  by  the  effective  “fiber  optic”  waveguide  [26],  is  propagated  through  a  series  of  oxide 
apertures.  The  energy  lost  to  diffraction  in  one  cavity  round-trip  is  calculated  and  then 
extrapolated,  based  on  mirror  losses  (e.g.,  emitted  light),  to  estimate  the  threshold  gain. 
This  model  is  popular  for  several  reasons:  (1)  it  yields  an  intuitive  semianalytic  field  so¬ 
lution,  the  VCSEL  analog  of  the  standard  laser  solution,  (2)  it  runs  quickly  and  does  not 
demand  the  computational  horsepower  of  numerical  techniques,  and  most  importantly  (3) 
it  qualitatively  predicts  many  of  the  experimental  observations.  Based  on  this  model,  Santa 
Barbara  derived  their  VCSEL  design:  thin  oxides  placed  at  nodes  of  the  field  standing  wave 
one  or  two  wavelengths  from  the  cavity.  However,  the  current  (repeatable)  threshold  record 
belongs  to  the  Texas  group,  whose  design,  as  discussed  above,  is  diametrically  opposed  to 
Santa  Barbara’s.  The  problem  with  the  unfolded  cavity  model  is  it  fails  to  incorporate 
the  fact  that  the  VCSEL  field  diffracts  into  a  dielectric  layered  structure,  similar  to  a  slab 
waveguide,  and  not  into  free  space  [19]. 
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D.G.  Deppe  (U.T.)  has  proposed  a  different  picture  of  VCSEL  diffraction,  considering 
it  as  energy  loss  to  “parasitic  modes.”  These  are  slab  waveguide  modes,  nearly  resonant 
with  the  lasing  mode,  propagating  in  the  transverse  direction^.  Based  on  this  picture,  his 
group  performed  qualitative  estimates  of  diffractive  loss  in  various  cavity  structures  [20]. 
Their  interpretation  of  the  loss  mechanisms  support  my  conclusions  from  the  FEM  analysis, 
suggesting  the  validity  of  the  parasitic  mode  concept.  Most  remarkably,  both  calculations 
(Texas’s  and  my  FEM)  consistently  explain  the  seemingly  opposite  U.C.S.B.  and  U.T.  op¬ 
timized  designs.  Following  the  Texas  success,  in  Chapter  V,  I  incorporate  parasitic  mode 
loss  into  the  WIM  to  complete  the  optical  model.  Since  the  WIM  is  primarily  intended  as 
a  fast  approximate  technique,  I  improve  the  efficiency  of  the  method — without  sacrificing 
significant  accuracy — by  adopting  the  linear  polarization  (or  paraxial)  approximation.  This 
allows  me  to  treat  the  lasing  mode  (but  not  the  parasitic  modes)  as  a  scalar  field,  greatly  sim¬ 
plifying  the  field  representation  and  cutting  run  times  by  roughly  a  factor  of  two.  I  validate 
this  approximation  through  comparisons  with  previously  obtained  vector  WIM  results. 

1.6  Complete  Laser  Modeling 

The  principal  threshold  metric  is  the  spontaneous  emission  factor  (I3) ,  describing  the 
percentage  of  photons  (spontaneously)  emitted  into  the  lasing  (or  any  other)  mode.  For 
traditional  lasers  with  cavity  dimensions  much  larger  than  the  optical  wavelength,  this  can 
be  accurately  estimated  using  the  three-dimensional  density  of  optical  modes — ^proportional 
to  In  VCSELs,  the  short,  high  Q  cavity  quantizes  the  propagation  constant  in  the  lon¬ 
gitudinal  direction  and  the  appropriate  density  of  modes  is  two-dimensional — proportional 
to  uj.  Furthermore,  the  quantum  well  emission  becomes  anisotropic  in  the  cavity,  increasing 
spontaneous  emission  in  some  directions  at  the  expense  of  others  [8, 10, 18].  The  combina¬ 
tion  of  these  two  phenomena  result  in  lasing  mode  /3^s  much  larger  than  traditional  lasers. 
Extrapolating  this  effect  to  lower  dimensional  systems  (i.e.,  more  degrees  of  confinement), 
the  density  of  modes  may  be  further  reduced  and  larger  j3^s  obtained.  In  the  theoretical, 
zero-dimensional  limit — where  all  cavity  dimensions  are  on  the  order  of  the  wavelength — the 

®More  precisely,  the  parasitic  modes  propagate  at  least  in  the  transverse  direction;  they  may  also  propa¬ 
gate  in  the  longitudinal  direction  as  well. 


1-10 


density  of  modes  becomes  a  series  of  delta  functions.  If  these  modes  are  spaced  sufficiently 
apart,  only  one  mode  may  interact  with  the  gain  media.  In  this  case,  all  the  spontaneous 
emission  will  be  channeled  into  this  mode,  ^  will  go  to  1/2  (due  to  polarization  degeneracy), 
and  the  threshold  will  disappear  [4,5].  In  realistic  VCSELs  this  limit  (currently)  can’t  be 
reached.  Parasitic  optical  and  electrical  losses  destroy  any  quantum  optical  benefits  well 
before  the  theoretical  minimum  cavity  size.  However,  as  evidenced  by  record  low  thresholds, 
oxide  apertured  devices  must  receive  some  quantum  optical  benefits  over  planar  devices. 
Since  the  oxide  confinement  is  incomplete,  one  could  envision  the  optical  density  of  modes 
as  a  set  of  discrete  cavity  modes,  similar  to  the  zero- dimensional  case,  superimposed  on  the 
background  two-dimensional  density  of  modes.  In  Chapter  VI,  I  employ  this  new  approach 
using  my  WIM  optical  modes  to  estimate  the  spontaneous  emission  factor  for  the  University 
of  Texas’s  low  threshold  VCSEL.  Integrating  these  results  with  J.P.  Loehr’s  gain  model  [42], 
the  photon  and  carrier  related  rate  equations  are  solved  to  estimate  the  threshold  current 
and  the  light-vs-current  characteristics — the  ultimate  goal  of  laser  modeling. 
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11.  Vector  Weighted  Index  Method 

In  this  chapter,  I  generalize  the  weighted  index  method  (WIM)— a  separable  approximation 
to  the  Helmholtz  equation — to  compute  cavity  modes  in  cylindrically-symmetric  dielectric 
VCSEL  structures  [50,51].  I  calculate  the  electric  and  magnetic  vector  potentials  and  use 
these  to  compute  the  resulting  fields.  Using  the  calculus  of  variations  in  Appendix  A,  I  show 
that  this  technique  provides  the  best  separable  solution  to  the  scalar  Helmholtz  equation.  The 
method  allows  me  to  approximate  the  spatial  profile,  optical  confinement  factor,  resonant 
frequency,  and  mirror  and  absorption  losses  of  cavity  modes  in  both  oxide-apertured  and 
etched-post  VCSELs.  The  method  explicitly  considers  complex  media,  allowing  me  to  include 
free  carrier  losses.  A  central  flaw  of  the  technique,  however,  is  an  exclusion  of  diffraction,  the 
dominant  loss  mechanism  in  recent  low  threshold  designs.  I  remedy  this  later  in  Chapter  V, 
by  perturbatively  including  these  losses  in  the  general  WIM  formalism  as  presented  below. 


2.1  Vector  Field  Equations 

I  want  to  find  the  electric  and  magnetic  field  profiles,  the  resonant  wave¬ 
length  (A),  and  the  mirror  and  absorption  losses  _j.  ^(absorption)^  £qj,  cavity  mode 

in  azimuthally-symmetric  VCSEL  structures.  For  this  I  must  solve  a  vector-wave  equation 
subject  to  appropriate  boundary  conditions  at  each  interface.  Because  there  are  several 
equivalent  electromagnetic  descriptions  of  any  system,  I  can  write  wave  equations  for  the 
electric  and  magnetic  fields,  scalar  potentials,  or  vector  potentials.  The  most  powerful  and 
convenient  method  for  this  problem  is  to  solve  for  the  magnetic  and  electric  vec¬ 
tor  potentials  and  use  them  to  compute  the  fields.  The  steady-state,  time-harmonic  vector 
potentials  A  and  F  satisfy  the  three-dimensional  vector  Helmholtz  equation  (in  MKS  units) 


V^  +  ^er{p,z) 


A(p,  <i>,  z) 
F{p,  (f),  z) 
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Here  A  and  F  depend  on  time  as  e*"*,  u  =  27rc/A,  and  I  have  assumed  a  cylindrically- 
symmetric,  complex  dielectric  function  .  Note  that 


y/Fi-  =  ?;  =  n  +  i/c, 


(2.2) 


where  t]  is  the  (complex)  refractive  index;  material  gain  is  incorporated  by  taking  «  positive 
in  the  active  region.  I  assume  =  1  in  all  regions. 

For  azimuthally  symmetric  structures,  the  vector  potentials  are  separable  in  (f),  de¬ 
pending  on  it  as  Furthermore,  expressing  the  potentials  in  cylindrical  coordinates,  the 

z  components  obey  a  scalar  Helmholtz  equation. 


^ 

dp^  ^  p  dp  ^  dz"^ 


t  \ 


(2.3) 


To  generate  transverse  magnetic  (to  z)  (TM)  modes  I  must  solve  for  while  to  generate 
transverse  electric  (TM)  modes  I  must  solve  for  [28].  Since  an  arbitrary  mode  can  be 
represented  as  a  superposition  of  TE  and  TM  modes,  I  need  only  solve  for  the  two  unknown 
scalar  functions  Az  and  Fz,  a  dramatic  simplification  over  solving  (2.1)  for  all  six  vector 
components.  Note  that  despite  the  fact  that  I  may  solve  a  scalar  equation  (2.3)  to  find 
modes,  the  vector  nature  of  the  solution  remains — manifested  in  the  boundary  conditions. 
The  power  of  the  vector  potential  approach  comes  from  the  fact  that  I  may  generate  all  the 
vector  field  components  by  taking  derivatives  of  the  z  components  of  the  vector  potentials, 
therefore  requiring  no  a-priori  knowledge  of  the  p  ox  z  dependence.  The  fields  themselves 
are  related  to  the  vector  potentials  by  [36] 


E 

H 


- X  V  X  (zAz)  -  V  X  (zFz), 

LOSaSr 

V  X  (zAz) - X  V  X  {zFz), 

ojpq 


(2.4) 

(2.5) 


2-2 


or,  more  explicitly, 


Ep{p,z)  = 
E^{p,z)  = 
Ez{p,z)  = 
Hp{p,z)  = 
H<t>{p,^)  = 

Hz{p,z)  = 


i  f  Id  .  /  X 

-  1  'a~2  — "fl  r  I  ■^)’ 

weoer  I dp^  pap  p^  } 


d'^ 


E i{pt  z), 


i_^ .  /  X _ i _ 

pd(j>  upodpdz 

-—A,{p,z)  - 

i  (  d"^  1  d  m?\  .  . 

-  I  'a~2  — « - T  I  Ezipj  z). 

upo  { dp^  pdp  p2  j 


(2.6) 

(2.7) 

(2.8) 

(2.9) 

(2.10) 

(2.11) 


In  separable  geometries,  the  two-dimensional,  azimuthally-symmetric  scalar  Helmholtz 
equation  (2.3)  may  be  solved  exactly  by  separation  of  variables,  yielding  the  potential  pro¬ 
files,  resonant  wavelength,  and  mirror  and  absorption  losses  for  each  cavity  mode.  Realistic 
VCSEL  structures,  however,  are  not  separable,  greatly  complicating  the  solution  of  (2.3). 
An  exact  semianalytic  solution  could  be  obtained  by  expanding  Az  and  Fz  in  terms  of  the 
general  solutions  in  each  region,  then  matching  boundary  conditions  to  determine  the  (infi¬ 
nite)  set  of  expansion  coefficients.  In  practice,  this  technique  requires  considerable  care  to 
implement,  though  it  does  have  the  advantage  of  incorporating  non-separable  behavior  in 
the  solutions  [12].  Below,  I  present  an  alternative  technique  to  generate  the  best  separable 
approximations  to  (2.3). 


2.2  The  Weighted  Index  Method 

Equation  (2.3)  represents  two  uncoupled  partial  differential  equations — one  each  for  Az 
and  Fz — which  are  quite  difficult  to  solve.  For  separable  geometries,  I  can  exactly  replace 
each  equation  in  (2.3)  with  two  independent  ordinary  differential  equations,  and  these  can 
be  solved  exactly.  For  non-separable  geometries,  I  approximate  the  solutions  to  (2.3).  In 
general,  there  are  two  possible  approximation  techniques.  The  most  common  approach  is 
to  maintain  the  exact  equations  (2.3)  and  construct  an  approximate  function  that  “almost” 
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solves  them.  An  alternative  approach  is  to  replace  the  exact  equations  (2.3)  with  approx¬ 
imate  equations,  and  solve  these  approximate  equation  exactly.  I  take  the  latter  approach, 
and  approximate  each  equation  in  (2.3)  with  two  coupled  ordinary  differential  equations.  I 
accomplish  this  by  extending  the  weighted  index  method  (WIM) — which  was  first  developed 
to  calculate  waveguide  modes  in  horizontal-cavity  ridge- waveguide  lasers  [39,61] — to  address 
the  eigenmodes  of  cylindrical  cavities.  This  technique  has  the  advantage  of  giving  the  best 
separable  solution  to  (2.3)  in  the  variational  sense,  and  allows  me  to  estimate  the  field  profile, 
optical  confinement  factor,  resonant  wavelength,  and  mirror  and  absorption  losses  for  each 
cavity  mode.  Below  I  derive  the  WIM  equations;  in  Appendix  A,  I  prove  that  their  solution 
yields  the  best  separable  approximation  to  (2.3). 

Proceeding  as  if  separable  solutions  to  (2.3)  exist,  I  take 


A^{p,z)  =  P{p)Q{z), 


F,{p,z)  =  R{p)S{z). 


(2.12) 


Substituting  either  of  these  into  (2.3)  gives 

C"W«2)  +  -  ^)  C(/-)«^)  =  0,  (2.13) 

where  (  =  P  or  R  and  ^  =  Q  or  S.  For  each  potential  in  (2.12),  I  can  separate  the 
resulting  equation  (2.13)  by  integrating  it  against  C*(/^)  or  ^*(;?).  This  procedure  yields  the 
axial  equation 

=  (2.14) 


and  the  radial  equation 


for  each  potential.  These  axial  and  radial  equations  are  coupled  by  the  weighted  axial  and 
radial  propagation  constants,  given  respectively  by 


,  ,-,2 _  (u.Vc^)(C|e.(o,^)K)  +  (CK")  +  (Clp-*IC')-m^(Clo-^K) 

Wz/W)  =  (<!(;) 


(2.16) 


and 


{k,mY  = 


UH) 


Here  (  |  )  denotes  an  inner-product  over  2;  or  p,  respectively  defined  by 


(2.17) 


/CO 

A*{z)B{z)  dz, 

•00 


(2.18) 


and 

poo 

( A(p)  \B(p))=  A\p)B(p)p dp.  (2.19) 

Jo 

Since  the  weighted  propagation  constants  depend  on  whether  I  solve  for  or  Fz,  I  have 
introduced  an  additional  superscript  a  =  TE  or  TM  to  distinguish  between  TE  modes 
resulting  from  Fz  (involving  averages  over  (  R  and  ^  =  S)  and  TM  modes  resulting  from 
Az  (involving  averages  over  C  =  P  and  C  =  Q).  Coupling  occurs  only  between  the  radial  and 
axial  equations  for  a  given  vector  potential  Az  or  Fg-.  the  two  potentials  remain  uncoupled 
in  (2.14)  and  (2.15).  But  Az  and  Fz  will  be  coupled  later  by  boundary  conditions  when  I 
solve  for  hybrid  modes. 


2.3  Weighted  Boundary  Conditions  and  Solutions 

In  this  section  I  discuss  the  solutions  of  (2.14)  and  (2.15)  in  piecewise-constant  re¬ 
fractive  index  profiles,  paying  particular  attention  to  the  interfacial  boundary  conditions. 
A  sample  structure  is  shown  in  Figure  2.1.  Since  I  will  work  exclusively  with  piecewise- 
constant  geometries,  I  simplify  my  notation  by  taking  er{p,z)  Sr-.i,j,  where  i  and  j  index 
the  radial  and  axial  regions,  respectively.  Thus  I  also  have 

and  (2.14)  and  (2.15)  reduce,  respectively,  to  the  one-dimensional  Helmholtz  equation  and 
Bessel’s  equation. 
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Figure  2.1  Illustration  of  the  piecewise  constant  permittivity  notation  for  an  axially  sym¬ 
metric  VCSEL. 
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To  solve  these  equations  subject  to  the  refractive  index  profile  £r:i,j  of  etched-post  or 
oxide-apertured  VCSEL  structures,  I  must  supplement  (2.14)  and  (2.15)  with  an  appropriate 
set  of  interface  and  endpoint  boundary  conditions.  The  interfacial  boundary  conditions 
are  the  usual  continuity  requirements  on  the  normal  and  tangential  components  of  various 
electromagnetic  fields.  Therefore  I  cannot  directly  enforce  boundary  conditions  on  P,  Q, 
P,  and  S,  but  must  perform  the  intermediate  step  of  computing  the  electric  and  magnetic 
fields  via  (2.6)-(2.11).  Furthermore,  since  the  underlying  partial  differential  equations  (2.3) 
do  not  separate,  these  boundary  conditions  cannot  be  satisfied  at  all  points  on  the  boundary 
surfaces — if  they  could,  (2.3)  would  be  separable. 

In  preparation  for  generating  approximate  boundary  conditions,  I  rewrite  the  weighted 
index  formulas  (2.16)  and  (2.17).  By  integrating  (2.15)  against  C*{p)  I  have 

(CIC">  +  (CI/>-|C')-'n'(C|p-"l<)  =  -{CIWflf).  (2.20) 


allowing  me  to  express  (2.16)  as 

(;?“)’  =  (‘^V<^)(ci^r,.rfic)-(ci(i^r)^!C)  ^  (2.21) 


where 


/(CKWIC)  (ooo. 

=  (flC)  “d  <*^  >  =  V  (CIO  • 


Similarly,  by  integrating  (2.14)  against  C(z)  I  find 


(2.23) 


allowing  me  to  express  (2.17)  as 


,  _  (u;Vc^)(f  |e.ulll)-(iM)llO 


UJ 


(eio 


(2.24) 


where 


I  a  V  _  ( ^  1  ^r:i,j  \  ^  ^  ,ga\  _  (  ^  I  (/^j  )^  I  0 

^^1^)  and(^)_W  , 


(2.25) 
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The  compact  expressions  (2.21)  and  (2.24)  allow  me  to  compute  weighted  variables  without 
using  the  derivatives  of  C  and 

I  specialize  these  expressions  to  VCSEL  lasing  modes  by  truncating  the  inner  product 
over  z,  replacing  (2.18)  with 

pZtop 

{ A{z)  I  B(z) )  =  /  A\z)B{z)  dz,  (2.26) 

^-^^bottom 

where  ^bottom  and  Ztop  denote  the  lower  and  upper  VCSEL  boundary  planes.  This  truncation 
is  necessary  to  force  (2.24)  to  converge  and  reflects  the  assumption  that  most  of  the  energy 
is  contained  inside  the  VCSEL  cavity.  In  contrast,  I  force  the  radial  wavefunctions  to  decay 
evanescently  to  zero — to  find  guided  modes — and  the  inner  product  defined  by  (2.19)  presents 
no  difficulty.  It  is  essential  to  realize  that  by  forcing  evanescent  behavior  in  the  radial 
direction,  I  explicitly  prohibit  diffraction  from  contributing  to  the  solution.  This  unavoidable 
drawback  is  repaired  perturbatively  in  Chapter  V. 

I  now  present  boundary  conditions  and  solutions  for  the  axial  and  radial  equations 
(2.14)  and  (2.15).  The  unknowns  in  this  formalism  are,  for  each  mode,  the  functions  P{p), 
Q(z),  R{p), and  S{z),  the  resonant  frequency  w,  and  the  mirror  and  absorption  losses, 

^(mirror)  ^(absorption)  =  r*°*47r/Cactive/A.  (2.27) 

Here  T*®*  is  the  total  confinement  factor  and  /Cactive  is  the  material  gain  in  the  active  quantum 
well  region  required  to  offset  the  mirror  and  absorption  losses.  Note  that  the  sum  of  the 
two  losses  equals  the  threshold  modal  gain  which  would  be  observed  if  diffraction  were  not 
present.  In  practice,  I  explicitly  solve  for  /Cactive?  not  +  o;(absorption)_ 

2.S,1  Axial  Boundary  Conditions  and  Solutions.  I  solve  the  axial  equation  (2.14) 
in  piecewise  constant  geometries,  such  as  in  Figure  2.2.  The  general  solutions  of  (2.14)  are 
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given  in  each  axial  region  [zj,Zj^i]  by 


Qj{z)  = 

(2.28) 

where  I  explicitly  denote  both  the  TE  and  TM  solutions  for  clarity.  Using  the  iterative 
solution  procedure  described  in  Section  2.4,  I  compute  and  from  (2.21);  assume 
for  now  that  they  are  known  constants.  These  solutions  must  be  joined  at  each  interface  Zj 
by  matching  the  tangential  electric  and  magnetic  fields.  Inserting  (2.12)  into  (2.6)-(2.11),  I 
compute  these  tangential  field  components  as 


Ep{p,z)  = 

- —PWiz)  -  -Rip)S{z), 

(2.29) 

UJeoSr  P 

E^{p,z)  = 

+  b:(p)s(z), 

(2.30) 

LOe^Srp 

Hp{p,z)  = 

^-^P{p)Qiz)  -  ^R\p)S\z), 

p  upq 

(2.31) 

H4,{p,z)  = 

-P'MQU)  + 

uPqP 

(2.32) 

I  consider,  in  turn,  two  distinct  cases:  m  =  0  and  m  ^  0. 

When  m  =  0  I  can  match  boundary  conditions  with  Q  =  0  (pure  TE  modes)  or 
5  =  0  (pure  TM  modes).  For  these  modes,  it  is  sufficient  to  force  just  two  tangential 
field  components  to  be  continuous:  demanding  continuity  of  the  other  component  gives  a 
redundant  condition.  For  TE  modes  I  require  and  Hp  to  be  continuous;  for  TM  modes 
I  require  Ep  and  to  be  continuous.  As  I  will  show  below,  for  pure  TE  and  TM  modes 
the  fields  themselves,  and  not  just  the  vector  potentials,  are  separable.  This  makes  it  easy 
to  generate  weighted  boundary  conditions  for  these  modes. 

For  pure  TE  modes,  I  have 

E^{p,z)  =  R'{p)S{z), 

Hp{p,z)  =  -^R\p)S'{z). 
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(2.33) 

(2.34) 


Figure  2.2  Illustration  of  the  weighted  permittivity  profile,  and  notation,  for  the  WIM 
axial  solution. 
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Since  I  assume  /[^^  =  1  in  all  regions,  both  and  Hp  depend  on  z  only  via  S[z)  and 
respectively.  Therefore  I  can  force  both  tangential  fields  to  be  continuous  by  setting  S{z) 
and  S'(z)  continuous  across  each  interface. 

For  pure  TM  modes,  on  the  other  hand,  I  have 

Ep{p,z)  = - i-F{p)Q'{z),  (2.35) 

H^{p,z)  =  P'{p)Q{z).  (2.36) 

Again,  I  can  make  continuous  by  forcing  Q{z)  to  be  continuous  across  each  interface. 
However,  Ep  depends  on  z  through  both  Q'(z)  and  e^.  Since  er{p,z)  is  not  a  separable 
function,  I  must  weight  the  relative  permittivity  to  obtain  average  boundary  conditions 
holding  for  all  p.  Therefore  I  require  Q\z)l{eJ^)  to  be  continuous  at  each  interface,  where 
(e™j)  has  already  been  defined  in  each  region  by  (2.25).  (I  could  also  have  generated  average 
boundary  conditions  using  (^),  but  this  approach  gave  inferior  results.) 

Inserting  the  functional  forms  (2.28)  and  applying  the  continuity  conditions  for  either 
a  =  TE  or  a  =  TM  modes,  I  relate  aj,  bf  to  aj+i,  at  each  axial  boundary  z  =  Zj 
through  the  transfer  matrices 


(2.37) 
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The  composite  transfer  matrix  [T"]  for  the  whole  system  is  formed  by  cascading  the  indi¬ 
vidual  transfer  matrices,  giving 


\T-]  =  w-ii- 


(2.40) 


Here  N  is  the  number  of  axial  regions  (including  substrate  and  air)  in  the  problem  geometry. 
Thus  I  relate  the  unknown  coefficients  in  the  j  =  \  region  (substrate)  to  the  coefficients  in 
the  j  =  N  region  (air)  via 


a? 

6“ 


ni  *12 

4-01 

t'21  ^22 


Oat 


(2.41) 


Since  I  am  searching  for  axially-emitting  (lasing)  modes,  I  permit  only  outgoing  radiation 
by  setting 

6“  =  =  0.  (2.42) 

Finally,  substituting  (2.42)  into  (2.41)  I  obtain  the  axial  threshold  condition^ 


^22(^5  ^active)  —  0- 


(2.43) 


Setting  the  real  and  imaginary  parts  of  zero  gives  two  independent  equations 

that  I  solve  to  obtain  the  modal  frequency  (w  =  27rc/A) ,  and  the  mirror  and  absorption 
losses  (q;(”“''°*')  -1-  Q;(absoipt>o“)  =  r*°*47r«active/A)  for  pure  TE  and  TM  modes.  The  expansion 
coefficients  a"  and  6"  for  each  region  are  found  by  back  substitution  through  (2.41)  and 
(2.37). 

In  order  to  generate  sensible  boundary  conditions  when  m  ^  0,  I  must  construct 
hybrid  modes  in  which  both  S  and  Q  axe  nonzero.  In  this  case  none  of  the  tangential  fields 
are  separable,  since  each  has  both  TE  and  TM  parts.  Each  TE  and  TM  part  is,  in  turn, 
a  sum  of  cylindrical  wave  terms  like  e*‘^'^Jm{kp)e*^^  (P  and  R  will  turn  out  to  be  Bessel 
functions).  An  exact  solution  would  require  me  to  include  a  superposition  of  cylindrical 

*As  shown  in  Appendix  B,  this  requirement  in  conjunction  with  the  treinsfer  matrix  solution  is  a  gener¬ 
alization  of  the  classic  Fabry-Perot  laser  round-trip  amplitude  and  phase  conditions. 
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waves  involving  all  values  of  k,  and  the  boundary  conditions  would  couple  all  terms  at  each 
interface.  Fortunately,  the  coupling  between  cylindrical  waves  with  different  k  is  small  enough 
to  ignore  [12,34,48].  Therefore  I  simultaneously  and  independently  enforce  the  continuity  of 
the  dominant  TE  and  TM  mode  components  of  each  tangential  field,  and  my  hybrid  mode 
boundary  conditions  become  the  same  as  those  for  pure  TE  and  TM  modes.  Although  these 
boundary  conditions  are  approximate,  I  feel  they  are  reasonable  within  the  spirit  of  the 
method.  Using  more  rigorous  boundary  conditions  would  require  a  very  complex  procedure 
of  questionable  value  given  the  fundamental  approximate  nature  of  the  solution. 


2.3.2  Radial  Boundary  Conditions  and  Solutions.  I  solve  the  radial  equation  (2.15) 
in  piecewise  constant  geometries,  such  as  in  Figure  2.3.  The  general  solutions  of  (2.15)  are 
given  in  each  radial  region  [ft,  pj+i)  by 


c™Jmik™p)  +  d[^Ym{kJ^p)  i  /  M 
cJ^Km(ikJ^p)  +  ^Imiik™p)  i  =  M 


Hp)  =  i 


CM^miikJpp)  +  <i^lm{ikjpp) 


i^M 
i  =  M 


(2.44) 


Here  J^.  and  Ym  are  m-th  order  Bessel  functions  of  the  first  and  second  kind,  Im  and  Km 
are  modified  m-th  order  Bessel  functions  of  the  first  and  second  kind,  and  i  =  1,2, ...  ,M 
indexes  the  inner  to  outer  radial  regions.  Using  the  iterative  solution  procedure  described 
in  section  2.4,  I  compute  kj^  and  k™  from  (2.24);  assume  for  now  that  they  are  known 
constants.  These  solutions  must  be  joined  at  each  interface  ft  by  matching  tangential  elec¬ 
tric  and  magnetic  fields.  Inserting  (2.12)  into  (2.6)-(2.11)  I  compute  these  tangential  field 
components  as 


Ez{p,z)  = 

{k™y  p{p)Q{z), 

(2.45) 

E4>iP,z)  = 

-^P(p)QXz)  R'{p)S(z), 

(2.46) 

H.{p,z)  = 

{kj^)  Rip)S(z), 

(2.47) 
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(2.48) 


777 

h4p,z)  =  -np)Q{^)+—Mi>)sV)- 

OJfMQP 

I  consider,  in  turn,  the  distinct  cases  m  =  0  and  m  ^  0. 

As  in  the  axial  problem,  when  m  =  0  I  can  match  boundary  conditions  with  P  =  0 
(pure  TE  modes)  or  =  0  (pure  TM  modes),  and  it  is  sufficient  to  force  just  two  tangential 
field  components  to  be  continuous.  For  TE  modes  I  require  and  Hz  to  be  continuous;  for 
TM  modes  I  require  Eg  and  to  be  continuous.  Following  the  same  arguments  as  in  the 
axial  problem,  I  construct  weighted  boundary  conditions  to  require  the  continuity  of 

R{p)  and  R'{p)  (for  TE  modes),  (2.49) 

/'jLTM\2 

/  ™  niodes),  (2.50) 

\^r:Q  ) 

at  each  radial  interface  p  =  p,-.  To  generate  radial  boundary  conditions  independent  of  I 
have  replaced  by  the  weighted  permittivity  {s™q)  defined  in  (2.25). 

Inserting  the  functional  forms  (2.44)  and  applying  the  continuity  conditions  for  either 
a  =  TE  or  TM  modes,  I  relate  cf,  df  to  cf+j,  at  each  radial  boundary  p  =  p,-  through 
the  transfer  matrices 

where 

aP 
BP 


bP 
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c" 

d’t 


Bf 


‘'<+1 

®t+l 


(2.51) 


krJL(l=7^Pi)  kJ^YiikT^p,) 

(kWiT 

■yk^-Lhi^WiPi) 

kTMj;^(kT^^,)  kT^Y^ik^P,) 
0LE„(jk™p,) 


(2.52) 

(2.53) 

(2.54) 

(2.55) 


Here  I  have  defined 


1  inner  regions 

^  5  ^ 

Jm  inner  regions 

1  ^  j 

^  ?  Ilm  —  ^ 

Ym  inner  regions 

1 

y/^  outer  region  J 

1  Km  outer  region  J 

Im  outer  region 

/n 

(2.56) 


As  in  (2.40),  I  form  a  composite  system  matrix 

[(7»]  =  W1  ■  ■  ■  (2.57) 


relating  the  unknown  coefficients  in  the  z  =  1  region  (core)  to  the  coefficients  in  the  i  M 
region  (cladding)  via 

c? 

Since  I  am  searching  for  longitudinally-propagating,  laterally-confined  VCSEL  modes  I  force 
regularity  at  the  origin  and  exponential  decay  as  p  ^  oo  by  setting 

d-  =  =  0.  (2.59) 


u 


11 


u 


u 


12 


21 


U 


22 


“M 


(2.58) 


Finally,  substituting  (2.59)  into  (2.58)  I  obtain  the  radial  threshold  condition 


^21(^5  ^2ictive)  —  0. 


(2.60) 


In  order  to  generate  sensible  boundary  conditions  when  m  7^  0, 1  must  construct  hybrid 
modes  in  which  both  R  and  P  are  nonzero.  Again,  none  of  the  tangential  field  components 
are  separable  and  it  is  impossible  to  match  them  for  all  2;  at  a  radial  interface.  Furthermore, 
if  I  mirror  my  axial  treatment  and  independently  force  the  TE  and  TM  components  to  be 
continuous  I  generate  inconsistent  boundary  conditions.  For  example,  independently  forcing 
Ez  and  the  TM  part  of  to  be  continuous  requires  [k™Y  T’(p)/(e™)  and  P{p)/ (e™)  be 
continuous  at  each  interface.  But  this  forces  k™  itself  to  be  continuous,  which  is  absurd.  I 
cannot  hope  to  work  with  TE  and  TM  modes  independently.  Therefore,  I  approximate  my 
problem  with  an  equivalent  cylindrical  dielectric  waveguide  problem,  which  admits  analytic 
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solutions  for  hybrid  modes.  These  approximations  axe  used  only  to  generate  boundary 
conditions,  not  as  a  substitution  for  the  actual  radial  and  axial  solutions  in  each  region. 

Cylindrical  waveguide  modes  depend  on  z  as  or  My  fields,  in  contrast,  depend 
on  z  via  (2.28).  There  are  two  differences  I  must  overcome.  First,  cylindrical  waveguide 
modes  are  characterized  by  a  single  axial  propagation  constant  jS,  whereas  my  fields  have  a 
different  axial  propagation  constant  /3“  for  each  axial  region  and  polarization.  I  can  easily 
remedy  this  by  replacing  with  (/?“),  as  defined  in  (2.25).  Second,  cylindrical  waveguide 
modes  depend  on  z  as  either  or  But  even  after  replacing  /?“  with  (^“)  my  fields 

have  a  different  linear  combination  of  and  in  each  region,  depending  on  the 

relative  values  of  a"  and  6".  Therefore,  for  the  purpose  of  constructing  boundary  conditions, 
I  assume  that  the  fields  approximate  a  “pure”  standing  wave  in  the  axial  direction,  with 


uf  =  and  =  -bj^. 


(2.61) 


I  assume  further  that  aj®  «  aj^,  and  approximate  the  z-dependence  of  my  fields  as 

Q{z)  =  sin((^™)z), 

S(z)  =  cos((i3'^^}z). 

Substituting  (2.62)  into  (2.46)  and  (2.48),  I  approximate  the  ^  field  components  as 


(2.62) 


B4p,z)  = 


m 


P(p)Q'(z)  +  B'(p)S(z), 


^^0(£r:Q}P 


(2.63) 


and 


m 


ff4p,z)  =  -P’(p)Q(z)  + —X(f,)S'(z), 

COPoP 


I  ^PoP  J 


(2.64) 
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Here  I  have  replaced  e^ij  by  the  appropriate  weighted  value,  and  have  selectively  assumed 
«  (^™)  «  (^)  to  factor  out  the  z  dependence.  Through  a-posteriori  comparison  with 
my  calculated  results,  I  find  all  these  assumptions  well  justified. 

These  approximate  expressions  for  and  are  separable,  as  are  expressions  (2.45) 
and  (2.47)  for  and  H^.  Therefore  I  can  immediately  construct  suitable  boundary  condi¬ 
tions  for  radial  hybrid  modes  by  requiring  the  continuity  of 


/„TM\ 

\^r:Q  1 

(2.65) 

R{p), 

(2.66) 

(2.67) 

(2.68) 

By  inserting  the  functional  forms  (2.44)  and  forcing  all  four  of  the  above  combinations  to  be 
continuous,  I  link  the  unknown  coefficients  at  each  radial  interface  p  =  pi  through  the  4  x  4 
transfer  matrix 

A.- 

where 
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(kTf.fSrr.i'rkTMp.)  (kj:^,fn^k-^,Pi)  Q  0 

'7°  'o’  (kW^f^ill^Ihl’i)  (*:5?.)'n„(7*:^®,A) 

"" 

'ik^.^uhki^d  7is^'„(7^5.■^/-.)  J 

>■  (2.71) 

These  matrices  are  mathematically  equivalent  to  those  for  a  cylindrical  dielectric  waveguide. 
Cascading  the  interface  transfer  matrices,  I  again  derive  a  composite  system  transfer  matrix 


Un  Ui2  Wl3  «14 

U21  ^22  U23  W24 

*^31  ^32  ^33  ^^34 

U41  U42  U43  ^44 


(2.72) 


relating  the  innermost  and  outermost  radial  coefficients.  The  endpoint  boundary  conditions 
(2.59)  remain  valid.  Applying  them  to  (2.72)  and  demanding  nontrivial  solutions  gives  the 
hybrid  threshold  condition 


U2l(w,  Kactive)  ^23(^5  ^active) 
^41  (^5  ^active)  ^43(^5 /Cactive) 


(2.73) 


which  I  solve  in  the  complex  plane  for  uj  and  /Cactive- 


2.4  Iterative  Solution  Procedure 

I  compute  the  longitudinal  and  transverse  mode  spectrum  by  self-consistently  solving 
the  radial  and  axial  problems.  The  modes  are  specified  by  the  longitudinal  mode  number, 
the  transverse  mode  number,  and  the  azimuthal  mode  number  m.  Different  longitudinal 
modes  correspond  to  successive  roots  of  the  axial  threshold  condition  (2.43),  while  different 
transverse  modes  correspond  to  successive  roots  of  the  radial  threshold  conditions  (2.60) 
(for  TE/TM  modes)  or  (2.73)  (for  hybrid  modes).  The  energy  spacing  between  longitudinal 
modes  is  much  greater  than  that  between  trajisverse  modes.  If  I  let  v  denote  a  generalized 
mode  index  corresponding  to  a  particular  TE,  TM,  or  hybrid  mode,  then  the  explicit  un- 
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knowns  for  each  mode  are  the  vector  potential  functions  P^(/)),  and/or  5^(2;), 

the  optical  mode  frequency  a;^,  and  the  material  gain  values  (Kactive)^-  I  iteratively  solve  the 
axial  and  radial  problems  as  follows. 

First  I  focus  attention  on  a  particular  family  of  modes  by  fixing  m:  for  m  =  0  I  can 
compute  TE  or  TM  modes,  while  for  m  >  0  only  hybrid  modes  are  possible.  Then  I  solve 
the  axial  problem,  taking  the  effective  indices  equal  to  the  corresponding  values  for  in 
the  innermost  radial  region:  this  solution  corresponds  to  the  standard  plane- wave  calculation 
appropriate  for  large-area  devices.  The  ordered  roots  of  the  axial  threshold  condition  (2.43) 
yield  initial  approximations  Quoi^)  and/or  and  ( mode  i/q  of 

interest.  For  VCSEL  lasing  mode  calculations  I  am  only  interested  in  the  first  root  of  (2.43), 
corresponding  to  the  fundamental  longitudinal  mode.  I  next  compute  the  kf  by  inserting 
Quoi^)  and/or  S^^{z)  into  (2.24),  giving  me  enough  information  to  address  the  radial  problem. 
Depending  on  whether  I  am  solving  for  TE,  TM,  or  hybrid  modes,  I  then  compute  the  roots 
of  the  radial  threshold  condition  (2.60)  or  (2.73),  yielding  approximations  for  Pp^ip)  and/or 
Rpoip),  and  ;  successive  roots  correspond  to  progressively  higher-order 

transverse  modes.  Then  I  alternate  between  solving  the  axial  and  radial  problems,  always 
updating  ,0"  and  kf  by  inserting  the  most  recent  wavefunctions  into  (2.21)  or  (2.24)  as 
appropriate.  In  this  way  I  generate  a  self-consistent  solution  to  the  coupled  WIM  equations 
(2.14)  and  (2.15),  terminating  when  and 

within  a  prescribed  tolerance.  The  procedure  converges  quite  rapidly,  allowing  me  to  solve 
for  a  large  number  of  cavity  modes. 

Finally,  I  note  that  the  character  of  the  modes  found  depends  on  both  the  differential 
equations  and  on  the  endpoint  boundary  conditions  enforced.  The  original  application  of 
the  WIM  to  rectangular  waveguide  geometries  assumed  propagating  behavior  in  the  z  di¬ 
rection  and  evanescent  decay  in  the  x  and  y  directions.  These  endpoint  conditions  resulted 
in  a  standard  eigenvalue  problem,  with  the  longitudinal  propagation  constant  ^  being  the 
eigenvalue.  In  this  case,  the  Rayleigh-Ritz  variational  principle  asserts  that  the  resulting 
approximation  of  /3  will  be  more  accurate  than  the  wavefunctions  themselves.  In  my  ap¬ 
plication,  I  have  enforced  evanescent  decay  in  the  radial  direction  and  have  permitted  only 
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outward  propagating  waves  in  the  z  direction.  This  constrains  the  mode  frequency  u)  and  the 
material  gain  k,  as  opposed  to  the  propagation  constant  ji,  and  these  unknown  parameters 
no  longer  appear  as  eigenvalues.  Therefore  the  Rayleigh-Ritz  principle  has  nothing  to  say 
about  the  relative  accuracy  of  a;  and  k.  Nevertheless,  general  variational  principles  dictate 
that  I  have  found  the  best  separable  solution,  and  therefore  I  expect  the  resulting  values  for 
w  and  K  to  be  reasonably  accurate.  Note  that  I  can  also  apply  the  general  WIM  formalism 
to  solve  for  radially  propagating  modes  by  relaxing  the  evanescent  decay  condition  on  the 
radial  functions.  As  it  turns  out,  these  “parasitic”  modes  are  precisely  what  is  required  to 
capture  diffraction  within  the  general  framework  of  the  weighted  index  method.  I  append 
the  current  WIM  treatment  to  include  these  modes  in  Chapter  V.  However,  before  this 
I  proceed  and  apply  the  current  version  of  the  method  to  several  popular  VCSELs  in  the 
following  chapter. 
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III.  Weighted  Index  Method  Applications  and  Results 

In  Chapter  II,  I  presented  a  new,  extremely  general  and  rapid  technique  for  estimating  the 
spatial  profile,  optical  confinement  factor,  resonant  frequency,  and  mirror  and  absorption 
losses  of  lasing  modes  in  cylindrically-symmetric  VCSEL  geometries.  In  this  chapter,  I 
apply  the  weighted  index  method  to  several  VCSEL  geometries  of  current  interest.  The 
chapter  is  broken  into  two  major  parts.  First,  I  compute  all  the  optical  laser  parameters, 
save  diffraction  loss,  for  several  lowest  order  modes  of  an  oxide  apertured  and  an  etched  post 
device.  I  compute  all  these  parameters  as  a  function  of  oxide  aperture  or  post  radius.  These 
general  results  showcase  the  capabilities  of  the  method.  Furthermore,  from  these  results  I 
ascertain  the  basic  changes  in  the  optical  mode  for  small  aperture  geometries.  These  include 
a  blueshift  in  mode  resonance,  a  loss  of  transverse  confinement,  and  a  heretofore  unreported 
polarization  dependent  change  in  DBR  reflectance. 

Second,  I  pit  the  weighted  index  method  against  a  simpler,  more  popular  effective  index 
method  for  blueshift  calculation.  Both  models  are  applied  to  two  different  VCSEL  designs 
emitting  near  780  nm  and  850  nm.  I  find  that  the  weighted  index  method  matches  the  data 
remarkably  well,  while  the  effective  index  calculation  underestimates  the  blueshift.  These 
results  show  that  the  weighted  index  model  may  be  used  to  match  the  cavity  resonance  of 
small-aperture  structures  to  the  gain  peak,  and  thereby  minimize  the  threshold  current. 

3.1  Etched  Post  and  Oxide  Apertured  VCSELs 

Using  the  methods  described  in  the  last  chapter,  I  calculate  lasing  modes  in  etched-post 
[6,76]  and  oxide-apertured  [15,32]  devices  fabricated  from  a  A  =980  nm,  1.5A-cavity  VCSEL. 
The  VCSEL  reflectors  consist  of  a  17.5  period  p-type  GaAs/Alo.92Gao.o8As  top  Distributed 
Bragg  Reflector  (DBR)  for  the  oxide-apertured  structure,  and  a  4  period  GaAs/AUO  top 
DBR  with  a  A/4  p"*"  GaAs  contact  layer  for  the  post  geometry.  The  bottom  reflector  in 
both  structures  is  a  22  period  n-type  GaAs/Alo.92Gao.o8As  DBR.  The  1.5A-cavity  is  step- 
tapered  with  (intrinsic)  layers  of  Al0.98Ga0.02As,  Alo.65Gao.35As,  Alo.30Gao.70As,  and  GaAs, 
culminating  in  a  single  Ino.2Gao.8As  quantum  well.  Both  structures  are  grown  on  a  GaAs 
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substrate.  The  etched-post  structure,  illustrated  in  Figure  3.1,  is  formed  by  etching  the 
top  GaAs/AlAs  DBR  down  to  the  A/4  GaAs  contact  layer,  then  oxidizing  the  AlAs  layers. 
The  oxide-apertured  structure,  illustrated  in  Figure  3.2,  is  formed  by  oxidizing  the  A/4 
Al.98Ga.02As  layers  in  the  cavity.  I  model  each  region  as  a  cylindrically-symmetric  layer  of 
constant  refractive  index,  assuming  the  material  parameters  in  Table  3.1.  Free  carrier  losses 
are  incorporated  by  taking  negative.  Material  gain  is  incorporated  by  taking  /Cactive 
positive  in  the  active  quantum  well  region. 


Material 

Index 

Doping  (10^*  cm“^) 

Loss  (cm 

Ino.2Gao.8As 

3.5691 

none 

N/A 

GaAs 

3.5256 

none 

N/A 

Alo.3Gao.7As 

3.3622 

none 

N/A 

Alo.65Gao.35As 

3.1637 

none 

N/A 

Alo.98Gao.O2As 

2.9713 

none 

N/A 

ALO 

1.55 

none 

N/A 

n-GaAs 

3.5256 

2 

10 

n-Alo.92Gao.08  As 

3.0067 

2 

10 

p-GaAs 

3.5256 

1 

11.5 

p-Alo.92Gao.08  As 

3.0067 

1 

11.5 

p^-GaAs 

3.5256 

5 

57.5 

Table  3.1  Material  Parameters  used  for  WIM  Simulations 


The  lowest  frequency,  or  fundamental,  VCSEL  lasing  mode  is  analogous  to  the  HEMn 
(or  HEii)  hybrid  waveguide  mode.  Here  the  first  and  second  subscripts  refer,  respectively, 
to  the  azimuthal  (m)  and  radial  mode  numbers.  The  HEMn  mode  is  the  most  plane-wave 
like  of  all  the  propagating  bound  modes,  despite  containing  both  and  [58,69].  It  is 
also  the  only  waveguide  mode  having  a  radial  intensity  distribution  with  a  maximum  at  the 
center.  Following  this  terminology,  I  refer  to  the  next  higher-order  VCSEL  modes  as  HEM21, 
TEoi,  and  TMqi  modes;  these  modes  make  up  the  degenerate  LPn  mode  under  the  “linear 
polarization”  approximation  [70].  All  three  of  these  modes  feature  a  radial  intensity  profile 
with  a  null  at  the  center. 

3.1.1  Field  Profile.  In  Figures  3.3  and  3.4,  I  plot  the  longitudinal  standing  wave 
intensity  profile  for  the  etched-post  and  oxide-apertured  VCSELs.  For  the  oxide-apertured 
structure,  the  contrast  between  DBR  layers  is  small  and  the  fields  penetrate  deeply  into 
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Figure  3.1  Quasi-3D  plot  of  the  etched-post  VCSEL  index  profile. 
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Figure  3.2  Quasi-3D  plot  of  the  oxide- apertured  VCSEL  index  profile. 
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the  top  DBR.  Therefore  the  large  index  difference  between  the  semiconductor  cavity  and 
surrounding  oxide  is  heavily  weighted  in  (2.25),  giving  a  large  discontinuity  in  (e"^)  between 
the  inner  and  outer  radial  regions.  As  a  result,  the  fields  within  the  oxide-aperture  VCSEL 
are  tightly  confined  to  the  core.  In  contrast,  for  the  etched-post  structure  the  large  index 
contrast  between  top  DBR  layers  allows  very  little  field  penetration,  resulting  in  a  smaller 
effective  index  difference  and  a  correspondingly  less  confined  field.  In  both  structures,  the 
fundamental  mode  is  more  confined  than  higher-order  modes,  and  smaller  VCSELs  exhibit 
less  confinement  than  larger  ones.  These  effects  can  also  be  seen  in  the  three-dimensional 
field  energy  distribution,  given  by  [7] 

If  -  2  -.21 

w  =  We  +  Wm  =  -^  SoSr  E  +  Hq  H  .  (3.1) 

I  plot  w  for  the  HEMn  and  TEoi  modes,  for  both  structures,  in  Figures  3.5  -  3.8.  Note  that 
the  familiar  standing  wave  pattern  does  not  appear  in  Figures  3.5  -  3.8.  This  is  a  result  of 
the  fact  that  these  energy  plots  include  both  an  electric  and  a,  magnetic  field  contribution. 
Nodes  of  E  (z)  correspond  to  antinodes  of  H  (z)  and  vice-versa,  effectively  “washing  out” 
the  standing  wave  and  leaving  only  the  envelope  function. 


S.1.2  Transverse  Confinement  Factor.  The  transverse  confinement  factor  F*"^  is 


usually  defined  as 


f\Efds  ’ 


where  the  integral  in  the  numerator  is  over  the  transverse  extent  of  the  active  region  and  the 
integral  in  the  denominator  is  over  the  entire  transverse  extent  of  the  field.  Figure  3.9  shows 
estimates  of  the  transverse  confinement  factor  for  the  first  two  modes  of  the  etched-post  and 


oxide-apertured  VCSELs  as  a  function  of  cavity  radius.  The  estimates  are  generated  from 
(3.2),  using  E  =  E{p)  E<j,{p)  for  TE  modes  (2.63)  and  E  — »•  H^{p)  for  TM  modes  (2.64); 
for  hybrid  modes  I  use  both  E^{p)  and  H^{p)  and  average  the  results.  I  use  the  effective  (f> 
components  since  they  are  already  averaged  over  z  and  are  representative  of  the  total  field 
intensity  profile.  As  the  cavity  radius  decreases,  more  of  the  field  intensity  leaks  out  of  the 
active  region  and  the  confinement  factors  drop  monotonically.  This  behavior  becomes  more 
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Figure  3.3  Index  and  standing  intensity  profile  along  the  axial  direction  for  a  1.4  (im  radius 
etched-post  VCSEL. 
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Figure  3.4  Index  and  standing  intensity  profile  along  the  axial  direction  for  a  1.4  nm  radius 
oxide-apertured  VCSEL. 
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Figure  3.5  HEMn  mode  energy  profile  for  a  1.4  jim  radius  etched-post  VCSEL.  The  energy 
distribution  on  the  top  surface  is  amplified  in  order  to  illustrated  the  emitted 
mode. 
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Figure  3.6 


TEoi  mode  energy  profile  for  a  1.4  nm  radius  etched-post  VCSEL.  The  energy 
distribution  on  the  top  surface  is  amplified  in  order  to  illustrated  the  emitted 
mode. 
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Figure  3.7  HEMn  mode  energy  profile  for  a  1.4  iim  radius  oxide-apertured  VCSEL.  The 
energy  distribution  on  the  top  surface  is  amplified  in  order  to  illustrated  the 
emitted  mode. 


3-10 


Figure  3.8  TEoi  mode  energy  profile  for  a  1.4  fim  radius  oxide-apertured  VCSEL.  The 
energy  distribution  on  the  top  surface  is  amplified  in  order  to  illustrated  the 
emitted  mode. 
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pronounced  for  higher-order  modes,  as  illustrated  in  Figure  3.10.  My  oxide-apertured  de¬ 
vices  confine  the  optical  mode  more  strongly  to  the  active  region  and  have  higher  transverse 
confinement  factors  than  my  etched-post  structures.  But  the  larger  rate  of  change  in  confine¬ 
ment  factor  for  the  etched-post  VCSEL,  as  shown  in  3.9,  yields  better  modal  discrimination 
via  For  example,  a  1  /xm  radius  device  has  AF*"^  =  r*’^(HEMll)  —  r*'^(TE01)  =  0.0759 
for  the  etched-post  VCSEL,  compared  to  only  0.0379  for  the  oxide-apertured  VCSEL.  This 
illustrates  the  effectiveness  of  employing  a  small  radial  index  difference  to  introduce  mode 
selective  losses  and  enhance  single  mode  lasing  [27].  However,  single  mode  lasing  may  also 
be  enhanced  by  exploiting  the  spectral  separation  between  transverse  modes.  I  will  show 
this  in  the  following  section. 

S.l.S  Resonant  Wavelength.  In  Figures  3.11  and  3.12  I  plot  the  resonant  wave¬ 
lengths  as  a  function  of  cavity  radius  for  various  modes  in  etched-post  and  oxide-apertured 
structures.  The  resonant  wavelength  blueshifts  as  the  oxide  or  post  diameter  shrinks,  a 
dramatic  departure  from  plane-wave  results.  This  wavelength  shift  can  be  easily  explained 
by  examining  the  weighted  dispersion  relations  (2.21)  and  (2.24),  both  of  which  take  the 
functional  form 

Although  both  {kj)  and  (/?^)  change  as  the  radius  shrinks,  (jSl)  remains  very  close  to  its  plane- 
wave  value.  Therefore,  as  {kj)  increases  from  its  plane-wave  value  of  zero,  w  increases,  leading 
to  the  blueshift  illustrated  in  Figures  3.11  and  3.12.  This  effect  has  been  previously  estimated 
using  simpler  approximations  [21],  but  the  self-consistent  WIM  results  are  more  accurate, 
as  shown  in  Figures  3.21  and  3.22.  To  construct  low-threshold  microcavity  VCSELs,  the 
quantum  well  emission  peak  must  be  matched  to  the  blueshifted  cavity  resonance  of  the 
desired  lasing  mode. 

The  resonant  wavelength  changes  more  quickly  with  radius  in  my  oxide-apertured  VC¬ 
SELs  as  compared  with  my  etched-post  structures.  This  occurs  because  the  oxide-apertured 
device  exhibits  a  larger  difference  in  effective  index  between  the  inner  and  outer  radial  re¬ 
gions,  resulting  in  a  larger  field  confinement  and  a  correspondingly  larger  value  for  (kj) 
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Figure  3.9  Transverse  confinement  factor  for  the  first  two  modes  of  the  etched-post  and 
oxide-apertured  VCSEL. 
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Figure  3.10  Transverse  confinement  factor  for  the  fundamental,  and  a  few  sample  higher 
order  modes  for  the  oxide-apertured  VCSEL. 
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Figure  3.11  Resonant  wavelength  for  the  first  two  modes  of  the  etched-post  and  oxide- 
apertured  VCSEL. 
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Figure  3.12  Resonant  wavelength  for  the  fundamental,  and  a  few  sample  higher  order 
modes  for  the  oxide-apertured  VCSEL. 


in  the  waveguide  core.  As  a  result,  the  oxide-apertured  structure  provides  more  spectral 
discrimination  between  the  VCSEL  resonant  modes.  For  example,  in  a  1  /im  VCSEL 
AA  =  A(HEMll)  —  A(TEOl)  =  68.9  A  for  the  oxide-apertured  structure  as  compared  to 
59.8  A  for  the  etched-post.  This  effect  might  be  exploited  to  obtain  low  threshold,  single 
mode  lasing.  As  discussed  in  my  patent  (pending)  [22,55],  if  the  gain  spectrum  is  sufficiently 
narrow,  spectral  mode  discrimination  may  be  used  to  select  the  fundamental  lasing  mode. 
It  is  possible  that  a  quantum  dot  gain  region  will  provide  the  necessary  narrow-band  gain. 
In  contrast,  as  discussed  in  the  last  section,  traditional  single  mode  designs  rely  on  thresh¬ 
old  gain  discrimination.  Here  the  difference  in  threshold  gain  between  the  fundamental  and 
higher  order  modes  is  large  and  for  low  powers,  the  VCSEL  will  lase  in  a  single  mode.  This 
type  of  design  is  often  referred  to  as  anti-guiding.  The  principal  drawback  of  this  method 
is  an  increased  threshold  for  all  the  modes.  By  using  the  spectral  separation,  this  problem 
may  be  avoided. 

Finally,  I  plot  the  components  of  the  quasi-degenerate  LPn  mode  in  Figure  3.13.  The 
results  show  a  small,  but  non-zero,  splitting  of  the  mode  for  aperture  radii  <«  0.85  /xm, 
indicating  where  the  LP  mode  approximation  begins  to  break  down.  Overall,  however,  the 
LP  mode  approximation  is  very  good  throughout  the  convergent  range  of  the  calculation.  I 
exploit  this  fact  to  simplify  the  weighted  index  method  to  a  scalar  technique  in  Chapter  V . 

S.1.4  Mirror  and  Absorption  Loss.  Within  the  WIM  framework,  I  approximate 
the  VCSEL  cavity  modes  as  superpositions  of  cylindrical  waves.  Each  of  these,  in  turn,  can 
be  viewed  as  a  superposition  of  TE  and  TM  plane  waves  propagating  at  an  angle 

{0)  =  arctan  (3-4) 

to  the  z  axis.  As  the  cavity  radius  decreases,  the  effective  transverse  propagation  constant 
{kp)  increases  for  all  modes  and,  consequently,  the  average  angle  of  incidence  for  the  compo¬ 
nent  plane  waves  impinging  on  each  DBR  interface  increases.  The  power  reflectivity  for  TE 
waves  increases  monotonically  with  angle  until  the  total  internal  reflection  angle  is  reached, 
while  the  power  reflectivity  for  TM  waves  decreases  monotonically  until  the  Brewster  angle 
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Figure  3.13  Resonant  wavelength  for  the  components  of  the  quasi-degenerate  LPn  mode, 
illustrating  the  point  at  which  the  degeneracy  is  broken. 
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is  reached  [65].  As  a  result,  the  TE  wave  components  encounter  more  reflective  DBR  mirrors 
as  the  cavity  radius  shrinks,  while  the  TM  wave  components  encounter  less  reflective  DBRs; 
Figure  3.14  illustrates  this  behavior  for  the  oxide-apertured  structure.  Therefore  the  TE 
wave  components  suffer  less  mirror  loss,  the  TM  wave  components  more,  and  the  mirror 
losses  for  pure  TE  VCSEL  modes  decreases  with  cavity  radius,  while  the  mirror  losses  for 
TM  modes  increases.  For  HEM  modes — comprised  of  both  TE  and  TM  components — the 
change  in  reflectivity  approximately  cancels  out.  All  of  these  trends  are  evident  in  Figures 
3.15  -  3.18,  which  show  the  mirror  and  absorption  losses  divided  by  the  total  confinement 
factor  r*®*, 

^  l^(mirror)  ^  ^(absorption) J  =  47r/Cactive/A,  (3.5) 

for  the  fundamental  and  several  higher  order  VCSEL  transverse  modes.  The  quantity  I  plot 
is  precisely  the  threshold  gain  in  the  absence  of  diffraction  loss. 

The  results  shown  in  Figures  3.14  -  3.16  agree  with  prior  in-plane  laser  studies  that 
show  a  higher  facet  reflectance  for  TE  than  TM  modes  and  a  similar  propagation  constant  {^) 
for  both  TE  and  TM  modes  [2,34,38,41].  For  both  VCSELs  and  in-plane  lasers,  these  trends 
in  mirror  reflectivity  become  more  pronounced  for  higher-order  modes.  The  minima  in  the 
TE  curves  in  Figures  3.15  and  3.16  occurs  at  the  point  where  the  change  in  confinement  factor 
balances  the  change  in  mirror  reflectance.  For  TM  and  HEM  modes,  no  minima  exists  since 
the  mirror  reflectance  is  approximately  constant  (HEM)  or  decreases  (TM)  with  decreasing 
radius.  The  radius  at  which  each  curve  terminates  indicates  the  “minimum”  aperture  size  for 
each  optical  mode.  This  size  effect  has  been  previously  estimated  by  simpler  calculations  [21], 
but  the  WIM  values  should  be  more  accurate.  Since  my  oxide-apertured  structure  confines 
the  fields  better  than  the  etched-post  structure,  it  supports  bound  modes  at  smaller  cavity 
radii.  This  completes  the  overview  of  the  WIM  features.  I  now  turn  to  a  more  challenging 
task  for  the  method,  prediction  of  experimental  VCSEL  data. 

3.8  Experimental  Verification  of  Blueshift  Predictions 

I  compare  measured  blueshifts  for  two  different  device  designs  against  predictions  from 
my  weighted  index  model  and  a  simpler  effective  index  model.  Both  models  require  only  one 
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Figure  3.15 


Mirror  and  absorption  loss  for  the  fundamental,  and  transverse  electric  and 
magnetic  modes  for  the  etched-post  VCSEL. 
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Figure  3.16 


Mirror  and  absorption  loss  for  the  fundamental,  and  transverse  electric  and 
magnetic  modes  for  the  oxide-apertured  VCSEL. 
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adjustable  parameter — the  cavity  thickness  is  adjusted  to  match  the  Fabry-Perot  resonance 
of  broad-area  devices.  I  find  that  the  weighted  index  calculation  reproduces  the  dependence 
of  lasing  wavelength  on  aperture  radius  remarkably  well,  while  the  effective  index  model 
underestimates  the  blueshift. 


3.2.1  Effective  Index  Method.  A  simple  and  popular  approximation  to  calculate 
the  blueshift  is  the  effective  index  method.  In  this  approach,  an  axial  wavevector  is  calcu¬ 
lated  for  each  layer  by  propagating  plane-waves  through  the  dielectric  layer  stack,  using  the 
refractive  indices  at  p  =  0  for  each  layer.  The  resulting  (scalar)  electric  field  profile  ^(2)  is 
then  used  to  compute  effective  indices 


^  I  dz  |^(z)|  n^(p  =  0,^:) 

Idz\az)f 


^  J  dz  |^(2r)|  n^jp  =  00,  z) 


(3.6) 


corresponding  to  the  core  and  cladding  indices  for  an  “effective”  cylindrical  dielectric  waveg¬ 
uide.  An  effective  axial  wavevector  {/3)  is  computed  similarly.  The  radial  wave  equation 
for  this  waveguide  is  then  solved,  taking  regular  Bessel  functions  inside  the  aperture  and 
evanescent  solutions  outside.  This  procedure  fixes  the  radial  wavevector  fe,  and  the  resonant 
frequency  may  then  be  calculated  from  a  weighted  dispersion  relation.  Augmentations  of 
this  model  and  applications  to  oxide-apertured  VCSELs  are  given  in  References  [26]  and  [21]. 


3.2.2  Experiment.  I  compare  both  effective-  and  weighted-index  calculations  of  the 
blueshift  against  measurements  performed  on  two  different  VCSEL  structures.  To  conform 
with  the  most  often  measured  modes^,  I  label  the  VCSEL  modes  in  accordance  with  the  linear 
polarization  (LP)  convention  [58].  Therefore,  the  fundamental  lasing  mode — characterized 
by  a  peak  in  intensity  at  the  mode  center — is  labeled  as  LPoi.  Similarly,  the  first  higher-order 
VCSEL  mode — characterized  by  a  null  at  the  mode  center — is  denoted  by  LPn.  These  LP 
modes  approach  the  exact  fiber  optic  modes  as  the  index  difference  between  the  core  and 


^The  fact  that  the  LP  modes,  not  the  vector  modes,  are  most  often  observed  is  a  result  of  two  factors:  (1) 
the  relatively  small  weighted  index  difference  between  the  inner  and  outer  radial  regions,  which  supports  the 
“weakly  guiding”  assumption  upon  which  the  LP  modes  are  based;  and  (2)  the  fact  that,  for  small  aperture 
radii,  diffraction  dominates  any  difference  in  threshold  between  the  quasi-degenerate  constituents  of  an  LP 
mode.  As  a  result  of  these  arguments,  I  change  the  WIM  to  solve  for  LP  modes  in  Chapter  V. 
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cladding  becomes  small:  the  LPoi  mode  includes  the  doubly  degenerate  HEMn  modes,  and 
the  LPii  mode  includes  the  TEoi,  TMoi,  and  doubly  degenerate  HEM21  modes. 

The  first  VCSEL  design  is  a  1-A  cavity,  single  oxide  aperture  device,  designed  for  780  nm 
emission  [66].  The  active  region  has  four  80  A  Alo.iiGao.89 As  quantum  wells  separated  by  80 
A  Alo.3Gao.7As  barriers.  On  either  side  of  the  active  region  lies  an  Al0.sGao.5As  spacer  layer, 
completing  the  1-A  cavity.  The  upper  and  lower  distributed  Bragg  reflectors  (DBRs)  consist 
of  26  and  40.5  periods  of  Alo.3Gao.7As/Alo.9Gao.1As,  respectively.  All  DBR  interfaces,  as  well 
as  the  spacer-to-DBR  interfaces,  are  linearly  graded.  To  facilitate  oxidation,  the  upper-DBR 
low-index  layer  adjacent  to  the  cavity  is  AlAs.  When  oxidized,  this  layer  contains  a  square 
aperture.  The  entire  structure  is  grown  on  a  GaAs  substrate.  The  structure  is  illustrated  in 
Figure  3.19. 

To  investigate  the  blueshift,  five  devices  were  fabricated  from  this  layer  structure:  four 
with  small  apertures  (<  5  fim  side  length)  and  one  broad  area  device  with  a  13  jwm  x  13  fxm 
aperture.  All  five  devices  were  located  within  a  0.5  mm  x  0.5  mm  area  of  the  wafer,  ensuring 
excellent  device  uniformity.  While  electrically  pumping  the  devices  at  half  threshold,  lumi¬ 
nescence  data  was  collected  with  an  optical  multichannel  analyzer.  The  electroluminescence 
data  from  the  broad-area  device  showed  a  fundamental  mode  at  782  nm,  which  I  interpret 
as  the  Fabry-Perot  resonance  of  the  VCSEL  layer  structure.  The  remaining  measurement 
results  for  the  four  small-aperture  devices  are  plotted  in  Figure  3.20;  a  semilog  scale  is  em¬ 
ployed  to  bring  out  the  details  in  the  electroluminescence  response.  For  each  aperture  size  at 
least  two  distinct  resonances  are  visible.  As  the  aperture  size  increases,  more  resonant  modes 
appear,  and  the  spacing  between  them  decreases.  To  precisely  locate  the  mode  resonances, 
the  peaks  in  Figure  3.20  were  fitted  with  Lorentzian  lineshape  functions.  Moving  from  right 
to  left,  the  first  two  resonances  represent  the  LPqi  and  LPn  modes. 

The  second  VCSEL  design  is  a  1-A  cavity,  double-oxide-aperture  device,  designed  for 
850  nm  emission  [14].  The  active  region  has  five  80  A  GaAs  quantum  wells  separated 
by  Alo.3Gao.7As  barriers.  The  upper  and  lower  DBRs  consist  of  26  and  36.5  periods  of 
Alo.i6Gao.84As/Alo.92Gao.08As,  respectively,  with  graded  interfaces  between  the  layers  and 
between  the  DBRs  and  the  cavity.  Square  apertures  are  formed  by  oxidizing  182  A  thick 
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Figure  3.19  Illustration  of  the  780  nm  test  VCSEL. 
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Electroluminescense  (a.u.) 


Wavelength  (nm) 

Figure  3.20  Electroluminescence  data  taken  at  0.5  It^  for  four  different  VCSELs,  each 
labeled  by  square-aperture  side  length. 
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Alo.98Gao.o2 As  layers,  inserted  in  the  third  low-index  layer  away  from  the  cavity  in  both  the 
upper  and  lower  DBRs.  Due  to  oxidation  into  the  graded  regions,  the  final  oxide  layer  is 
approximately  300  A  thick.  The  entire  structure  is  grown  on  a  GaAs  substrate.  Oxide- 
apertured  devices  of  various  sizes  were  fabricated  as  described  in  Reference  [14],  and  the 
resonant  frequency  of  the  fundamental  cavity  mode  was  measured  as  a  function  of  aperture 
size.  The  fundamental  resonance  of  broad-area  devices  occurred  at  851  nm. 


S.2.S  Blueshift  Comparison.  Figures  3.21  and  3.22  compare  both  effective-index 
and  weighted- index  blueshift  calculations  against  the  measurements  for  the  LPqi  and  LPn 
modes  [57].  Since  the  measurements  are  for  a  square-aperture  device,  while  the  calculation 
is  based  on  a  circular  aperture,  I  compare  the  results  by  matching  cross-sectional  area. 
Therefore,  the  abscissa  in  Figures  3.21  and  3.22  is  effective  aperture  radius: 


reff  = 


< 


a 

r 


measurement 

calculation 


(3.7) 


where  a  is  the  side  length  of  the  square  aperture  and  r  is  the  device  radius  in  the  calculation. 
Although  this  substitution  cannot  be  rigorously  justified,  it  does  provide  excellent  agreement 
with  the  data.  The  theoretical  results  shown  in  Figures  3.21  and  3.22  were  generated  by 
analyzing  device  structures  almost  the  same  as  those  measured.  (I  approximated  the  graded 
interfaces  with  abrupt  heterojunctions.)  To  calibrate  the  models  to  the  actual  structures 
grown,  I  adjusted  the  cavity  thicknesses  in  both  VCSEL  designs  so  that  the  calculated 
Fabry- Perot  resonance  matched  measurements  performed  on  broad  axea  devices  782  nm 
for  the  first  VCSEL  design,  «  851  nm  for  the  second).  Apart  from  these  layer  thicknesses, 
my  calculations  included  no  adjustable  parameters.  Note  that  adjusting  these  layers  serves 
primarily  to  shift  the  calculated  curves  up  or  down  in  wavelength:  it  does  not  affect  the 
shape  as  a  function  of  radius.  All  calculated  curves  terminate  when  a  radially-bound  mode 
can  no  longer  be  found.  This  provides  an  estimate  of  the  minimum  aperture  radius  required 
to  laterally  confine  the  mode.  The  effective  index  step  ncore  —  «clad?  as  calculated  from  (3.6), 
was  equal  to  0.033  for  the  first  VCSEL  design  and  0.061  for  the  second. 
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Figure  3.21  Calculated  and  measured  resonant  wavelengths  of  the  fundamental  (LPoi, 
upper  curves)  and  first  higher-order  (LPn,  lower  curves)  modes  as  a  function 
of  aperture  size  for  the  first  VCSEL  design;  discrete  points  represent  measured 
values.  The  inset  shows  the  longitudinal  refractive  index  and  field  intensity 
profiles  near  the  cavity. 
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Figure  3.22  Calculated  and  measured  resonant  wavelengths  of  the  fundamental  (LPoi) 
mode  as  a  function  of  aperture  size  for  the  second  VCSEL  design.  The  inset 
shows  the  longitudinal  refractive  index  and  field  intensity  profiles  near  the 
cavity. 
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From  Figures  3.21  and  3.22,  I  see  that  the  WIM  calculation  accurately  predicts  the 
dependence  of  wavelength  on  aperture  size,  matching  the  measurements  remarkably  well. 
The  simple  effective  index  calculation,  on  the  other  hand,  significantly  underestimates  the 
blueshift  for  both  the  fundamental  (LPqi)  and  first  higher-order  (LPn)  mode.  The  small  dif¬ 
ferences  between  the  WIM  calculations  and  the  measurements  in  Figures  3.21  and  3.22  could 
result  from  several  factors.  First,  the  model  is  restricted  to  piecewise-constant  refractive  in¬ 
dices,  while  the  actual  structures  have  graded  interfaces.  Second,  the  modeling  of  a  square 
aperture  with  a  circular  one  could  introduce  error.  Third,  the  weighted  index  method  is  itself 
not  exact,  since  it  assumes  separable  solutions.  Non-separable  effects — such  as  diffractive 
loss — are  not  included,  and  these  effects  could  influence  the  resonant  wavelengths.  Last,  the 
measured  aperture  sizes  are  only  accurate  to  within  ±  0.2  ^m. 

I  have  now  completed  an  “initial  validation”  of  the  weighted  index  method,  verifying 
the  resonance  predictions  against  experiment.  The  ability  to  predict  mode  blueshifts  is 
certainly  valuable  in  device  design,  allowing  for  accurate  alignment  of  the  cavity  resonance 
to  the  gain  peak.  However,  the  principal  optical  parameter  of  interest  is  the  threshold  gain. 
Since  the  WIM  is  inherently  limited  by  the  separability  assumption,  it  may  not  be  used 
directly  to  compute  diffraction  effects.  Therefore,  I  forfeit  the  speed  and  efficiency  of  the 
method  and  turn  to  a  computational  technique.  I  present  my  ntw  finite  element  model — 
capable  of  computing  all  optical  VCSEL  parameters — in  the  next  chapter. 
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IV.  Vector  Finite  Element  Method 

In  this  chapter,  I  present  a  quasi-exact  finite  element  method  (FEM)  model  for  analyzing 
the  optical  modes  of  microcavity  VCSELs  [53,54,56].  Although  computationally  intensive, 
this  model  has  the  distinct  advantage  that  it  calculates  the  total  optical  loss,  including 
diffraction.  This  model  is  based  on  a  variational  solution  of  the  vector  Helmholtz  equation 
in  microcavity  geometries.  The  results  of  this  model  not  only  allow  for  direct  calculation  of 
lasing  mode  parameters,  but  also  a  better  understanding  of  the  underlying  physics  associated 
with  VCSEL  oxides.  To  derive  basic  oxide  design  rules,  I  apply  the  model  to  several  versions 
of  a  novel  oxide-aperture,  oxide  mirror  VCSEL  [67].  I  find  that  the  diffraction  loss  may 
be  qualitatively  understood  a.s  a  coupling  between  the  lasing  eigenmode  and  the  parasitic 
mode  continuum  [20].  Within  this  interpretation,  the  diffraction  is  a  result  of  three  basic 
physical  effects:  (1)  the  amount  of  transverse  confinement  afforded  to  the  optical  mode,  (2) 
the  relative  alignment  between  the  propagation  vectors  of  the  eigenmode  and  the  radially 
radiating  parasitic  modes,  and  (3)  the  density  of  parasitic  modes  in  the  spectral  vicinity  of 
the  eigenmode. 

4-1  Vector  Field  Equations  and  the  Variational  Form 

I  want  to  find  the  electric  and  magnetic  field  profiles,  the  resonant  wave¬ 

length  ( Ao),  and  the  threshold  gain  (gth)  for  each  cavity  mode  in  azimuthally-symmetric 
VCSEL  structures.  For  this  I  must  solve  Maxwell’s  equations  subject  to  appropriate  bound¬ 
ary  conditions  at  each  material  interface.  The  steady-state,  time-harmonic  electric  field 
satisfies  the  vector  Helmholtz  equation  (in  MRS  units) 

V  X  —V  X  E  —  klsrE  =  —iuju^j,  (4.1) 

f^r 

where  E,  H,  and  the  electric  current  J  depend  on  time  as  e®"*  (u  =  27rc/A);  eo)  Mo 5 

are  the  free  space  permittivity,  permeability,  and  propagation  constant,  respectively.  It  can 

be  shown  [37]  that  a  weak  solution  to  (4.1)  may  be  obtained  by  extremizing  the  functional 


S{E)  =  (4.2) 

+ioJfio  jj{n.E).H,.^jjl  E  •  J  dv. 

(By  extremize,  I  mean  forcing  =  0,  where  6  is  the  variational  operator.)  Here  U  and  F 
are  the  problem  domain  and  boundary,  respectively.  For  azimuthally-symmetric  structures, 
the  material  parameters  ( and  Sr )  are  functions  of  p  and  z  only,  and  I  may  separate  out 
the  ^  dependence  in  (4.2)  by  assuming 

E^  cos  (m^) ,  (4.3) 

and 

Ep,  Ez  ^  sin  {m4>) ,  (4.4) 

where  m  is  the  azimuthal  mode  number.  For  integer  m,  the  integrals  over  =  [0, 27r]  yield 
a  constant  factor  which  may  be  ignored,  eifectively  reducing  the  dimension  of  the  problem 
from  three  to  two.  In  addition,  since  different  azimuthal  modes  are  orthogonal,  I  may  deal 
with  each  value  of  m  independently. 

For  lasing  mode  analysis,  I  set  the  source  current  J  to  zero,  making  (4.1)  a  source-free 
eigenmode  problem.  By  assuming  perfect  conducting  boundary  conditions  •  E  =  0^  on 
r,  which  I  justify  in  Section  4.3,  the  surface  integral  in  (4.2)  drops  out  and  I  am  left  with 

S{E)  =  JJ  xSy  (y  X  E^dv-kl  JJ  SrE-  Edv.  (4.5) 

Here  ft  represents  the  two-dimensional  domain  (of  the  VCSEL)  over  the  p  —  z  plane.  My  task 
is  then  to  find  the  E  field  which  extremizes  (4.5).  That  is,  I  must  find  the  E  that  satisfies 

11  —  (vxE\-(vx6E^dv-kllferE-  8Edv  =  0,  (4.6) 

./  V  n  Mt*  j 
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for  all  field  variations  6E.  Equation  (4.6)  is  a  generalized  eigenvalue  problem, 

S(E)-mE)  =  0.  (4.7) 


The  S  and  T  operators  are  defined  by 


5(^)  =  //  —  (v  X  •  (v  X  51®)  dv, 

(4.8) 

and 

T(^)  =  jj  6^  dv. 

(4.9) 

The  eigenvalue  ^  is  defined  as 

o 

(4.10) 

All  the  desired  mode  information  may  be  found  by  solving  (4.6).  For  each  mode,  the  eigen¬ 
value  is  the  square  of  the  (generally  complex)  free-space  propagation  constant,  which  is 
related  to  the  modal  wavelength  (Aq)  (4.22)  and  the  total  optical  loss  or  threshold  gain 
{9th)  (4.26).  The  eigenvectors  are  simply  the  time-harmonic  (vector)  electric  fields  ^E^  .  I 
approximate  the  solution  to  (4.6)  using  the  finite  element  method. 

Finite  Elements  and  Matrices 

The  formalism  developed  in  the  last  section  is  general  and  could  apply  to  any  of  several 
variational  approaches  to  solving  Maxwell’s  equations.  In  this  section,  I  narrow  my  attention 
to  the  finite  element  method  (FEM).  In  the  finite  element  method,  the  solution  to  (4.6)  is 
approximated  by  limiting  the  space  of  admissible  functions  E  to  the  linear  superposition  of 
a  finite  set  of  basis  functions  (generally  characterized  by  the  fact  that  they  are  non-zero  only 
over  a  subdomain  fie,  the  domain  of  mesh  element  e).  In  my  vector  FEM,  I  expand  the  fields 
over  a  basis  of  vector  functions  [49], 

N 

£=5^1*.  (4.11) 

«=1 
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Here  N  is  the  total  number  of  basis  functions  in  the  expansion,  x,-  are  unknown  coefficients, 
and  are  the  vector  basis  functions.  These  functions  are  second  order  node,  edge,  and 
face  element  functions,  given  in  Appendix  C.  They  are  specifically  designed  to  model  m  =  1 
modes.  Substituting  (4.11)  into  (4.6)  and  exchanging  the  order  of  summation  and  integration 
yields 

^  ^*  {//  ^  ^  ^  ^  I Iq  =  0-  (4.12) 

Ensuring  (4.12)  holds  for  all  (same  set  of  functions  as  tj;-),  ensures  it  will  hold  for  any 
linear  superposition  of  the  finite  basis  analogy  of  8E.  The  N  equations  represented  by 

(4.12)  are  exactly  (4.7)-(4.9)  taken  over  a  finite  basis,  written  conveniently  in  linear  algebra 
notation  as 

SX-^TX  =  0,  (4.13) 

where, 

Sij  =  jj  $,)  ■  (v  X  dv,  (4.14) 

kj  =  //  dv,  (4.15) 

J  Jn 

and  the  eigenvalue  (^  definition  (4.10)  remains  unchanged. 

I  define  my  basis  set  over  a  triangular  mesh  in  the  p  —  z  plane,  as  illustrated  in  Figure 
4.1.  By  using  (randomly  shaped)  triangular  elements,  I  can  accommodate  general  VCSEL 
designs — including  for  example  tapered  oxides — and  avoid  creating  an  artificial,  mesh-driven 
predisposition  to  any  given  vector  field  component.  Due  to  the  form  of  my  basis  expansion 
and  the  use  of  absorbing  regions  (discussed  below),  S  and  T  will  be  very  large  {N  ~  50, 000), 
sparse,  non-Hermitian  matrices.  As  a  result,  special  matrix  techniques  are  required  to  solve 

(4.13) .  In  practice  this  is  the  most  challenging  part  of  the  finite  element  solution,  and 
certainly  the  most  time  consuming.  I  chose  to  solve  the  eigenvalue  problem  using  an  iterative 
Arnoldi  algorithm  [46, 62]  with  spectral  transformation;  this  algorithm  allows  me  to  search 
for  “mildly  complex”  (e.g.,  Re(^)  »  Ini(^))  eigenvalues  over  a  given  range  of  the  real  axis. 
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Figure  4.1  Sample  finite  element  mesh  for  an  oxide-apertured,  oxide  DBR  VCSEL. 
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4-3  Mesh  Termination 

By  far  the  greatest  potential  error  source  in  the  finite  element  VCSEL  analysis  is 
mesh  termination.  Due  to  the  incomplete  optical  confinement  of  both  the  oxide  apertures 
and  the  VCSEL  mirrors,  the  true  domain  of  a  VCSEL  field  solution  extends  to  infinity  in 
all  directions.  Hence,  an  artificial  mesh  termination  is  required  for  FEM  application.  A 
properly  designed  mesh  termination  for  this  problem  must  mimic  the  unbounded  nature  of 
the  domain,  eliminating  nonphysical  reflections  at  the  mesh  edge.  In  addition,  since  the 
FEM  computational  demand — ^both  CPU  time  and  memory — scales  super-linearly  with  N,  I 
prefer  to  place  my  termination  as  close  as  possible  to  the  primary  domain  of  interest,  thereby 
minimizing  the  amount  of  “wasted”  mesh  space. 

To  terminate  my  mesh,  I  insert  an  artificial  absorbing  layer  (AL)  between  the  principle 
problem  (VCSEL)  domain  (Dy)  and  the  problem  boundary  (F)  [37].  This  layer  allows  me  to 
use  perfect  conductor  boundary  conditions  on  F,  cis  I  assumed  in  Section  4.1,  eliminating  the 
surface  integral  term  in  the  variational  form  and  dramatically  simplifying  the  FEM  analysis. 
I  define  fly  by  a  rectangular  region  in  the  p  —  z  plane,  bounded  by  the  bottom  mirror- 
to-substrate  plane  (Fb)  and  the  top  mirror- to- air  plane  (Fj)  in  the  z  direction,  and  the 
transverse  lasing  mode  size  (F5)  in  the  p  direction.  I  determine  Fs  a-posteriori,  as  discussed 
in  Appendix  D.3.  I  surround  fly  by  the  AL  as  illustrated  in  Figure  4.2,  where  the  AL  domain 
is  the  union  of  the  top,  side,  and  bottom  AL  regions  {CIal  =  flrFlflsHflB).  Due  to  the  high 
reflectance  of  the  distributed  Bragg  reflectors  (DBRs),  my  main  concern  is  absorbing  any 
radiation  incident  on  the  radial  boundary  F5.  Therefore,  I  focus  my  analysis  on  the  radial 
AL  (D5);  the  optimal  axial  AL  design  (Dr  and  Q,b)  falls  out  of  the  radial  analysis  as  I  show 
in  Appendix  D. 

A  basic  requirement  for  minimizing  reflection  is  that  the  impedance  of  the 

absorbing  layer  in  each  region  must  match  the  radially  adjacent  VCSEL  region.  I  enforce 
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Figure  4.2  Illustration  of  the  FEM  VCSEL  problem  domain:  is  the  VCSEL  domain, 

and  Cls,  fir?  and  Hb  are  the  side,  top,  and  bottom  absorbing  layer  (AL)  do¬ 
mains,  respectively.  The  VCSEL  and  AL  domains  are  separated  by  the  bound¬ 
aries  Fs,  Fr,  and  Fb,  and  the  entire  problem  domain  is  bounded  by  the  closed 
cylinder  F. 
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(4.16) 


this  condition  by  defining  the  absorbing  layer  material  parameters  as^ 

4!j{p)  =  and  =  a{p)pr,r 

Here  and  are  the  VCSEL  material  parameters  in  axial  region  j  (defined  by  Zj-i  < 
z  <  Zj)  radially  adjacent  to  the  AL  (e.g.,  just  to  the  left  of  F^),  and 

aip)  =  1  -  ib{p).  (4.17) 

The  absorbing  layer  performance  depends  entirely  on  the  function  b{p).  I  arrive  at  a 
suitable  function  b{p)  through  a  detailed  design  optimization.  First,  I  obtain  a  rough  estimate 
using  an  asymptotic  description  of  the  fields  and  reflections,  as  described  in  Appendix  D.l. 
Then  I  fine-tune  the  layer  by  minimizing  the  exact  reflection  values  as  obtained  by  a  rigorous 
transfer  matrix  calculation,  described  in  Appendix  D.2. 

4-4  Results  and  Analysis 

To  illustrate  the  application  of  the  FEM  for  VCSELs,  I  analyze  several  versions  of 
a  basic  870  nm  oxide-apertured,  oxide  DBR  VCSEL  [67].  The  VCSEL  has  five  and  a  half 
periods  of  Alo.3Gao.7As/Ala;03,  in  the  bottom  DBR  and  four  in  the  top  DBR.  The  cavity  is  lA 
thick  and  contains  a  single  600  A  GaAs  bulk  gain  region  centered  between  two  Alo.3Gao.7As 
barrier  layers.  A  300  A  AlAs  layer  is  included  in  the  top  barrier  to  form  an  oxide  aperture. 
Although  an  actual  oxide  aperture  formed  in  AlAs  would  have  a  square  cross-section,  I  treat 
it  as  circular  to  maintain  the  azimuthal  symmetry.  The  entire  structure  including  the  GaAs 
substrate  is  illustrated  in  Figure  4.1.  The  versions  of  this  structure  I  examined  are: 

•  lA-lTHIN  the  structure  as  described  above, 

•  lA-lTHICK  same  structure  with  a  600  A  thick  oxide  aperture, 

^Before  choosing  (4.16),  I  considered  a  diagonal  anisotropic  a(p)  for  a  perfectly  matched  layer  (PML) 
design  [45,63].  However,  I  found  that  the  extra  degrees  of  freedom  provided  no  advantage  for  absorbing 

general  outward  propagating  cylindrical  waves  of  the  form  (h^  (kp)  . 
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•  1A-2THICK  same  structure  with  a  600  A  thick  oxide  aperture  in  both  the  top  and 
bottom  barrier  regions,  and 

•  A/2-1THICK  same  DBRs  but  with  a  A/2  cavity  and  a  600  A  thick  oxide  aperture  in 
the  top  barrier  region. 

All  materials  indices  are  assumed  to  be  real  so  that  the  only  source  of  field  loss  is 
through  absorption  in  the  ALs.  For  this  particular  structure,  the  absorbing  layer  function 
h{p)  is  derived  in  Appendix  D.3. 


44.1  Lasing  Mode  Fields.  In  Figures  4.3  and  4.4  I  plot  |  E4, \  and  the  time  averaged 
electromagnetic  energy  density, 


(4.18) 


for  the  sample  lA-lTHIN,  =  0.4  pm  VCSEL;  these  plots  are  representative  of  the  general 
lasing  mode  profiles  for  all  the  VCSELs  I  tested.  The  fields  are  found  by  substituting  the 
lasing  mode  eigenvector  (X)  into  the  field  expansion  given  in  (4.11).  The  profile  (Figure 
4.3)  is  similar  to  the  familiar  standing  wave  profile  obtained  via  simple  scalar  field  techniques, 
however,  I  anticipate  my  result  is  more  accurate  due  to  the  full  vector  solution.  The  energy 
density  (4.4)  is  found  by  estimating  H  from  the  E  field  expansion  (4.11)  and  Faraday’s  law. 
Although  this  is  less  familiar  than  the  E^,  profile,  it  is  a  more  accurate  representation  of  the 
spatial  mode  energy  distribution  throughout  the  VCSEL. 


4.4-^  Confinement  Factor.  Using  the  spatial  mode  profile,  I  estimate  the  total  and 
transverse  confinement  factors  as 


and 


(4.19) 


(4.20) 
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Figure  4.3  Example  plot  of  \E^\  for  the  lA-lTHIN  structure  with  oxide  aperture  radius 
Poj.  =  0.4  jum.  To  increase  clarity,  the  figure  domain  is  smaller  than  the  calcu¬ 
lation  domain,  and  the  background  intensity  is  set  to  white. 
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Figure  4.4 
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Example  plot  of  the  stored  energy  density  w  for  the  lA-lTHIN  structure  with 
oxide  aperture  radius  =  0.4  pm.  To  increase  clarity,  the  figure  domain  is 
smaller  than  the  calculation  domain,  and  the  background  intensity  is  set  to 
white. 


respectively;  fly  is  the  VCSEL  volume,  and  flpump  and  (Igain  are  the  pumped  and  total  volume 
of  the  gain  region,  respectively.  For  purpose  of  calculation,  I  estimate  flpump  as  the  volume 
of  the  active  layer  inside  the  oxide  aperture  (e.g.,  with  p  <  pox)-  The  total  confinement 
factor  represents  the  percentage  of  the  mode  energy  overlapping  the  active  gain  region^.  I 
use  this  later  in  the  Fabry-Perot  laser  equation  to  estimate  material  threshold  gain  from  the 
total  modal  loss  (4.26).  The  transverse  confinement  factor  has  less  quantitative  application 
but  is  a  nice  indicator  of  how  well  the  lasing  mode  is  confined  in  the  transverse  dimension. 

I  plot  r*''  for  each  of  the  VCSEL  structures  in  Figure  4.5.  The  results  for  the  lA  cavity 
are  somewhat  intuitive:  for  any  given  radius  more  and  thicker  oxides  yield  greater  confine¬ 
ment.  For  the  A/2  case,  the  oxide  aperture  does  not  overlap  well  (in  the  z  direction)  with 
the  standing  wave.  Therefore,  the  transverse  confinement  is  relatively  small,  however,  the 
total  confinement  factor — due  to  the  short  cavity — is  very  large,  as  shown  in  Figure  4.6.  The 
high  contrast  DBRs  allow  less  field  penetration  and  therefore  a  very  high  total  confinement 
factor  for  all  four  test  cases,  compared  to  analogous  semiconductor  DBR  VCSELs. 

4.4-3  Resonant  Wavelength,  Total  Optical  Mode  Loss,  and  Threshold  Gain.  Due 
to  the  absorbing  layers,  the  eigenvalues  of  (4.10)  are  complex,  and  take  the  general  form 


=  ki=^  =  {ri  +  iqi?- 


*  c2 


(4.21) 


Here  and  qi  axe  the  real  and  complex  parts  of  the  (total)  propagation  constant  ki  of  mode 
i.  From  (4.21)  I  immediately  recognize  the  mode  resonance  as 


By  definition,  the  fields  vary  harmonically  as 


(4.22) 


exp  (  iuit)  =  exp  ( ikict)  =  exp  ( inct)  exp  {—qict) 


(4.23) 


^This  is  a  somewhat  different  definition  than  the  more  standard  definition  using  the  standing  wave 
intensity  \e\  V  However,  the  two  definitions  should  yield  similar  results  in  most  cases  , 
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Transverse  Confinement  Factor 


1 


Figure  4.5 


Transverse  confinement  factor  verses  oxide  aperture  radius  for  the  fundamental 
lasing  mode.  The  lines  are  cubic  spline  fits  of  the  discrete  calculation  data. 
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Figure  4.6  Total  confinement  factor  verses  oxide  aperture  radius  for  the  fundamental  lasing 
mode.  The  lines  are  cubic  spline  fits  of  the  discrete  calculation  data. 
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where  the  mode  frequency  is  Re(c<;t)  =  r,c  and  the  total  cavity  loss  rate  (1/seconds)  is 
li  =  qiC.  These  results  describe  a  three-dimensional  leaky  cavity,  losing  energy  at  a  rate  of  k. 
I  would  prefer  to  express  my  results  in  terms  of  the  more  familiar  Fabry-P erot  (F-P )  laser 
model.  To  do  this,  I  convert  ix>i  to  a  complex  propagation  constant 


n 


gain? 


(4.24) 


where  *  is  the  complex  conjugate,  required  such  that  a  lossy  Wj  maps  to  a  lossy  13^.  The 
imaginary  part  of  gives  the  field  loss  rate  (for  Im(;di)  >  0)  in  cm“^  Since  energy  is 
proportional  to  the  square  of  the  field,  the  intensity  lost  per  unit  length  (a,)  is  given  by 


ai  = 


(4.25) 


In  deriving  Oj  all  I  have  done  is  convert  energy  lost  per  time  to  energy  lost  per  propagation 
length.  To  apply  this  loss  rate  to  a  “textbook”  Fabry-Perot  laser,  the  cavity  length  would 
have  to  be  adjusted  in  accordance  with  the  blueshift  [26].  Finally,  I  calculate  the  threshold 
gain  (gth)  from  the  F-P  lasing  condition. 


(4.26) 


In  Figures  4.7  -  4.9  I  plot  the  modal  resonance,  total  optical  loss,  and  threshold  gain 
as  a  function  of  oxide  aperture  radius  for  all  four  test  VCSELs.  The  resonance  results  show 
the  now  familiar  blueshift,  with  more  and  thicker  oxide  apertures  yielding  a  larger  shift. 
Interpretation  of  the  loss  and  threshold  gain  curves  is  more  complicated.  I  discuss  these 
results,  in  detail,  in  Section  4.5. 


444  Cavity  Q  and  Spectral  Mode  Width.  From  w,-  and  a,,  I  estimate  the  cavity  Q 
and  spectral  half  width  at  half  maximum  for  each  of  the  eigenmodes  as  [65] 


Re(a;,) 

ri 

VpCXi 

2qi 

(4.27) 
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Figure  4.7  Lasing  mode  resonance  verses  oxide  aperture  radius.  The  lines  are  cubic  spline 
fits  of  the  discrete  calculation  data. 
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Figure  4.8 


Lasing  mode  total  optical  loss  verses  oxide  aperture  radius.  The  lines  are  cubic 
spline  fits  of  the  discrete  calculation  data. 


4-17 


Threshold  Material  Gain  (cm 
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Figure  4.9  Threshold  (material)  gain  verses  oxide  aperture  radius.  The  lines  axe  cubic 
spline  fits  of  the  discrete  calculation  data. 


4-18 


and 


8(jOi  «  lupail  =  Im(a;i), 


(4.28) 


respectively.  This  estimate,  along  with  the  a;,-  and  a,  results,  may  be  used  to  estimate 
the  discrete  optical  mode  spectrum  within  a  spontaneous  emission  rate  calculation.  This 
should  yield  results  more  accurate  than  the  typical  plane  wave  mode  expansion,  particularly 
for  the  case  of  highly  confined  VCSELs  where  the  resonant  modes  deviate  substantially  from 
plane  waves.  To  illustrate  this  application,  I  plot  the  eigenmodes  for  the  lA-lTHIN  VCSEL 
in  Figure  4.10.  Each  of  the  modes  has  an  associated  field  distribution  which  may  be  used  in 
the  emission  rate  equations^. 

4.4.5  Radial  and  Axial  Mode  Loss.  To  better  understand  the  optical  loss  sources, 
I  would  like  to  divide  the  total  optical  mode  loss  (a,)  into  a  mirror  loss  ,  due  to 

emission  out  the  ends  of  the  VCSEL,  and  a  diffraction  loss  ^  due  to  loss  out  the 

VCSEL  side.  To  do  this,  I  use  conservation  of  energy  and  the  relationship  between  a  and 
the  radiated  power  and  stored  energy  [7], 


Pv 

VpWv' 


(4.29) 


Here  Vp  is  the  field  phase  velocity. 


ExR 


)■*■} 


(4.30) 


is  the  total  time  averaged  (real)  power  exiting  the  VCSEL  (flv)  through  Ty,  and 


Wv  =  W^  + 


(4.31) 


^Note  that  in  some  cases  (not  shown)  the  lasing  mode  (the  mode  with  the  smallest  Oj)  does  not  have  the 
longest  wavelength.  This  is  not  a  problem,  since  the  longer  wavelength,  higher  loss  eigenmodes  have  very 
little  spatial  overlap  with  the  gain  region,  and  therefore  will  not  contribute  significantly  to  the  emission  rate. 


4-19 


Mode  Q 


5000 


4000 


3000 


2000 


1000 


0 


Figure  4.10  Illustration  of  the  discrete  m  =  1  mode  spectrum  for  the  lA-lTHIN  structure. 

The  modes  are  depicted  as  Gaussian  with  height  equal  to  the  mode  Q  and 
full  width  half  maximum  equal  to  the  line  width  defined  in  Section  4.4.4.  The 
mode  search  range  was  limited  to  co  values  corresponding  to  Aq  ranging  from 
870  nm  to  820  nm. 
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is  the  total  time  averaged  stored  energy.  This  consists  of  the  electric  and  magnetic  field 
energies, 


=  ^o£r, ja  e\v, 

(4.32) 

and 

(4.33) 

respectively.  Equation  (4.29)  is  derived  from  the  relationship  between  cavity  Q  and  loss  rate, 
adapted  to  the  F-P  laser  model.  I  find  in  practice  that  numerical  application  of  (4.29)  gives 
results  exactly  matching  those  obtained  with  (4.25),  verifying  the  general  approach.  It  is 
then  a  simple  matter  to  break  a  into  radial  and  axial  components. 

(diffraction) 

~  VpWv' 

(4.34) 

and 

(mirror)  _  “I"  Pb 

VpWv  ■ 

(4.35) 

Here 

Pv  =  Ps  +  Pt  +  Pbi 

(4.36) 

breaks  the  total  power  exiting  the  VCSEL  into  that  leaving  through  the  side  {Ps),  and  that 
leaving  through  the  top  and  bottom  surfaces  (Pj)  and  (Pb),  respectively.  These  powers  are 
calculated  by  restricting  the  integral  in  (4.30)  to  the  corresponding  surfaces.  Note  that  for 
these  calculations,  I  assume  the  VCSEL  materials  are  lossless  and  calculate  “cold  cavity” 
radiative  loss  parameters.  Although  absorptive  loss  may  be  significant  in  operating  VCSELs, 
this  calculation  should  give  a  good  estimate  of  the  diffraction  loss.  The  calculation  can  be 
easily  modified  to  fold  in  absorptive  losses  if  required. 

In  Figures  4.11  -  4.13,  I  plot  the  radial  mode  loss,  the  axial  mode  loss,  and  the  radial 
percentage  of  the  total  mode  loss,  respectively.  Comparing  the  scales  in  Figures  4.11  -  4.13, 1 
see  that  is  relatively  constant,  and  that  most  of  the  change  in  a  with  aperture  radius 

is  due  to  This  is  not  surprising,  since  I  expect  the  change  in  mirror  reflectivity 

to  be  relatively  small  and  the  diffraction  to  increase  as  the  aperture  size  decreases. 
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Figure  4.11  Radial  mode  loss  verses  oxide  aperture  radius  for  the  fundamental  lasing  mode. 
The  lines  are  cubic  spline  fits  of  the  discrete  calculation  data. 
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Figure  4.12  Axial  mode  loss  verses  oxide  aperture  radius  for  the  fundamental  lasing  mode. 
The  lines  are  cubic  spline  fits  of  the  discrete  calculation  data. 
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Radial  Loss  Percentage 


Figure  4.13  Percentage  of  the  total  mode  loss  due  to  radial  losses  for  the  fundamental 
lasing  mode.  The  lines  are  cubic  spline  fits  of  the  discrete  calculation  data. 
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4-5  Diffraction  Mechanisms 

By  examining  the  loss  characteristics  for  each  VCSEL  design,  I  may  deduce  the  physical 
mechanisms  governing  diffraction.  My  results  support  the  idea  that  diffraction  may  be  viewed 
as  a  coupling  between  the  bound  eigenmode  and  the  continuum  of  parasitic  modes  [19]. 
These  are  radially  propagating  slab  modes  in  the  unapertured  (cladding)  region.  Deppe 
explains  in  [21]  that  for  a  cavity  bound  by  perfectly  conducting  mirrors,  the  parasitic  mode 
density  will  resemble  a  slanted  staircase  following  the  three  dimensional  density  of  modes. 
The  jumps  in  the  staircase  occur  at  each  vertical  resonance  in  the  cladding  region.  Deppe’s 
results  suggest  that  shorter  cavities  are  superior  since  the  eigenmode  overlaps  with  a  smaller 
density  of  parasitic  modes.  This  idea  is  illustrated  in  Figure  4.14. 

For  my  high  contrast  DBRs,  I  expect  very  similar  behavior  to  this  ideal  mirror  case. 
Comparing  the  threshold  gains  of  each  structure  (Figure  4.9),  the  A/2  cavity  threshold  is 
consistently  lower  than  all  the  lA  cases.  Moreover,  comparing  the  percentage  of  the  total 
mode  energy  lost  to  diffraction  (Figure  4.13),  the  A/2  VCSEL  is  again  superior.  These  results 
can’t  be  attributed  to  transverse  confinement,  since  the  A/2  VCSEL  has  the  smallest  T*”  of 
all  four  cases.  However,  these  results  could  be  attributed  to  the  larger  total  confinement 
factor  (Figure  4.6),  or  to  a  smaller  density  of  parasitic  modes  interacting  with  the  A/2 
cavity  eigenmode.  It  turns  out  that  these  two  factors  are  closely  related:  the  longitudinal 
confinement  factor — or  effective  cavity  length  in  the  cladding  region — determines  the  location 
of  the  steps  in  the  parasitic  mode  density  (Figure  4.14). 

Adopting  Deppe’s  parasitic  mode  density  interpretation,  I  can  explain  the  difference  in 
threshold  between  the  A/2  and  lA  cavity  VCSELs,  but  I  cannot  easily  explain  the  disparity 
among  the  three  lA  cases,  which  should  have  very  similar  parasitic  mode  densities.  Further¬ 
more,  the  parasitic  mode  density  is  solely  a  function  of  (effective)  cavity  length,  and  therefore 
cannot  address  radius- dependent  changes  in  diffraction.  To  capture  these  radius-dependent 
effects,  and  to  distinguish  the  various  lA  cavity  structures,  I  propose  that  the  diffractive  loss 
is  a  function  of  both  the  density  of  parasitic  modes  and  the  coupling  strength  between  the 
eigenmode  and  the  parasitic  modes.  Moreover,  this  coupling  strength  is  a  function  of  two 
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Figure  4.14  Illustration  of  the  “slanted  staircase”  parasitic  mode  density.  The  density 
follows  the  three-dimensional,  free-space  density  of  optical  modes  (dotted  line). 
The  jumps  in  the  staircase  occur  at  each  vertical  resonance  in  the  cladding 
region.  The  two  discrete  modes  are  representative  of  a  low  threshold  (A),  and 
a  higher  threshold  (B)  lasing  mode. 
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factors:  (1)  the  relative  alignment  of  the  eigenmode  and  parasitic  mode  propagation  vectors, 
and  (2)  the  eigenmode  penetration  into  the  cladding  region. 

The  lA-lTHIN  VCSEL  exhibits  the  lowest  threshold  and  the  lowest  diffraction  loss 
(Figure  4.11)  of  all  three  lA  VCSELs.  Because  of  the  weak  transverse  confinement  (Fig¬ 
ure  4.5)  and  the  resulting  mode  spread,  the  field  in  this  structure  is  more  planar  and  the 
propagation  vector  is  paraxial  with  respect  to  normal  to  the  planar  interfaces.  As  a  result, 
the  eigenmode  wavevector  is  nearly  orthogonal  to  the  parasitic  mode  wavevectors,  which  lie 
principally  in  the  plane  of  the  layer  structure.  The  parasitic  mode  wavevectors  can  have 
components  out  of  the  plane,  but  the  larger  these  components  are  the  less  energy  they  will 
carry  away  in  the  radial  direction.  I  attribute  the  low  threshold  for  this  design  to  this 
misalignment  between  the  eigenmode  and  parasitic  mode  wavevectors.  This  hypothesis  is 
further  supported  by  the  peculiar  results  for  the  double- apertured  1A-2THICK  VCSEL. 

The  double-apertured  VCSEL  has  the  largest  transverse  confinement  factor  of  all  three 
lA  cavity  VCSELs.  Using  the  weighted  index  method  [50],  I  showed  that  this  transverse 
confinement  has  two  primary  effects: 

1.  It  confines  the  mode  energy  within  the  aperture  region. 

2.  It  causes  the  eigenmode  propagation  vector  to  tilt  away  from  normal  to  the  DBRs. 

The  first  effect  acts  to  decrease  the  parasitic  mode  (or  diffractive)  loss  by  containing 
the  mode  energy,  while  the  second  acts  to  increase  it  through  stronger  coupling  (better  prop¬ 
agation  vector  alignment)  to  the  parasitic  modes.  These  processes  compete,  and,  depending 
on  VCSEL  design  and  aperture  radius,  either  effect  may  dominate.  For  the  1A-2THICK 
VCSEL,  the  transverse  confinement  is  >  0.95  for  aperture  radii  from  1.0  to  0.6  )«m;  for  these 
radii,  the  first  effect  dominates.  However,  somewhere  between  0.6  and  0.4  //m  enough  mode 
energy  exists  outside  the  oxide  aperture  such  that  the  strong  parasitic  mode  coupling  causes 
the  optical  loss  to  rapidly  increase.  In  essence,  the  cavity  in  the  unapertured  region  {p  >  p^^) 
forms  a  waveguide  whose  source  is  the  lasing  eigenmode.  When  the  eigenmode  penetration 
into  the  waveguide  region  becomes  large  enough,  the  waveguide  appears  to  “siphon”  energy 
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away  from  the  eigenmode.  This  effect  is  clearly  illustrated  in  the  plots  of  \E^\  for  the  0.6 
and  0.4  fim  cases  in  Figures  4.15  and  4.16,  respectively. 

In  this  chapter,  I  have  ascertained  the  physical  mechanisms  governing  diffraction  in 
oxide  apertured  VCSELs.  The  diffraction  is  well  described  by  the  parasitic  mode  picture.  In 
the  next  chapter,  I  augment  my  weighted  index  method  to  perturbatively  include  coupling 
to  the  parasitic  modes.  The  result  is  a  fast,  semi-analytic  method  for  estimating  all  the 
relevant  optical  laser  parameters  for  microcavity  VCSELs. 
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Figure  4.15  \E^\  for  the  1A-2THICK  structure  with  oxide  aperture  radius  =  0.6  pm. 

The  background  intensity  has  been  set  to  white  to  increase  contrast. 
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Figure  4.16  \E^\  for  the  1A-2THICK  structure  with  oxide  aperture  radius  =  0.4  pm. 

Note  the  field  leakage  into  the  waveguide  formed  by  the  oxide  apertures  and 
the  cavity.  The  background  intensity  has  been  set  to  white  to  increase  con¬ 
trast. 


V.  Diffractive  or  Parasitic  Mode  Loss 

“Things  should  be  made  as  simple  as  possible,  but  not  any  simpler.” 

. . .  Albert  Einstein 

The  weighted  index  method  applied  to  the  vector  Helmholtz  equation  is  a  powerful 
tool  for  analyzing  the  optical  properties  of  azimuthally  symmetric,  but  otherwise  arbitrarily 
complex  VCSELs.  From  this  simple  variational  model  I  was  able  to  predict  the  lasing  mode 
blueshift,  confinement  factors,  spatial  field  profile,  and  mirror  and  absorption  losses.  This  led 
to  the  discovery  of  a  new  phenomenon,  a  polarization  dependence  of  the  mirror  loss  due  to 
the  change  in  the  DBR  reflectance  with  aperture  radius  [52].  This  effect  is  well  documented 
in  edge  emitters  but  had  not  previously  been  predicted  or  observed  in  VCSELs.  The  reason 
it  hasn’t  been  observed  results  from  the  fact  that — ^for  the  current  VCSEL  state-of-the-art — 
diffraction  is  the  dominant  loss  mechanism.  It  is  orders  of  magnitude  larger  than  the  change 
in  mirror  loss  in  the  regime  where  the  DBR  polarization  dependence  is  significant. 

In  the  last  chapter,  I  was  able  to  ascertain  the  fundamental  physical  parameters  de¬ 
termining  diffraction  loss  in  apertured  VCSELs.  The  results  suggested  that  the  parasitic 
mode  concept  proposed  by  Deppe  [21]  accurately  depicts  the  diffractive  loss.  In  this  chapter 
I  present  a  new  perturbative  technique  for  incorporating  the  parasitic  modes  and  calculating 
diffractive  loss  within  the  WIM  formalism. 

Since  the  change  in  reflection  with  radius  is  the  only  polarization  dependent  effect 
observed  for  the  vector  WIM,  and  since  the  resulting  change  in  mirror  loss  is  much  smaller 
than  the  diffraction  loss,  I  simplify  my  WIM  implementation  to  a  scalar  solution  using  the 
well  established  lineax  polarization  (LP)  approximation.  LP  modes  are  much  simpler  in 
form  and  yield  a  much  more  robust  and  rapid  calculation:  roughly  a  factor  of  two  faster 
than  the  full  vector  WIM  for  HEM  modes.  This  increased  efficiency  comes  with  virtually 
no  degradation  in  accuracy.  Furthermore,  the  simpler  form  of  the  LP  modes  facilitates 
my  diffraction  calculation,  which  requires  rigorous  coupling  of  the  lasing  modes  to  a  set  of 
resonant  parasitic  modes. 
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5.1  Linear  Polarization  Approximation  for  the  Lasing  Eigenmodes 

To  solve  for  a  hybrid  (m  ^  0)  mode  in  the  full  vector  WIM,  I  must  iteratively  solve  a 
2x2  TM,  a  2  X  2  TE,  and  a  4  x  4  radial  transfer  matrix  equation.  Since  mode  convergence 
requires  convergence  of  the  resonant  wavelength  and  threshold  gain  for  all  three  equations, 
the  HEM  mode  calculation  is  significantly  slower  (by  roughly  an  order  of  magnitude)  than 
pure  TE  and  TM  mode  calculations.  Furthermore,  despite  the  fact  and  are  separable, 
the  HEM  field  components  are  not,  requiring  a  non-trivial  sum  of  Az  and  Fz  contributions. 
Based  on  the  importance  of  the  HEM  modes — the  HEMu  mode  is  the  fundamental  lasing 
mode — it  would  be  convenient  if  a  simpler  technique  could  be  derived  for  calculation  of  these 
modes.  I  do  this  using  the  linearly  polarized  (LP)  mode  approximation,  originally  applied 
to  the  weakly  guiding  fiber. 

The  LP  modes  are  an  approximation  to  the  full  waveguide  modes  in  the  limit  of 
vanishing  difference  between  core  and  cladding  index  of  refraction.  Equivalently,  the  LP 
modes  are  the  scalar  modes  of  the  cylindrical  waveguide.  Most  often,  the  LP  modes  are 
derived  by  assuming  within  the  fiber-optic  eigenmode  equation  [58] 


,  KM' 

uJn(u)  wKn{w)_ 


I  KH 
uJn  (u)  wKn  (w) 


gclad  ^2 


(5.1) 


where  u  =  fecoreO?  ^  =  ^cladO^  and  a  is  the  fiber  radius.  This  assumption  simplifies  (5.1)  to 


the  LP  eigenmode  equation 


Jm  (w)  Km  (^) 


(6.2) 


which  is  equivalent  to  enforcing  continuity  of  the  scalar  field  and  its  first  derivative  at  p  =  a. 
The  azimuthal  index  m  ^  n  because  the  LP  modes  describe  the  Cartesian/scalar  field  com¬ 
ponents,  while  the  exact  fiber-optic  modes  are  given  in  cylindrical  components.  Furthermore, 
within  the  basic^  LP  mode  approximation  Ez  =  Hz  =  0. 


^As  is  the  case  with  most  perturbational  techniques,  several  variations  of  the  LP  mode  approximation 
exist.  Some  authors  include  z  field  components,  by  either  estimating  them  from  the  in-plane  Cartesian 
components  [24]  or  by  deriving  the  LP  modes  using  a  power  series  in  A  «  {ncore  —  Uciad)  lucore  and  keeping 
higher  order  terms  [70].  In  any  case,  the  basic  approximation  is  the  most  popular  and  is  quite  adequate  in 
most  situations. 
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LP-mode 

Constituent  Fiber  Modes 

Degeneracy 

LPoi 

HEMii 

2 

LPii 

TEqi,  TMoi,  HEM21 

4 

LP21 

HEMii,  HEM31 

4 

LP02 

HEM12 

2 

LP31 

HEM21,  HEM41 

4 

LPx2 

TE02,  TM02,  HEM22 

4 

Table  5.1  First  six  lowest  order  LP  modes,  their  constituent  cylindrical  modes,  and  their 
total  (polarization  and  azimuthal)  degeneracy 

The  LP  modes  are  waveguide  modes  with  no  polarization  effects  included.  Therefore, 
when  the  LP  mode  approximation  is  valid  the  constituent  fiber  modes  are  nearly  degener¬ 
ate,  sharing  the  same  propagation  constant  and  cross-sectional  intensity  distribution.  For 
example,  the  LPn  mode  consists  of  the  doubly  degenerate  HEM21  mode,  and  the  TEqi  and 
TMoi  modes,  all  characterized  by  an  intensity  pattern  with  a  single  null  at  the  center.  Table 
5.1  lists  the  six  lowest  order  LP  modes,  their  corresponding  cylindrical  (fiber)  modes,  and 
their  total  degeneracy  due  to  polarization  and  azimuthal  field  dependency:  either  sin  (m^) 
or  cos  {m(f>)  [58].  For  VCSELs,  I  expect  the  same  criteria  will  apply,  and  the  LP  approxima¬ 
tion  should  be  valid  for  nearly  degenerate  modes.  The  vector  WIM  LPn  resonance  results 
(Figure  3.13)  show  very  little  splitting  throughout  the  convergent  range  of  the  calculation, 
suggesting  the  validity  of  the  LP  approximation  for  this  case.  This  conclusion  is  further 
supported  by  the  DBR  reflectance  calculation  which  displayed  a  relatively  small  difference 
between  modes  over  the  entire  span  of  radii  investigated.  Lastly,  I  note  that  the  LP  approx¬ 
imation  is  equivalent  to  a  paraxial  approximation,  suggesting  it  should  be  most  valid  for  the 
lowest  order  WIM  modes.  Overall,  I  suspect  the  error  due  to  the  LP  approximation  will  be 
no  worse  than  the  error  introduced  by  the  fundamental  separability  assumption  upon  which 
the  WIM  is  based. 

Proceeding  with  the  “standard”  LP  modes,  I  write  the  lasing  mode  fields  as 

E  w  ExX  =  F  (p)  G  {z)  cos  {mcf))  X ,  (5.3) 

H  «  HyV  =  —F{p)G' {z)  cos  {m(f>)y.  (5.4) 

U}P,Q 
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Here 


GW 


Cl  J„  (kip)  p<p 

OX 

C2Km  {kip)  p>  p  ox 


(5.5) 

(5.6) 


and  I  have  limited  the  geometry  to  two  radial  regions:  core  (p  <  p^^)  and  cladding  {p  >  Pox)- 
To  solve  for  modes,  I  enforce  continuity  of  the  field  and  its  first  derivative  at  each  axial  and 
radial  interface.  The  remaining  analysis  parallels  the  vector  TE  mode  solution  presented  in 
Chapter  II,  with  the  exception  that  I  now  solve  directly  for  the  fields  rather  than  the  vector 
potentials. 

For  the  axial  problem,  I  find  roots  of  the  TE  threshold  condition  (2.43).  For  each 
resonant  mode,  the  Pj  are  found  from  the  axial  dispersion  relation 


Here  {er-.F,j)  and  (k)  are  the  weighted  relative  permittivity  and  radial  propagation  constants, 
respectively  defined  by 


{F\^rUj\F) 

{F\F) 


and 


w  = 


{F\F) 


(5.8) 


Since  there  can  be  no  confusion  with  the  polarization  of  the  weighting  functions  for  LP 
modes,  no  a  superscript  is  required.  It  is  then  a  simple  matter  to  back  solve  the  individual 
transfer  matrices  (2.37)  to  find  the  aj  and  bj.  For  the  radial  problem  I  find  the  roots  of  the 
LP  eigenmode  equation  (5.2).  This  is  precisely  the  same  as  solving  the  radial  TE  threshold 
condition  (2.60)  for  two  regions.  The  core  and  cladding  propagation  constants  are  given  by 
the  radial  dispersion  relation 

hi  =  -  (fif.  (5.9) 
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Here  (erri.G)  and  (^)  are  defined  as 


{^r:i,G) 


{G\ei,s\G) 

{G\G) 


and 


(G|G) 


(5.10) 


respectively.  For  each  mode,  the  coefficients  Ci  and  C2  may  be  found  by  the  field  continuity 
condition.  The  axial  and  radial  problems  are  alternately  solved,  iterating  until  both  the 
resonant  wavelength  and  material  gain  converge.  The  result  of  this  analysis  is  a  set  of  LP 
lasing  modes,  all  bound  in  the  radial  direction  and  radiative  in  the  axial  direction.  In  order 
to  capture  diffractive  effects  I  must  augment  these  modes.  I  accomplish  this  by  coupling 
each  lasing  mode  to  a  set  of  resonant  parasitic  modes. 


5.2  Parasitic  Modes 

To  calculate  lasing  modes  with  the  weighted  index  method,  I  enforce  radiation  condi¬ 
tions  in  the  axial  direction  and  evanescent  decay  in  the  radial  direction  to  form  constraint 
equations  on  the  mode  frequency  and  material  gain.  To  calculate  parasitic  modes  with 
the  weighted  index  method,  I  relax  the  evanescent  condition  and  allow  for  (outward)  radial 
propagation.  Alternatively,  I  may  derive  the  parasitic  modes  as  exact  field  solutions  in  the 
cladding  region  (p  >  p^^).  As  I  will  show  below,  these  two  versions  of  the  parasitic  modes 
are  equivalent. 

Since  parasitic  modes  may  propagate  at  any  angle  with  respect  to  the  planar  interface 
normal  (1),  I  do  not  make  the  LP  approximation  and  instead  solve  explicitly  for  TE  and 
TM  cylindrical  modes.  By  dividing  the  modes  into  TE  and  TM,  they  may  be  conveniently 
described  by  a  single  longitudinal  field  component:  Hz  for  TE  modes  and  Ez  for  TM  modes. 
For  HEM  lasing  modes,  the  relative  strength  of  the  TE  and  TM  components  was  important 
in  order  to  match  the  tangential  field  components  across  radial  boundaries.  But,  for  the 
parasitic  modes  there  is  no  such  requirement:  I  am  only  concerned  with  their  behavior 
in  the  cladding  region.  Furthermore,  the  TE  and  TM  parasitic  modes  each  obey  the  axial 
boundary  conditions  on  their  own,  implying  they  may  be  dealt  with  independently  regardless 
of  the  azimuthal  mode  order. 
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Based  on  the  azimuthal  form  of  the  LP  modes,  I  choose^ 


If™  =  P(p)Q{z)  sin  (5.11) 

E™  =  P(p)Q  (z)  cos  (n^.),  (5.12) 

where 


PM  = 

ClJnihp)  P<Pox 

C2HI  {k2p)  P  >  Pox 

(6.13) 

QM  = 

(5.14) 

Since  the  TE  and  TM  modes  are  dealt  with  independently,  no  polarization  superscript  is 
required,  and  I  am  free  to  use  the  same  functions  (P,Q)  for  and  if™.  The  choice  of 
Ef^  ~  cos  (n^)  and  ~  sin(n<^)  is  necessary,  due  to  azimuthal  orthogonality,  for  the 

lasing  mode  to  couple  with  the  parasitic  modes.  For  each  lasing  mode,  I  calculate  a  set  of 
resonant  parasitic  modes  following  a  procedure  analogous  to  the  first  axial  iteration  of  the 
TE/TM  vector  mode  solution  given  in  Section  2.4. 

To  solve  the  axial  problem,  I  must  first  calculate  a  set  of  weighted  material  parameters 
{^r:P,j)i  found  via 

M  ^  (5-15) 

Since  the  parasitic  modes  are  radiative  in  the  radial  direction,  the  inner  product  over  p 
(2.19), 

poo 

{P\P)=  \P{p)fpdp,  (5.16) 

Jo 

will  not  converge.  However,  employing  an  appropriate  limiting  procedure,  (5.15)  will  con¬ 
verge^  and  gives  the  material  values  in  the  cladding  region, 

^Note,  to  avoid  excessive  sub/superscripting,  I  have  recycled  the  P  and  Q  functions,  redefining  them  from 
their  original  designation  in  Chapter  II. 
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Similarly,  the  weighted  radial  propagation  constant  will  equal  the  cladding  value, 


{k) 


{P\kUP) 

{P\P) 


=  ^2- 


(5.18) 


This  implies  that  the  parasitic  mode  solution  obtained  via  the  WIM  procedure  is  ■precisely 
the  same  as  the  exact  cylindrical  mode  solution  in  the  cladding  region.  The  dispersion 
relation  for  the  axial  problem  is 


(5.19) 


Since  I  am  looking  for  parasitic  modes  resonant  with  a  particular  lasing  mode,  I  set 
oj  =  wiase  and  solve  for  the  complex  ^2  roots  of  the  axial  threshold  condition 


^22(^2)  —  O'  (5.20) 

Equation  (5.20)  is  the  same  as  the  TE/TM  threshold  condition  given  in  Chapter  II  (2.43), 
the  only  difference  is  the  argument  for  which  the  roots  are  found.  I  allow  k2  to  be  complex, 
rather  than  solve  for  a  real  and  Kactivo  since  material  gain  does  not  apply  to  the  parasitic 
modes.  The  source  of  the  parasitic  modes  is  the  lasing  mode.  Furthermore,  the  complex  k2 
term  compensates  for  the  axial  radiation  of  the  parasitic  modes.  As  I  will  show  in  Section 
5.4,  the  Im(jb2)  has  no  impact  on  the  energy  carried  away  from  the  lasing  mode.  Since  this 
energy  “siphoning”  is  the  only  effect  I  am  after,  I  allow  k2  to  be  complex  as  a  matter  of 
practical  convenience  in  my  root  search  algorithm,  ignoring  the  Im(A;2)  for  all  subsequent 
calculations. 

Once  (5.20)  is  solved,  I  could  continue  with  the  WIM  procedure  and  solve  the  radial 
problem,  determining  weighted  (/?)  and  {er:i,Q),  and  finding  the  ki  and  ci  that  satisfy  the 
radial  TE/TM  (or  HEM)  boundary  conditions.  However,  these  results  will  have  no  bearing 
on  Q  (z),  implying  the  axial  solution  is  uncoupled  from  the  radial  solution.  Moreover,  these 
results  will  have  no  effect  on  P{p)  for  p  >  p^^,  since  k2  is  completely  determined  by  the 
axial  solution  (5.20)  and  the  radial  boundary  conditions  only  give  ci  in  terms  of  C2  (i.e.,  one 
degree  of  freedom  exists  in  the  radial  solution,  typically  this  is  used  to  normalize  P{p))-  As 
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I  will  demonstrate  in  Section  5.3,  only  the  cladding  region  behavior  is  required  to  calculate 
the  diffractive  loss.  Therefore,  for  each  parasitic  mode  I  only  solve  (5.20)  and  ignore  P{p) 
for  p  <  poj.,  substantiating  my  earlier  claim  that  TE  and  TM  modes  may  be  solved  for 
independently.  Based  on  this  discussion,  and  to  facilitate  the  following  analysis,  I  simplify 
the  parasitic  mode  definition  to 

^PM  _  H^\kp)Q{z)sm{n<t>),  (5.21) 

(kp)  Q  {z)  cos  {n<f)) .  (5.22) 

Here  k  =  k2  from  (5.13),  and  the  fact  that  these  only  apply  for  the  cladding  region  is  im¬ 
plicitly  understood.  An  illustration  of  the  standing  wave  intensity  profile  for  a  TM  parasitic 
mode  is  given  in  Figure  5.1. 

5.3  Lasing  Mode-to-Parasitic  Mode  Coupling 

The  WIM  lasing  modes  are  not  eigenmodes  of  the  non-separable  VCSEL  cavity,  but 
are  instead  a  product  of  (one-dimensional)  eigensolutions  of  the  weighted  axial  and  radial 
problems.  One  could  envision  a  separable  cavity  with  Sr(p,2:)  =  er{p)  +  £r(z)  for  which 
the  WIM  lasing  mode  was  a  true  eigenmode.  In  this  case,  each  WIM  lasing  mode  would  be 
independent  and  exact,  and  there  would  be  no  mechanism  for  coupling  to  the  parasitic  modes. 
This  may  in  fact  be  a  viable  new  design  approach.  Theoretically,  if  such  a  device  could  be 
fabricated,  the  (WIM)  lasing  modes  would  suffer  no  diffraction  loss.  An  illustration  of  this 
type  of  device  is  shown  in  Figure  5.2.  The  diffraction  or  parasitic  mode  loss  results  from  the 
non-separability  of  the  VCSEL  structure^,  which  effectively  couples  the  lasing  mode  to  the 
parasitic  modes.  To  calculate  this  coupling,  I  employ  two  basic  electromagnetic  theorems: 
surface  induction  and  reciprocity. 

®In  analogy  to  quantum  optics  theory,  the  WIM  lasing  mode  and  the  resonant  parasitic  modes  are  “bare” 
field  modes.  The  true  three-dimensional,  lossy  cavity  modes  are  “dressed”  field  modes.  The  coupling  between 
the  WIM  lasing  mode  and  the  parasitic  modes  determines  the  weighting  coelficients  of  the  various  bare  modes 
in  the  representation  of  a  single  dressed  mode. 
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Intens 
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Figure  5.2  Illustration  of  an  azimuthally  symmetric,  separable  VCSEL  geometry.  The 
change  in  relative  permittivity  as  a  function  of  z  produces  the  Fabry-Perot 
cavity  structure.  The  change  in  relative  permittivity  as  a  function  of  p  defines 
the  lasing  modes  analogous  to  fiber-optic  waveguide  modes.  If  such  a  VCSEL 
could  be  fabricated,  theoretically  the  lasing  modes  would  suffer  no  diffraction 
loss. 
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5.3.1  Surface  Induction.  For  the  first  step  in  my  coupling  calculation,  I  derive 
fictitious  electric  and  magnetic  surface  currents  in  the  cladding  region  {p  > 

pg^).  These  currents  replace  the  eigenmode  in  the  cladding  region  and  act  as  a  source  for 
the  parasitic  modes.  The  sole  purpose  for  performing  this  substitution  is  to  mathematically 
isolate  the  parasitic  mode  fields  from  the  eigenmode  fields,  making  it  much  easier  to  solve 
for  the  coupling  strength  to  each  parasitic  mode.  It  is  crucial  to  realize  that  all  I  am 
doing  is  restating  the  field  problem  in  a  different,  but  equivalent,  form.  The  combination 
of  the  currents  and  the  parasitic  modes  in  the  cladding  region  is  exactly  equivalent  to  the 
real  electromagnetic  problem  involving  the  eigenmode  and  parasitic  modes.  The  surface 
induction  theorem^  [7]  describes  how  to  represent  the  fields  with  equivalent  currents. 

To  derive  equivalent  currents,  I  assume  the  total  field  in  the  core  region  [p  <  p^j.)  is 
accurately  represented  by  the  WIM  lasing  mode, 

<l>,  z)  «  E™{p  <  <^,  z),  (5.23) 

z)  «  H™{p  <  <t>,  z).  (5.24) 

In  the  cladding  region,  I  assume  the  total  field  may  be  written  as  the  sum  of  the  lasing  mode 
and  the  resonant  parasitic  modes, 

#*'^(p,  <l),  z)  «  E^“^(p  >  p„^,  (j),z)  +  '^  AiEf^ip  >  P<,^,  <j>,  z),  (5.25) 

;=i 

^PM 

«  H™(P>  p„A,z).  (5.26) 

l=\ 

Here  iVpM  is  the  total  number  of  peirasitic  modes  which  couple  to  the  lasing  mode  (includes 
both  TE  and  TM  parasitic  modes)  and  A/  is  the  amplitude  of  the  /’th  parasitic  mode. 
These  assumptions  are  required  for  a  “clean”  electromagnetic  problem  but  are  otherwise 
not  directly  related  to  surface  induction.  These  fields  constitute  the  “actual  problem,” 

^The  inductance  theorem  is  closely  tied  to  the  uniqueness  theorem  which  states  that  the  fields  in  a  closed 
region  are  uniquely  specified  by  the  sources  with  the  region  plus  the  tangential  field  components  on  the 
boundary. 
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illustrated  in  Figure  5.3.  In  the  actual  problem,  no  sources  exist  and  the  total  tangential 
fields  are  continuous  across  the  p  =  p^^  interface. 

To  apply  surface  induction,  I  consider  the  “closed”  surface  Fdad  which  encompasses 
the  cladding  region  volume®, 


ficiad  =  [pox^  oo]  X  [0» 27r]  X  [-00,  oo].  (5.27) 

Within  Ociad,  I  replace  the  WIM  eigenmode  fields  with  the  equivalent  electric  and  magnetic 
currents.  Therefore,  in  the  “surface  inductance  equivalent  problem,”  illustrated  in  Figure 
5.3,  the  total  fields  in  the  cladding  region  are  given  by 

(5.28) 

/=1 

^PM 

PoxA,^)-  (5-29) 

1=1 

Here  SI  denotes  surface  inductance  equivalent  fields. 

In  the  actual  problem,  no  currents  exist  and  the  fields  obey  the  tangential  field  con¬ 
tinuity  relations.  In  the  surface  induction  equivalent  problem,  the  equivalent  currents  are 
required  to  enforce  continuity  of  the  tangential  fields.  The  electric  and  magnetic  currents  are 
defined  by 


=  n  X  ,  (5.30) 

^  ^^core  _  ^clad^  ^  (5  3I) 

where  n  is  the  outward  normal  to  Fdad-  The  difference  in  form  between  (5.30)  and  (5.31) 
is  strictly  due  to  convention.  Since  the  lasing  mode  is  evanescent  in  the  cladding  region, 
the  equivalent  currents  are  only  non-zero  along  the  boundary®  p  =  p^^.  Hence,  the  outward 

®Here  I  have  assumed  that  Fdad  may  be  closed  at  infinity,  a  common  practice  within  theoretical 
electromagnetics . 

^Consistent  with  the  general  WIM  procedure,  I  assume  the  currents  go  to  zero  for  z  >  ztop  and  z  <  Zbottom, 
the  upper  and  lower  VCSEL  boundary  planes 
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normal  n  =  —p  and  the  equivalent  surface  currents  are 


^PM 


= 


= 


(-?)  X  AiSrip  >p„A, 

I  /=i 

{^PM  ^ 

^^{p:..<i>.z)  -  £  AiEj^ip  >  p,,A,z)\  X  {-p) 
=  px  E'^^^{p'f^,<f>,z). 


(5.32) 


(5.33) 


Here  and  p^^.  lie  just  outside  and  inside  of  Fdad,  respectively. 

Note  that  I  could  have  represented  the  WIM  eigenmode  in  the  cladding  region  using 
equivalent  volume  currents.  In  this  approach,  the  equivalent  currents  arise  from  the  difference 
in  the  weighted  material  parameters — for  which  the  WIM  solution  is  exact — and  the  actual 
material  parameters.  Unlike  the  surface  induction  solution,  there  is  no  magnetic  volume 
current,  since  p  =  Pq  everywhere.  However,  the  electric  current  has  both  a  p  and  a  ^ 
component,  and  therefore  couples  to  both  the  TE  and  TM  parasitic  modes.  This  approach  is 
attractive,  since  it  allows  for  an  easier  interpretation  of  the  source  of  the  parasitic  mode  loss; 
the  difference  between  the  weighted  and  real  geometries.  I  chose  to  use  equivalent  surface 
currents,  since  they  allow  for  an  easier  implementation  and  yield  the  same  diffraction  loss 
results. 

The  surface  induction  picture  is  now  complete:  inside  fidad  the  total  fields  are  given  by 
the  weighted  sum  of  the  parasitic  modes  ((5.28)  and  (5.29)),  whose  source  is  the  magnetic 
and  electric  equivalent  currents  ((5.32)  and  (5.33)).  The  surface  inductance  problem,  which 
is  equivalent  to  the  actual  problem  in  the  cladding  region,  is  illustrated  in  Figure  5.3.  All 
that  remains  is  to  couple  the  currents  to  the  parasitic  modes  to  find  the  amplitude  (A/) — or 
coupling  strength — of  each  parasitic  mode.  For  this,  I  use  reciprocity. 

5.S.2  Reciprocity.  The  reciprocity  theorem  is  a  simple  rearrangement  of  the 
Maxwell  curl  equations  using  appropriate  vector  identities.  It  is  commonly  used  in  waveg¬ 
uides  problems  to  calculate  the  coupling  between  two  guides  or  the  radiation  from  a  source 
into  a  guide.  I  employ  reciprocity  in  a  similar  capacity,  to  compute  the  At  coefficients. 
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Figure  5.3  Illustration  of  the  actual  VCSEL  electromagnetic  problem  (left)  and  the  surface 
inductance  equivalent  problem  (right).  In  the  surface  inductance  problem  the 
WIM  fields  in  the  cladding  region,  denoted  the  dashed  box,  are  replaced  by 
equivalent  surface  currents  J®®  (5.32)  and  M®®  (5.33).  The  total  fields  in  the 
cladding  region  in  the  equivalent  problem,  and  are  given  by  the 

weighted  sum  of  parasitic  modes  (5.28)  and  (5.29). 
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Reciprocity  dictates  that  any  two  fields  in  the  volume  ficiad  must  satisfy 


JJ  X  -  E’’  X  ■  ds  =  jJJ  1^'’  •  J“  -  •  M“}  dV. 

Fclad  ^clad 

(5.34) 

For  my  application,  I  assign  the  a  fields  and  sources  to  the  surface  induction  cladding  fields 
((5.28)  and  (5.29))  and  equivalent  currents  ((5.32)  and  (5.33)).  For  the  b  fields,  I  use  a 
fictitious  “test”  mode^:  a  single  parasitic  mode  propagating  radially  inward.  The  test  mode 
is  a  non-physical — since  it  propagates  inward — parasitic  mode,  which  otherwise  obeys  all  of 
Maxwell’s  equations  and  the  parasitic  mode  orthogonality  relation  (5.49).  Since  the  strength 
of  the  parasitic  modes  is  completely  determined  by  the  equivalent  currents  (and  therefore, 
the  eigenmode),  I  solve  for  TE  and  TM  parasitic  mode  Ai  separately,  even  for  m  ^  0  modes. 

First,  I  consider  the  TE  parasitic  modes.  I  define  the  test  mode  or  b  fields  (5.34)  by 


{k(p)  Qi  {z)  sin  {n<f) . 


(5.35) 


h;  is  an  exact  duplicate  of  one  of  the  real  parasitic  modes,  with  the  exception  that  it  travels 
in  the  opposite  direction.  To  apply  reciprocity,  I  assign  [60] 


ATpM 


_  £iclad  _ 

=  ^si  -- 

Y, 

(5.36) 

/=1 

•^PM 

ga 

_  frclad  _ 

=  ^si  — 

EM™, 

(5.37) 

/=1 

= 

(5.38) 

= 

(5.39) 

=  E-, 

(5.40) 

=  H-, 

(5.41) 

=  0, 

(5.42) 

M” 

=  0. 

(5.43) 

^The  use  of  a  “test”  mode  is  standard  practice  in  waveguide  solutions;  see,  for  example,  [60]. 
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Beginning  with  the  surface  integral  in  (5.34),  I  note  that 


ds  =  {-p)po^d<l)dz 


(5.44) 


on  p  =  Pax-  Substituting  in  ds  and  performing  the  cross  and  dot  products,  the  surface 
integral  reduces  to 

/  /  {ElH;-ElHl]  p^d^dz,  (5.45) 

assuming  the  surface  region  [0,2x]  x  [^bottom, -^top]  is  the  only  region  contributing  to  the 
integral®.  Calculating  and  via  Ampere’s  law,  (shown  here  tailored  for  TE  modes) 


yields 


(5.46) 


El  =  '  (k(p)  Q(  {z)  sm  (n^) , 

(5.47) 

^PM 

u°H^^'{hp)Qi{z)sin{n^), 

/=i  ' 

(5.48) 

where  '  denotes  differentiation  with  respect  to  the  argument.  Before  proceeding,  I  make  the 
crucial  assumption  that  the  parasitic  modes  are  orthogonal,  obeying  the  relation 


X  -ds  =0  for  a^h. 

Tclad 


(5.49) 


Equation  (5.49)  sifts  the  /  =  /  mode  out  of  the  summation  (5.48),  reducing  the  surface 
integral  (5.45)  to 

{k,p„)  n  lQ,(z)fdz.  (5.60) 

(  ♦'^bottom 


®This  assumption  may  be  justified  by  including  a  very  small  amount  of  absorption  in  the  cladding  region, 
such  that  the  parasitic  modes  go  to  zero  at  infinity. 
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Here  W  {z)  is  the  Hankel  function  Wronskian  defined  by  [1] 


w  (^)  =  (^)  (")  -  w  (")  = 


(5.51) 


and  the  azimuthal  integral  equals  tt.  Finally,  substituting  in  (5.51)  and  rewriting  the  integral 


over  in  Dirac  notation  gives 


in 


E^xH’’  -E^x  •  ds  =  -Ac 


4u}fio 

kj 


Q{  I  ’ 


(5.52) 


where  the  indicates  the  inner  product  is  unconjugated. 

Turning  now  to  the  right  hand  side  of  (5.34),  I  begin  by  noting  that  the  LP  modes 
have  no  z  field  component.  Therefore,  the  effective  currents  calculated  in  (5.32)  and  (5.33) 
will  be  strictly  z  directed.  Furthermore,  and  are  surface  currents,  non-zero  only  on 
the  boundary  p  =  reducing  the  volume  integral  in  (5.34)  to  a  surface  integral.  Making 
these  simplifications,  the  right  hand  side  of  the  reciprocity  relation  (5.34)  becomes 


r^top  /*27r 


HlMtPo.d<f>dz, 


(5.53) 


where  the  only  coupling  occurs  through  M®®  (conversely,  the  TM  parasitic  modes  will  only 
couple  through  J®«).  Explicitly  calculating  the  magnetic  current  due  to  the  LP  lasing  mode 


gives 


•  (p  X  «)  F{p)G{z)  cos  (m^) 

=  ^F{p)G{z)  [sin  ((m  -  1)  ^)  -  sin  ((m  -|- 1)  4)] 


(5.54) 


Substituting  M®®  and  into  (5.53)  yields 


jjj  {£*  M‘}  dV  =  (k,p„)  (  G I  Q,)  (^) .  (5,55) 
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Here 


{<!>) 


—  1)  <f>)  —  sin  ((m  +  1)  (f>)]  sin  {n<f))  d<f> 


—2%  m  =  0,n  =  m  +  l 

<  —TT  m>l,n  =  m  +  l  » 

X  m>l,  n  =  m  —  1 


(5.56) 


determines  to  which  parasitic  modes  the  lasing  mode  will  couple.  Setting  the  left  hand  side 
(5.52)  equal  to  the  right  hand  side  (5.55),  and  solving  for  the  coupling  constant  gives 

^  -  «„(„«)  (1  +  i^o)]  .  (5.57) 

{Qi\Qi) 

Here  I  have  replaced  {(j))  with  the  corresponding  sum  of  Kronecker  delta  functions  (5a6), 
and  added  the  TE  superscript  on  A(  for  clarity. 

To  calculate  the  TM  parasitic  mode  coupling  coefficient,  I  use  the  same  a  fields  but 
consider  a  radially  inward  propagating  TM  mode  for  the  b  fields  in  (5.34), 


K  = -^^i^n%p)<?/'(^)cos(ra(^).  (5.58) 


Following  an  analogous  procedure,  I  apply  reciprocity  using  the  same  field  assignments  given 
in  (5.36)  -  (5.43).  For  this  case,  the  surface  integral  is  given  by 

p^top  /*27r 

/  /  f,„d,j,dz.  (5.59) 

'^•Sbottom  ''0 

Calculating  and  via  Faraday’s  law,  (shown  here  tailored  for  TM  modes) 


(5.60) 
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yields 


•  '  JT  IVA 

HI  =  -£a, 


HI  = 


H^^'i^lp)  <3/'(2;)cos  (n(^). 


(5.61) 


(5.62) 


Substituting  the  fields  into  (5.59),  and  minding  the  orthogonality  relation,  gives 

ff  n  eAQi{^)fdz 

y^rclad  ^  «^^bottom 

=  A,^(<3,KTi3()(1+M-  (5-63) 


The  TM  volume  integral  couples  the  parasitic  modes  to  the  electric  current. 


r^top  /*27r 
*'  ^bottom  *'0 


Hz^z  Pox 


Explicitly  calculating  the  current  due  to  the  LP  lasing  mode  gives 


=  z  -{y  X  p) - F{p)G'{z)  cos  {mcj)) 

OJPQ 


%jpo 


F{p)G'{z)  [cos  ((m  -!)</')  +  cos  ((m  +  1)  4>)] 


Substituting  J®®  and  into  (5.64)  yields 


where 


(5.64) 


(5.65) 


fjj  it’’  lH^}  iV=  (kKj  {  G'  I  Ql)  *™  W  ’  (5-63) 


p2'K 

(<^)  =  /  [cos  ((m  —  1)  ^)  +  cos  ((m  +  1)  <^)]  cos  (n(^) 

Jo 
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27r  m  =  0,  n  =  m  +  l 

TT  m>l,  n  =  m  +  l 

27r  m  =  l,  n  =  m— 1 
TT  m>l,  n  =  m  —  1 


(5.67) 


Setting  the  left  hand  side  (5.63)  equal  to  the  right  hand  side  (5.66),  replacing  {<i>)  with 
the  corresponding  sum  of  Kronecker  delta  functions,  and  solving  for  the  coupling  constant 
gives  (5.68), 


^  (Qi  \  Sri  Qf)  L  1+^nO 

(5.68) 

The  appearance  of  {kfPox)  (5.57)  and  (5.68)  is  at  first  disconcerting.  One  would 
assume  the  form  of  the  test  field  would  have  no  influence  on  the  coupling.  However,  the 
analysis  is  not  yet  complete.  To  have  any  physical  meaning  Aj^  and  A™  must  be  applied 
to  the  parasitic  modes  and  the  resulting  radiative  loss  computed,  as  I  show  below. 


5.4  Parasitic  Mode  Loss 

The  ultimate  goal  of  the  parasitic  mode  analysis  is  to  determine  the  change  in  threshold 
gain  from  the  WIM  value  due  to  diffraction.  To  do  this,  I  define  diffraction  loss  as  the 
eigenmode  power  coupled  to  the  parasitic  modes  and  radiated  in  the  transverse  direction. 

In  the  WIM  solution,  the  calculated  threshold  gain  compensates  for  the  ab¬ 
sorption  (fjfabsorption)^  mirror  losses. 


^WIM 

9th 


1 

ptot 


|Q,(absorption)  _|_  ^(mirror) 


(5.69) 


where  is  the  total  confinement  factor.  The  total  threshold  gain,  on  the  other  hand, 
incorporates  the  absorption,  mirror,  and  diffraction  (afdiffraction)^  losses. 


^  I  ^(absorption)  ^  ^(mirror)  _|_  ^(diftaction)| 


(5.70) 
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As  shown  in  Section  4.4.5  (4.29),  the  losses  may  be  written  in  terms  of  the  power  dissipated 
^p(absoiption)^  p(inirror)^  p(diffraction)^  ^  stored  field  energy  (We  +  W,„),  and  the  phase 


velocity  (vp)  [7]  as 


F(  source) 


,,  (5.71) 

where  source  =  absorption,  mirror,  or  diffraction.  By  substituting  the  expressions  for 
Q;(absorption)^  Q,(miiTor)^  gi,(<liffraction)  rearranging,  I  may  write  the  total  threshold 

as 

^(diffraction)  \ 


_  /  pldiffraction)  N 

yth  9th  I  *  ^(absorption)  J3(niirror)  J 


(5.72) 


Equation  (5.72)  gives  the  perturbative  adjustment  to  the  WIM  threshold  gain  required  to 
capture  diffractive  loss®.  The  only  remaining  task  is  to  calculate  the  three  dissipated  power 
terms. 


To  find  the  power  terms,  I  follow  a  procedure  similar  to  the  one  presented  in  Section 
4.4.5.  In  this  approach,  I  assume  the  total  threshold  gain  compensates  for  the  total  power 
exiting  the  cylindrical  volume  of  the  VCSEL  core. 


f^core  —  [0? ^  [0,  27r]  X  [-S^bottom? -^^top]  ? 


(5.73) 


either  through  absorption  or  radiation  out  the  ends  or  the  side.  The  power  absorbed  due 
to  material  loss  (for  example,  free  carrier  losses)  results  directly  from  the  imaginary  part  of 
the  field  energy  distribution.  For  a  lossy  region,  is  complex  and  the  total  (time  averaged) 
power  absorbed  from  the  electric  field  is  given  by 


^(absorption) 


-nil 


eo  Im  (cr) 


E 


WIM 


dV, 


(5.74) 


where  the  Im(e:r)  <  0  for  absorption.  There  is  no  corresponding  magnetic  field  term  since 
fi  =  Hq  throughout  the  structure.  Substituting  in  the  explicit  form  for  the  LP  mode  (5.3), 

®In  the  derivation  of  (5.72),  I  have  implicitly  assumed  that  the  parasitic  modes  only  contribute  to 
jj(diffraction)^  and  that  the  absorption  and  mirror  loss  are  well  described  by  the  WIM  results. 
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p(absorption)  ‘g  gjy^n  in  Diiac  notation  as 


p(.wi».)  =  f  I F  )^  { G I  Im  (e,)  I G )  (1  +  «™o) .  (5.75) 


Here 

rpox 

{P\F)^=  \F?pip,  (5.76) 

Jo 

is  the  standard  inner  product  over  p  in  the  core  region,  and  (G|  Im(er)  I  is  found  via 
the  standard  inner  product  over  (2.26).  The  two  radiated  power  terms  are  calculated  from 
the  (time  averaged)  integrated  Poynting  vector, 


(radiation)  _  ^ 

2 


y/r 


ds 


(5.77) 


where  the  *  indicates  the  complex  conjugate.  For  the  mirror  loss,  the  fields  are  the  WIM 
lasing  mode  fields  and  ds  =  iz  p  d<f>  dp:  +  for  the  top  surface,  and  —  for  the  bottom. 
Substituting  in  the  LP  modes  gives 

p^)  ^  --J:_(F|F)^Re[iG(z.,p){G'(^.op)r](l+i„o),  (6.78) 

ZLOIXq 

TftZ'’  =  2;^(J'I^')„„R*  [;G(^b...o«){G'(.^.»«))'](l  +  <S».o),  (5.79) 

and  =  P.*””'*  +  PfiiZ’- 

For  the  diffraction  loss,  the  fields  in  (5.77)  are  the  parasitic  mode  fields  and  ds  = 
ppd<j>dz.  The  unconjugated  parasitic  mode  orthogonality  relation  (5.49)  is  generaP®,  and 
is  equally  valid  in  the  conjugated  form.  As  a  result,  I  neglect  the  cross  terms  and  pull  the 
summation  of  modes  outside  the  integral  in  (5.77)  to  give 


^PM  1  r  ^27r  /  /  \  *\ 

p(diffraction)  =  ^  L  \Aif  Re  /  X  Y  p  p  dz  d<j> 

I— I  ^  l-•'0  Jz}iottom 


^®The  conjugated  form  of  the  orthogonality  relation  holds  for  nonabsorbing  materials.  The  unconjugated 
form  holds  for  nonabsorbing  and  absorbing  materials.  This  implies  the  unconjugated  form  is  a  more  general 
form,  encompassing  the  conjugated  form.  See  Chapter  31  of  [70]  for  a  detailed  discussion  of  the  relationship 
between  the  two  forms. 
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(6.80) 


/=1  /=! 

where  Npu  =  -^pm  +  -^PM  •  ^or  the  TE  parasitic  modes  the  integral  kernel  is  E™ 
Substituting  in  the  TE  parasitic  mode  fields  given  in  (5.21)  and  (5.46)  yields 


pTe  ^  |^TE|2  {nm  (  <?,  |  Ol  >  (1  -  «no)  •  (5.81) 


Substituting  in  (5.57)  for  Aj^  gives 


(5-82) 


128 


X  - - -  .  {^IPox)  [^n(m-l)  ^n(m+l)  (1  +  ^mo)]  • 

{Qi\Qi) 

Here  I  have  used  the  identity  Re  (z)  =  Im  (iz).  Similarly,  for  TM  parasitic  modes  the  integral 
kernel  is  where  is  given  in  (5.22)  and  is  found  via  Faraday’s  law, 


Substituting  the  TM  parasitic  mode  fields  into  (5.80)  yields 


(5.83) 


pTM  ^  J5^  y  ( Q,  I  I >]  (1  +  «„,) .  (5.84) 


Substituting  (5.68)  for  A™  gives 


128  (cu/io)  cjeo  L^/  J 


P{Pox)Hn^  i^lPox)  [^n(m+l)  (1  +  ^mo)  +  ^n(m-l)  (1  +  ^ml)]  • 


(g'IQi) 

Ql\£r\ 
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I  now  have  all  the  elements  required  to  calculate  the  diffraction  loss  contribution  to  the  total 
threshold  gain.  To  gain  a  better  physical  insight  into  the  diffraction  effects,  it  is  useful  to 
first  interpret  the  TE  and  TM  diffraction  power  terms. 

and  P™  may  be  simplified  by  approximating  the  Hankel  functions  with  their 
asymptotic  forms, 

p-iz 

and  (5.86) 

Substituting  into  the  two  TE  terms  gives. 


Im 


\  /  1  \  6  j  g  ^^iPox 

_k  V  ^kPox)  \/kPox  [  V^lPox 

g2inx{ki)p„^  f  Re(^/)  Im(A:/)  I 

Re  (A:,) 

""  W^\  ’ 


♦1 


(5.87) 


and 


2 


g-2Im(fcj)Po* 

I  ^iPox  I 


(5.88) 


The  exact  same  result  is  obtained  for  the  two  TM  terms.  Plugging  these  results  into 
and  P™  yields 


pTE 


7r^Re(A;;)/>o 

128t(;/io 


{QilQi) 


(g|9i) 

(oTie,) 


F{Pox) 


X  (1  ^no)  [^n(m— 1)  <^n(m+l)  (1  "H  ^mo)]  ) 


(5.89) 


pTM 


7r^Re{ki)p^ 


( Q;  I  Re  (e,.)  I  Qi ) 


(g'IQi) 


128  {upo)  U}£o  (  <?/  I  Sr  I  Ql) 

X  [^n(TO+l)  (1  “f  ^mo)  “t"  ^n(m—l)  (1  d"  *^ml)] 


F(Pox) 


(5.90) 


From  (5.89)  and  (5.90),  the  role  of  the  p  and  z  functions  becomes  more  clear.  The  Hankel 
functions  effectively  cancel,  assuaging  my  earlier  concern  that  the  results  may  be  incorrectly 
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influenced  by  my  choice  of  “test”  function  .  The  only  remnant  of  the  radial  parasitic 

mode  functions  is  the  Re  (k)  term,  which  is  a  normal  weighting  term  in  the  Poynting  vector^^ 
Moreover,  this  supports  my  earlier  thesis  that  the  Im(fc;)  does  not  affect  the  radiated  power 
calculation. 

The  F{Pox)  f^rm  weights  the  power  by  the  lasing  mode  penetration  into  the  cladding 
region,  thereby  explicitly  incorporating  the  confinement  factor.  The  axial  inner  products  give 
an  intuitive  weighting  factor  due  to  the  parcisitic  mode  overlap  with  the  lasing  mode.  Note 
that  there  is  no  need  to  normalize  the  parasitic  modes,  since  any  constant  factor  will  cancel 
in  the  radiated  power  expressions.  Similarly,  the  same  comment  holds  for  the  lasing  mode,  as 
evidenced  by  (5.72).  This  is  only  logical:  any  normalization  just  scales  the  mode  amplitude. 
Finally,  the  constants  and  delta  functions  are  simply  dimensional  scaling  factors  necessary  to 
couple  the  various  LP,„  lasing  and  TE„/TM„  parasitic  modes.  In  the  next  chapter,  I  apply 
this  approach  to  calculate  threshold  gains  (Figure  6.3)  and  currents  in  actual  VCSELs;  I 
find  the  results  agree  qualitatively  with  experiments  [33]. 


^^This  can  be  understood  by  considering  a  plane  wave:  E 
Poynting  vector  S  ^ 


-ikz 


,  H 


-j^e  which  implies  the 
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VI.  Spontaneous  Emission  and  Laser  Analysis 

Ultimately,  an  optical  analysis  is  incomplete  until  it  is  coupled  to  a  semiconductor  gain 
model  to  calculate  measurable  lasing  parameters.  These  are:  light  verses  current,  threshold 
current,  and  the  percentage  of  the  total  spontaneous  emission  that  couples  into  the  lasing 
mode  [P).  In  this  chapter,  I  combine  my  weighted  index  LP  mode  calculation — including 
diffraction — with  J.P.  Loehr’s  gain  model  [42]  to  calculate  the  laser  operating  characteristics 
for  a  low  threshold  VCSEL  [33]  structure  for  various  oxide  aperture  radii. 

To  estimate  the  spontaneous  emission  properties,  I  compute  a  set  of  three-dimensional 
cavity  eigenmodes.  For  each  mode,  I  find  the  energy  (^),  spatial  field  profile 
confinement  factor  (r‘°*),  and  threshold  gain  (gth).  These  modes  are  then  applied  to  a 
new  spontaneous  emission  model,  tailored  to  incorporate  the  partial  transverse  confinement 
afforded  by  oxide  apertures.  In  past  approaches,  the  optical  modes  were  either  dealt  with 
as  a  two  dimensional  continuum  [10,18,44,75],  which  is  appropriate  for  VCSELs  with  no 
transverse  confinement,  or  as  a  discrete  set  of  modes  of  a  fully  enclosed  three-dimensional 
cavity  [4,5,73,77],  which  is  most  relevant  for  etched  posts.  Since  the  oxide  aperture  only 
results  in  partial  transverse  field  confinement,  I  model  the  mode  structure  as  the  median  of 
the  two  extremes.  I  use  a  discrete  mode  sum,  due  to  the  WIM  eigenmodes,  superimposed 
on  a  background  two-  and  three-dimensional  density  of  (plane  wave)  modes  to  represent 
the  complete  set  of  photon  states  which  interact  with  the  gain  media.  This  is  analogous  to 
modeling  the  continuum  states  of  a  quantum  well  as  a  three  dimensional  density  of  states 
plus  a  discrete  set  of  resonances. 

Steady-state  light  verses  current  curves  are  calculated  by  solving  the  photon  and  car¬ 
rier  multimode  related  rate  equations  using  the  LP  eigenmodes,  Loehr’s  gain  model,  and 
the  spontaneous  emission  results.  By  extrapolating  these  curves,  the  threshold  current  is 
estimated  and  compared  with  measurement. 
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6.1  Test  Structure 


I  perform  the  lasing  mode  analysis  on  the  University  of  Texas’s  semiconductor-dielectric 
composite  980  nm  VCSEL  [31,33].  The  A/2  cavity  is  formed  by  Alo.97Gao.03As  surrounding 
a  single  80  A  Ino.20Gao.80As  quantum  well  sandwiched  between  pairs  of  100  A  GaAs  and  100 
A  Alo.67Gao.33As  barrier  layers.  The  lower  DBR  is  a  26  period  n-type  GaAs/AlAs  DBR.  The 
top  DBR  is  a  post-growth  electron-beam  deposited  5  period  CaF-ZnSe  dielectric  DBR.  A  p"*’- 
GaAs  quarter  wave  cap  layer  is  included  between  the  cavity  and  the  upper  dielectric  DBR  for 
electrical  contacting.  The  entire  structure  and  the  standing  intensity  pattern  are  illustrated 
in  Figure  6.1.  In  my  model,  I  used  n  =  1.35  and  2.54  for  CaF  and  ZnSe,  respectively. 
The  index  and  free  carrier  loss  data  for  the  semiconductor  materials  are  listed  in  Table  3.1, 
Chapter  III. 

Circular  apertures  are  formed  in  this  VCSEL  by  oxidizing  the  top  Al0.97Ga0.03 As  cavity 
layer.  The  fact  the  oxides  lay  within  the  cavity,  a  unique  design  approach,  has  been  credited 
with  the  low  thresholds  obtained  from  this  device.  The  nearness  of  the  oxide  to  the  quantum 
well  suggests  the  aperture  will  have  a  strong  impact  on  the  emitted  fields,  making  this 
structure  an  ideal  test  subject  for  my  microcavity  VCSEL  analysis. 

6.2  Optical  Modes 

Using  the  LP  mode  analysis  of  the  previous  chapter,  I  calculated  the  eigenmode  data 
for  the  test  VCSEL  with  oxide  aperture  radii  ranging  from  0.5  ^m  to  4.0  fim  in  0.5  pm  steps. 
For  each  radius  I  found  the  pertinent  optical  mode  parameters  for  each  discrete  eigenmode. 
In  Figures  6.2  -  6.4, 1  plot  the  mode  energy,  radial  confinement  factor,  and  threshold  gain  for 
the  LPoi  and  LPu  eigenmodes  as  a  function  of  aperture  radius.  The  longitudinal  confinement 
factor  is  roughly  0.02  for  all  cases. 

As  I  have  now  come  to  expect,  the  mode  energy  blueshifts  and  the  transverse  confine¬ 
ment  decreases  with  decreasing  aperture  radius.  The  higher  order  modes  are  less  confined 
and  blueshift  more  than  the  lower  order  modes.  The  order  of  magnitude  and  behavior  of  the 
calculated  threshold  gain  with  aperture  radius  is  in  good  qualitative  agreement  with  exper¬ 
iment.  However,  to  determine  the  threshold  current — the  true  threshold  metric — I  need  to 
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Figure  6.1  Index  and  standing  intensity  profile  for  the  University  of  Texas’s  low  threshold 
A/2  cavity,  dielectric  DBR  VCSEL. 
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Figure  6.3  LPoi  and  LPn  mode  threshold  gain  vs  radius.  The  solid  and  dashed  curves  are 
a  cubic  spline  fit  of  the  discrete  data. 
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Figure  6.4  LPoi  and  LPn  mode  transverse  confinement  factor  vs  radius.  The  solid  and 
dashed  curves  are  a  cubic  spline  fit  of  the  discrete  data. 
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incorporate  the  eigenmodes  with  Loehr’s  quantum  well  gain  results  [42,43]  into  a  full  laser 
model.  I  do  this  in  the  next  two  sections. 


6.S  Spontaneous  Emission 

Before  I  solve  the  related  rate  equations  and  compute  the  light  verses  current,  I  must 
calculate  the  spontaneous  emission  rate  {Rsp}-  To  do  this,  I  combine  my  WIM  eigenmodes 
with  Loehr’s  bandstructure  and  gain  model.  Spontaneous  emission  rates  are  then  calculated 
using  a  new  hybrid  approach  that  combines  emission  into  the  photon  continuum  and  emission 
into  localized  cavity  eigenmodes.  In  both  cases,  the  spontaneous  emission  rate  is  calculated 
using  Fermi’s  golden  rule  [16,47]. 

To  calculate  optical  emission  and  absorption  rates,  only  transitions  between  the  con¬ 
duction  subbands  and  the  heavy-  and  light-hole  states  are  considered;  transitions  to  and  from 
split-off  bands  are  ignored.  The  quantum  well  material  gain  for  a  cavity  eigenmode  (m)  (not 
to  be  confused  with  the  field  azimuthal  mode  number)  with  energy  hwm  and  (approximate) 
plane  polarization  em  is  obtained  using  the  standard  single-particle  gain  expression  [42] 


g{€m,TkJm) 


VcmlLHu,^  (2ir)  ^ 
rf  (27r)  Jq 


(6.1) 


Here  rj  is  the  refractive  index  of  the  cavity,  c  is  the  speed  of  light,  mo  is  the  free  electron 
mass,  e  is  the  electronic  charge,  L  is  the  quantum  well  width,  fe(E,  p)  is  the  Fermi-Dirac 
occupation  function  for  electrons  with  energy  E  and  chemical  potential  /x,  and  ko  is  defined 
by  the  energy- conservation  condition 


(6.2) 
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The  optical  matrix  elements  P^^'(A;o,^6)  are  given  by 


nn 


'(kt)=  X)  (/nkJ/i°®^)KIPKB) 


(6.3) 


where 

(/s,i/r’> = / (6-4) 

The  functions  and  fn^  refer  to  the  envelope  functions  for  electron  and  holes,  respec- 
tively,  in  the  quantum  well;  the  corresponding  energies  are  E\  ’  and  En  •  The  conduction 
and  valence  {p^)  quasi-Fermi  levels  are  set  under  high-level  injection  by  forcing  charge 
neutrality  n  =  p  in  the  undoped  quantum  well:  occupation  of  the  continuum  states  above 
the  barriers  is  also  included. 

The  gain  results  are  broadened  by  convolution  with  a  5  meV  half-width  Lorentzian 
lineshape  function  to  obtain  realistic  gain  spectra.  Figure  6.5  shows  the  calculated  material 
gain  for  plane-polarized  (TE)  modes.  Since  it  is  very  difficult  to  calculate  the  exact  peak 
gain  wavelength,  the  peaks  of  the  gain  curves  in  Figure  6.5  are  aligned  with  the  energy  of 
the  lasing  eigenmode  by  a-posteriori  shifting  of  the  gain  spectrum. 

The  spontaneous  emission  rate  for  coupling  to  a  three-dimensional  photon  continuum 
of  plane  waves  is  given  by  [42] 


^Tje  nuj  1  Aco 

3m§»V  (2^)  ^  I  -  iE^'‘\ko) 

1  /*27r 

/t/t/  -»'=1  9  a  V  /  *'0 


(6.5) 


Similarly,  the  spontaneous  emission  rate  for  coupling  to  a  ^ wo- dimensional  photon  continuum 
for  emission  in  the  x  —  y  plane,  with  the  photon  wavefunction  localized  to  one  wavelength 
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Figure  6.5  Material  gain  for  an  80  A  Ino.2Gao,8As/Alo.2Gao.8As  quantum  well  for  various 
two-dimensional  carrier  densities  (n). 


along  the  z  axis,  may  be  calculated  as 


rje^huij 
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(6.6) 


Essentially,  this  represents  an  estimate  of  the  spontaneous  emission  into  all  parasitic  modes, 
whether  or  not  they  are  resonant  with  the  eigenmode.  Since  the  transitions  in  strained 
quantum  wells  occur  predominantly  between  heavy  hole  bands,  and  since  z-polarized  light 
does  not  couple  to  pure  heavy-hole  states,  the  two  rates  may  be  related  by  the  approximate 
ratio 


(6.7) 

Finally,  I  consider  emission  into  the  localized  cavity  eigenmodes  (m)  calculated  in 
Section  6.2.  Since  I  have  already  expended  great  effort  to  compute  the  electric  field  profiles 
of  these  modes,  it  is  most  useful  to  express  the  rate  in  terms  of  the  average  electric  field 
(jErnCO))  at  the  center  of  the  active  region.  In  the  dipole  approximation,  this  rate  can  be 
calculated  as 


Note  that  this  transition  strength  depends  strongly  on  the  “volume”  |(£^to(0))  1”^  of  the  cavity 
eigenmode  (m):  smaller  cavities  will  more  tightly  confine  the  mode  and  result  in  stronger 
transitions. 
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The  total  spontaneous  emission  rate  will  be  a  combination  of  3D,  2D  and  localized  mode 
emission.  I  assume  that  2D  photons  are  always  present,  since  the  vertical  layer  structure 
supports  many  (parasitic)  guided  modes  that  emit  radially.  Also,  I  assume  that  the  3D 
photon  continuum  begins  once  the  photon  energy  exceeds  that  of  the  highest-energy  cavity 
mode,  which  I  denote  by  Ezd-  Thus,  I  approximate  the  total  spontaneous  emission  rate  as 

/*-^3D  /*00 

^  0  JEzd 

Here  I  have  explicitly  indicated  that  depends  on  the  quantum  well  carrier  concentra¬ 

tion  n  through  the  electron  occupation  function  /e. 

6.4  Related  Rate  Equations 

The  optical  gain,  spontaneous  emission,  stimulated  emission,  injected  radiative  current 
(J,ad),  and  carrier  concentration  are  linked  by  the  multimode  related  rate  equations  [11,42], 


=  -  fi^m)  -  "m]  Sm  -  R^\  (6.11) 

V 

=  —  -  igg^'^^(n)  -  -  (6.12) 

Here  Sm  is  the  photon  density  inside  the  cavity  for  a  localized  mode  m.  Since  the  two-or 
three-dimensional  continuum  photons  immediately  leave  the  cavity  after  they  are  emitted, 
I  assume  that  none  of  these  can  build  up  significant  density.  Because  I  consider  very  small 
cavity  volumes,  at  most  a  few  dozen  localized  cavity  modes  (m)  are  considered  in  (6.11)  and 
(6.12). 

I  solve  these  rate  equations  only  in  steady-state.  In  this  case,  the  left-hand-sides  of 
(6.11)  and  (6.12)  are  both  zero,  and  I  immediately  solve  (6.11)  for  Sm,  obtaining 


dSm 

dt 

dn 

dt 


Sm  = 


rj  P^m 


(6.13) 
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By  substituting  this  expression — for  each  mode  m — into  the  steady-state  equation  (6.12),  I 

obtain  .. 

...(n)  =  e  {<-)(n)  .  .  (6.14) 


By  inverting  (6.14)  I  compute  n(Jrad)  and,  in  turn,  -S'„^(Jrad)  via  (6.13).  Finally,  I  convert 
the  modal  photon  densities  inside  the  cavity  to  the  total  optical  power  Pout  emitted  from 
the  cavity  by  computing 


/>„„.  =  (6.15) 

^  7] 

m 

where  r  is  the  cavity  radius.  (I  include  emission  out  of  both  the  top  and  bottom  mirrors.) 
To  relate  the  laser  current  I  to  the  current  density  Jradj  I  consider  the  best-case  scenario, 
ignoring  carrier  leakage,  nonradiative  recombination,  and  current  spreading,  and  write 

I  =  (6.16) 

By  combining  (6.15)  and  (6.16),  I  plot  the  light  (Pout)  verses  current  (I)  for  a  variety  of 
cavity  radii  in  Figure  6.6. 

My  calculated  light  verses  current  curves  give  good  agreement  with  measured  thresh¬ 
olds.  The  measured  threshold  current  in  a  2  fim  square  aperture  VCSEL  fabricated  from 
this  layer  structure  is  91  fiA  [33].  My  results  underestimate  this  number  by  roughly  a  factor 
of  three.  My  “best  case”  analysis  is  most  likely  to  blame  for  this  disparity;  I  may  also  have 
underestimated  the  free-carrier  absorption  losses.  Note  that  the  curvature  near  threshold 
becomes  more  pronounced  as  the  radius  decreases,  indicating  that  a  greater  fraction  of  spon¬ 
taneous  emission  is  coupled  into  the  lasing  mode.  Likewise,  the  external  quantum  efficiency 
decreases  sharply  for  the  smallest  devices,  indicating  that  much  of  the  optical  power  is  lost 
to  diffraction. 

By  evaluating  the  fraction  of  the  total  emission  that  enters  the  lasing  mode  at  “threshold” 
which  can  be  ambiguous  [9] — I  compute  the  spontaneous  emission  factor  ^  as 
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Figure  6.6  Output  power  (from  both  mirrors)  verses  current  for  various  device  radii. 
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(6.17) 


^  ~  ■ 

In  Figure  6.7  I  plot  as  a  function  of  cavity  radius.  Note  that  even  though  /?  increases 
significantly  as  the  cavity  radius  shrinks  from  4  to  1.5  fiva.,  it  still  does  not  approach 
unity  even  for  r  =  1.5  /zm.  The  reason  for  this  is  that  there  is  always  significant  coupling 
to  the  two-dimensional  “parasitic”  continuum  of  photons  radiating  out  through  the  DBR 
layers.  The  only  way  to  eliminate  this  coupling  would  be  to  fabricate  a  completely  separable 
structure  with  very  high  transverse  confinement,  such  as  the  design  discussed  in  Section  5.3. 
Then  the  coupling  to  the  two-dimensional  modes  would  be  eliminated,  and  |  operation 
could  become  possible.  (If  the  fundamental  mode  is  degenerate,  as  is  usually  the  case,  then 
I  is  an  upper  bound  for  p.) 
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Spontaneous  Emission  Factor  (P) 


Figure  6.7  Spontaneous  emission  factor  (/?)  verses  device  radius. 
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VIL  Conclusions  and  Recommendations 

7.1  Review  of  Present  Work 

The  first  vertical  cavity  surface  emitting  laser  was  reported  in  1979  [71].  Since  then, 
work  has  increased  exponentially  [17],  with  whole  conference  sessions  and  journal  special 
editions  dedicated  to  these  devices.  Much  of  the  appeal  of  these  devices  over  their  in-plane 
counterparts  results  from  their  radically  different,  and  in  several  respects  improved,  operating 
characteristics,  better  suiting  them  to  many  applications.  This  fact  has  only  recently  been 
recognized  by  the  commercial  community,  which  is  just  now  adopting  VCSELs  for  products 
such  as  print-heads  and  local  communication  fiber  sources. 

Despite  their  potential,  these  devices  have  not  yet  been  optimized,  and  the  ideal  device 
has  not  yet  been  designed,  let  alone  built.  In  fact,  the  VCSEL  fabrication  technology  is  just 
now  maturing  to  the  point  where  a  discussion  on  ideal  device  design  is  reasonable.  The 
latest  (and  perhaps  last?)  major  advancement  occurred  when  the  University  of  Texas  [33] 
researchers  integrated  an  oxide  aperture  layer  into  the  structure.  Electrically,  the  oxide  is 
insulating,  and  can  provide  excellent  current  confinement.  Optically,  the  oxide  index  is  lower 
than  its  host  crystal.  This  index  difference  may  provide  index  guiding,  or  perhaps  optical 
scattering  loss  to  the  lasing  mode.  This  dual  nature  of  the  oxide — providing  waveguiding 
and/or  scattering — has  created  substantial  confusion  with  respect  to  VCSEL  design.  The 
primary  goal  of  this  thesis  is  to  resolve  this  issue  by  understanding  the  physics  of  the  oxide 
aperture  within  the  VCSEL  cavity.  By  constructing  several  optical  models,  I  have  explained 
the  seemingly  opposite  reports  on  oxide  effects  and  optimal  oxide  design,  and  have  derived 
simple  rules  for  low  threshold  VCSEL  design. 

Mathematically,  the  optical  problem  is  a  complex,  second  order,  linear,  boundary  value 
problem.  The  principal  difficulty  arises  from  the  fact  that  the  governing  differential  equa¬ 
tions  are  not  separable  in  any  coordinate  system,  greatly  frustrating  analytic  solutions. 
Simple  Fabry- Perot  models,  the  staple  of  most  semiconductor  laser  analysis,  break  down  in 
the  interesting  realm  of  small,  highly  confining  oxide  apertures.  For  my  first  approach  at 
solving  this  problem  accurately,  but  economically,  I  employed  a  variational  technique:  the 
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weighted  index  method.  In  this  technique  the  best  separable  solution  is  found  in  a  vari¬ 
ational  sense.  This  iterative  method  converges  on  an  “eigenmode”  by  alternately  solving 
two  coupled  one-dimensional  integro-differential  equations.  The  seed  for  the  search  is  the 
standard  Fabry- Perot  solution.  Limiting  my  analysis  to  azimuthally  symmetric  structures, 
I  used  this  technique  to  solve  the  full  vector  Helmholtz  equation  to  calculate  the  basic  opti¬ 
cal  lasing  parameters:  wavelength,  confinement  factor,  and  spatial  field  distribution  for  the 
fundamental  and  higher  order  lasing  modes. 

From  the  full  vector  weighted  index  analysis  I  was  able  to  uncover  two  new  effects:  a 
polarization  dependent  change  in  mirror  reflectance — previously  only  noted  in  edge-emitting 
lasers — and  a  change  in  spectral  separation  of  the  lasing  modes  with  oxide  aperture  size 
and  radius.  The  spectral  mode  separation  results  from  the  dispersion  in  blueshift  rate  with 
aperture  size  between  the  various  cavity  modes.  I  confirmed  the  blueshift  by  comparison  with 
experiment  for  two  different  devices.  My  results  matched  the  measurements  better  than  any 
previous  calculations.  The  spectral  separation  led  me  to  a  new  concept  for  achieving  low 
threshold,  single  mode  lasing.  In  my  approach,  a  spectrally  sharp  gain  region  is  coupled  to 
the  modes  of  a  highly  confined  VCSEL.  For  a  sufficiently  small  and  highly  confined  VCSEL, 
all  but  the  lowest  order  lasing  mode  will  be  blueshifted  away  from  the  gain  region,  and 
only  a  single  mode  will  lase.  Ideally,  the  device  would  be  made  as  small  as  possible  for 
maximum  mode  separation.  However,  such  small  devices  suffer  enormous  diffraction  losses, 
unless  a  truly  separable  device  geometry  could  be  fabricated.  In  a  separable  structure,  the 
cavity  eigenmodes  would  theoretically  suffer  no  diffraction  and,  ignoring  absorption  losses, 
the  ultimate  quantum  optical  threshold  limit  could  be  achieved.  In  this  case,  the  polarization 
dependent  change  in  mirror  reflectivity  predicted  by  the  weighted  index  method  may  actually 
be  observed. 

Since  the  weighted  index  method  is  incapable  of  addressing  diffraction  on  its  own,  I 
also  employed  a  quasi- exact  computational  electromagnetic  technique — the  finite  element 
method.  My  implementation  of  the  finite  element  method  employed  new  body-of-revolution, 
second-order  vector  edge  and  node  elements.  With  the  finite  element  method  I  was  able 
to  capture  all  of  the  VCSEL  optical  parameters,  including  the  elusive  diffractive  loss.  By 
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applying  my  model  to  a  novel  high  index  confinement  VCSEL,  I  was  able  to  discover  the 
fundamental  mechanisms  governing  diffraction  in  apertured  VCSELs.  These  are: 

1.  The  density  of  parasitic  modes  in  the  spectral  vicinity  of  the  eigenmode. 

2.  The  relative  alignment  between  the  eigenmode  and  parasitic  mode  propagation  vectors. 

3.  The  transverse  confinement  of  the  optical  mode  to  the  core  region. 

These  three  factors  may  be  used  to  qualitatively  explain  both  the  University  of  Cal¬ 
ifornia  at  Santa  Barbara  and  University  of  Texas  group’s  oxide  design  results.  The  Santa 
Barbara  group  has  suggested  that  thinner  oxides  placed  at  standing  wave  nodes  are  supe¬ 
rior  [29,30].  They  consider  semiconductor  DBR  VCSELs  with  longer  effective  cavity  lengths 
and  relatively  laxge  parasitic  mode  densities,  approaching  the  three-dimensional  density  of 
modes.  As  a  result,  even  weak  coupling  to  the  parasitic  mode  continuum  will  yield  a  large 
diffractive  loss,  implying  that  thinner  oxides  are  better.  The  Texas  group  works  with  a 
VCSEL  with  one  semiconductor  and  one  dielectric  DBR  [20].  In  their  case,  the  cavity  and 
corresponding  optical  density  of  modes  is  much  smaller,  and  the  confinement  aspect  of  their 
oxide  aperture  is  the  dominant  effect.  Hence,  they  prefer  thick  oxides  at  standing  wave  peaks 
near  the  gain  region. 

The  success  of  the  parasitic  mode  picture  in  describing  diffraction  loss  prompted  me 
to  try  and  incorporate  these  modes  within  the  weighted  index  formalism.  I  successfully 
accomplished  this  by  coupling  a  set  of  resonant  parasitic  modes  to  the  eigenmode  using  a 
rigorous  electromagnetic  treatment  of  the  coupling.  In  this  treatment  the  implicit  source  of 
the  parasitic  modes  is  the  weighted  index  method  mode  error.  In  the  process  of  incorporating 
the  parasitic  modes,  I  simplified  my  weighted  index  model  to  a  scalar  solution.  This  resulted 
in  a  huge  savings  in  complexity  and  computational  effort  with  little  cost  in  accuracy.  The 
only  real  cost  was  that  the  scalar  model  forfeits  the  change  in  mirror  reflectivity.  But  since 
this  has  not  yet  been  observed,  this  was  not  much  of  a  sacrifice.  The  advantages  of  this 
new  approach  are  immense.  The  weighted  index/parasitic  mode  calculation  should  provide 
the  best  semianalytic  optical  solution  for  apertured  VCSELs,  incorporating  all  the  pertinent 
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physics  within  the  field  solution.  Furthermore,  it  is  fast  and  does  not  require  a  supercomputer 
to  run,  a  distinct  advantage  over  my  finite  element  calculation. 

Of  course,  you  can  model  optical  physics  until  you  are  blue  in  the  face,  but  you  won’t 
raise  an  experimentalist’s  eyebrow  until  you  can  compute  measurable  device  parameters.  For 
this,  I  combined  my  optical  data  with  Loehr’s  gain  model  and  a  new  emission  model  con¬ 
taining  both  a  sum  of  discrete  modes  and  an  integral  over  both  a  two-  and  three-dimensional 
continuum.  This  approach  is  designed  to  capture  the  incomplete  confinement  effects  of  oxide 
apertured  devices.  The  full  laser  model  was  applied  to  the  University  of  Texas’s  low  thresh¬ 
old  VCSEL  [33]  to  compute  all  the  standard  laser  parameters:  light  verses  current,  lasing 
threshold,  and  beta  factors.  My  results  are  in  good  qualitative  agreement  with  their  mea¬ 
surements,  differing  by  no  more  than  a  factor  of  three.  They  provide  a  good  initial  validation 
of  the  new  emission  model  and  my  new  weighted  index/parasitic  mode  computation. 

7.2  Suggestions  for  Future  Study 

As  is  often  the  case  in  research,  for  every  door  you  shut  two  more  open.  Indeed,  staying 
focused  on  a  specific  task  may  be  the  most  challenging  aspect  of  research.  I  certainly  found 
this  to  be  the  case,  with  many  research  avenues  presenting  themselves  to  me  through  the 
course  of  my  work.  Having  completed  my  dissertation,  I  feel  that  I  now  have  the  perspective 
to  see  where  my  work  fits  into  the  general  body  of  knowledge,  and  to  recommend  future 
studies. 

A  model  should  pass  through  three  distinct  phases:  derivation,  validation,  and  ap¬ 
plication;  I  have  completed  the  first  phase  and  have  just  begun  the  second.  Therefore,  I 
suggest  the  following  course  of  action.  First,  my  two  models  should  be  comprehensively 
tested  against  one  another.  Before  this  can  be  accomplished,  however,  the  finite  element 
linear  algebra  solver  must  be  revamped  to  address  larger  problems.  Second,  both  models 
should  be  tested  against  experiment  for  threshold  predictions.  Along  the  same  lines,  the 
complete  VCSEL  model  must  be  tested  and  refined  to  improve  its  overall  accuracy.  Third, 
the  models  should  all  be  applied  to  optimize  the  VCSEL  threshold.  This  is  pretty  much  the 
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whole  point  of  the  work:  the  ultimate  goal  is  to  design  a  better  device  and  take  another  step 
towards  the  “thresholdless  grail.” 
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Appendix  A.  Variational  Justification  of  the  WIM 

The  weighted  index  method  (WIM)  allows  me  to  compute  a  separable,  approximate  solution 
to  a  two  dimensional  partial  differential  equation.  In  this  appendix,  I  use  the  calculus  of 
variations  to  demonstrate  that  the  WIM  solution  represents  the  best  separable  approximation 
to  the  two-dimensional,  azimuthally-synometric  Helmholtz  equation. 

I  wish  to  solve 

where  ^  denotes  the  Laplacian  with  respect  to  the  radial  and  axial  coordinates  and  er{p,  z) 
is  the  relative  permittivity.  Equation  (A.l)  is  a  Sturm- Liouville  differential  equation  [3] 
of  the  form  Cij>  =  0  with  zero-Dirichlet  boundary  conditions  at  /)  =  oo,  zero-Dirichlet 
(m  ^  0)  or  Neumann  (m  =  0)  boundary  conditions  at  />  =  0,  and  zero-Robin  (mixed 
Dirichlet/ Neumann)  boundary  conditions  at  z  =  ±oo  (radiation  conditions)  [7]. 

I  solve  (A.l)  variationally  by  forming  a  functional  J[V’]  that,  when  extremized  over  the 
space  of  admissible  functions  V’(P5  ^r),  minimizes  the  Hilbert  space  distance 

/CX>  AOO 

/  li’ip,  z)  -  V’exact(p,  I V  dpdz.  (A.2) 

•oo  Jo 

Here  V’exact  denotes  the  exact  solution  to  (A.l)  [64].  I  refer  to  the  particular  function  mini¬ 
mizing  (A.2)  as  the  best  variational  solution.  The  variation  6J[ip]  is  given  by 

=  {vL+ = 

(  Si;  I  \  i;)  +  ‘^{6i;  \  £r{p,  z)\  i;)  -  {Si;  \  p~'^  \  i;)  =  0,  (A.3) 

c 

where  the  appropriately-weighted  inner-product  for  this  geometry  is 

/OO  1*00 

/  {A*B)pdpdz.  (A.4) 
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Integrating  the  first  term  on  the  second  line  of  (A. 3)  by  parts  yields 

=  ^  ( I £r(p, z)\i’) -rn^  I  (A.5) 

where  the  surface  terms  drop  out  for  bound  modes.  Then  J[^]  is  found  by  bringing  the 
variational  operator  outside  the  inner-products  to  yield 


m  =  \\%W  z)  I  V-  >  -  <  V’  I  I  V„v- )]  (A.6) 

Zt  C- 


I  now  search  for  the  optimal  function  0  that  extremizes  J,  i.e.  that  gives  ^J[^]  =  0.  This 
procedure  yields  differential  equations  of  constraint — the  Euler  equations — on  V*-  If  I  allowed 
i})  to  range  over  the  complete  space  of  complex  two  dimensional  functions,  the  constraint 
equation  would  be  (A.l)  and  the  whole  exercise  would  be  useless.  However,  by  restricting 
the  functional  domain  of  V*  I  can  generate  different  constraint  equations.  In  particular,  I  can 
find  the  best  separable  solution  by  taking  the  variation  of  (A.6)  while  restricting  V’  to  the 
space  of  all  separable  functions 

i,(p,z)=F{p)G{z).  (A.7) 

To  generate  constraint  equations  for  F  and  G,  I  take  the  total  variation  AJ[i^,  G\.  This 
is  given  by  [23] 

AJ[F,  G]  =  S[F  +  f,G  +  g]- J[F,  G],  (A.8) 

where  /  and  g  are  incremental  deviations  from  F  and  G.  The  first  variation  (6J)  is  the 
principal  part  of  the  total  variation  given  above.  Substituting  F-f/  and  G-1-^  appropriately 
into  (A.6)  and  then  taking  the  principal  part  (the  part  involving  single  terms  of  /  and  g  to 
the  first  power)  yields 

alf,G]  =  !^^^e,(p,z)-y^{Ffa^  +  F^Gg)- 

(F‘f'CP  +  FVsf  +  F/(G')^  +  }  pipi’:, 


(A.9) 


where  the  prime  denotes  differentiation  with  respect  to  the  argument.  Integrating  by  parts 
and  applying  the  appropriate  Dirichlet  and  Neumann  boundary  conditions  on  F{p)  and 
G{z),  I  remove  the  derivatives  from  /  and  g  to  yield 


6S[F,G]  = 

f[/' 


oor  poo  r  r^2 


e,{p,z)-y  FG^-F{G'y  +  ^{pF')'G^^dz  fpdp 


(A.IO) 


r  \  r  f 

J-oo  [Jo  I  L  ^ 


er{p,  z)-^  F^G  -  {F'fG  +  F^G"  \pdp\  gdz. 


Since  either  f{p)  or  g{z)  may  equal  zero,  this  implies  that  both  the  first  and  second  terms 
must  independently  equal  zero  to  extremize  J[^].  Thus  I  generate  from  (A. 11)  the  two 
independent  variational  equations: 

F  r  {\‘4Mp,z)-'^]a\z)-{G'{z)f]dz  +  -(pF')'  r  G\z)dz  =  0,  (A.ll) 
J-oo  Pi  )  P  J-oo 


aj"  I  cn"}  pip  +  [F‘‘pdp  =  0.  (A.12) 


Employing  the  axial  and  radial  inner-products 


/OO 

■OO 


(A.13) 


poo 

{F\F)=  IFfpdp, 
Jo 


I  can  rewrite  (A.ll)  and  (A.12)  as 


t(/’-P')+( - (gTg) - ^1^-® 


(A.14) 


(A.15) 


(A.16) 
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Introducing  further  notation  and  applying  the  product  rule  to  (A.  15),  I  obtain  the  weighted 
index  equations 

no)  +  ^no)  +  ol  -  y)  Hp)  =  0  ( A.17) 

and 

G"(z)+f,,j[F,z]G{z)  =  0,  (A.18) 


where, 


,,  ,r  ,{‘pVp'‘){G\e,\G)  +  (G\a") 

KiflG^  P]  = - (OjG) - 


(A.19) 


and 


iu;yc‘^){F\er\F)  +  {F\F>')  +  {F\p-^\r)-m^{F\p-^\F) 

UW)  ' 


(A.20) 


These  variational  constraint  equations  are  exactly  the  same  as  (2.14)  and  (2.15)  from  Section 
2.2.  Therefore,  the  solution  of  these  equations  represents  the  best  separable  approximation 
to  (A.l). 
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Appendix  B.  Transfer  Matrix  Generalization  of  Fabry-Perot  Laser  Theory 

In  this  appendix,  I  show  that  the  lasing  conditions  generated  via  transfer  matrix  techniques 
are  in  fact  a  generalization  of  the  classic  Fabry-Perot  lasing  conditions.  Specifically,  I  show 
that  for  the  simple  plane- wave  case  the  transfer  matrix  lasing  conditions  reduce  to  the  Fabry- 
Perot  conditions. 

In  classical  Fabry-Perot  laser  theory  the  threshold  gain  and  lasing  mode  wavelengths  are 
given  by  a  condition  on  the  round-trip  field  amplitude  and  phase.  In  order  for  a  lasing  mode 
to  be  supported  the  round-trip  change  in  amplitude  must  be  >  1.  This  is  the  gain/amplitude 
condition 

=  1,  (B.l) 


or 

Here  is  the  total  confinement  factor,  L  is  the  cavity  length,  is  the  mirror  loss,  Ri 

and  i?2  are  the  mirror  power  reflectivities,  and  I  have  ignored  any  losses  within  the  cavity.  In 
order  for  a  mode  to  be  supported  the  round-trip  change  in  phase  must  be  an  integer  multiple 
of  27r  or 

or 

2kL  =  m27r  A  =  — ,  (B-3) 

m 

where  m  is  an  integer.  These  conditions  are  illustrated  in  Figure  B.l. 


In  transfer  matrix  theory  the  fields  on  either  side  of  the  VCSEL  are  related  to  one 
another  by  a  2  x  2  transfer  matrix.  In  Figure  B.2,  the  traveling  wave  amplitudes,  Ei+  and 
El-,  in  the  substrate  region  at  the  substrate-to- VCSEL  boundary  are  related  to  the  traveling 
wave  amplitudes,  £'3+  and  E3-,  in  free-space  at  the  VCSEL-to-air  boundary  through  a 
transfer  matrix. 
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(B.4) 


A  lasing  condition  is  generated  by  requiring  no  incident  radiation  on  either  end  of  the 
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Figure  B.l  Illustration  of  classical  laser  amplitude  and  phase  conditions  generated  via 
Fabry-Perot  theory. 
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VCSEL,  i.e.  E3-  =  Ei+  =  0.  Hence,  for  a  non-trivial  solution  to  (B.4)  I  require  m22  =  0, 
where  m22  is  complex,  implying  Re(m22)  =  Im(m22)  =  0  (two  conditions). 

Based  on  the  effective  material  parameters  shown  in  Figure  B.2, 

m22  =  (^32r2i  +  ,  (B.5) 

132*21 

where  <32,  <21  and  r32,  r2i  are  the  effective  Fresnel  transmission  and  reflection  coefficients  and 
(^2  =  is  the  effective  phase  change  for  a  single  pass  through  the  VCSEL  of  effective 

length  L2  with  effective  propagation  constant  '  realistic  VCSELs,  the  first  term  in 
(B.5),  can  not  equal  zero  or  infinity.  Hence,  the  lasing  condition  reduces  to 

g»2^2  =  -r32r2i  =  r23r’2i-  (B.6) 


The  Fabry-Perot  amplitude  condition  can  be  recovered  by  taking  the  absolute  value 
squared  of  (B.6), 

=  |r23*’2ip  =  R23R21,  (B.7) 

where  Kth  is  the  effective  threshold  material  gain,  and  R21  and  R23  are  the  effective  mirror 
power  reflectance  coefficients  for  a  field  impinging  from  inside  the  VCSEL.  Recognizing  that 
=  2jKth,  (B.7)  may  be  rewritten  as 


ptot 


R21R23 


(B.8) 


which  is  equivalent  to  (B.2)  with  the  substitution  of  the  appropriate  effective  parameters. 

The  Fabry-Perot  phase  condition  is  found  by  equating  the  phase  terms  on  either  side 
of  (B.6).  Expanding  each  side  of  the  equation  yields 


gi2^n2L2  =  +  iKth  -  n3\  / +  ^ 

\n2  +  iKth  +  nsj  V?*2  +  iKth  +  ^i/  ’ 


(B.9) 


where  °Z  denotes  the  complex  phase  angle.  For  realistic  VCSEL  structures  Kth  ~  0-01  while 
ni  —  n2  >  0.5  and  n2  —  >  1.0,  hence  I  may  assume  the  phase  of  the  right  hand  side  of 
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(B.9)  is  approximately  zero.  Thus, 

2— ^2^2  =  2m7r  X  =  (B.IO) 

c  Tn 

which  is  equivalent  to  (B.3)  with  the  substitution  of  the  effective  cavity  length  L2  for  L. 
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Figure  B.2  Illustration  of  transfer  matrix  solution  to  the  plane-wave  VCSEL  problem. 


Appendix  C.  Vector  Basis  Functions  for  m=l  Modes 

In  this  appendix,  I  review  the  vector  basis  functions  employed  in  my  finite  element  model. 
Based  on  the  literature,  I  choose  to  span  the  mesh  using  second  order  node,  edge,  and 
face  element  functions  [25,59].  These  functions  have  been  shown  to  be  a  good  compromise 
between  mesh  density  and  function  order,  roughly  minimizing  the  total  number  of  unknowns 
required  to  obtain  a  given  solution  accuracy.  For  each  triangle  in  the  mesh  there  are  14 
basis  functions:  6  node  based  (iVo  -  iVg) ,  6  edge  based  (Wq  -  ,  and  two  face  based 

(#6,  Wt)  ,  as  illustrated  in  Figure  C.l  [35,  59].  To  define  the  14  element  functions,  I 
use  simplex  (or  barycentric)  coordinates,  defined  over  the  triangular  element  via  the  affine 
transformation  [13,68], 

[  1 1  r  1  1  1 1  r  Co 

P  =  Po  Pi  P2  Cl  • 

Z  Zo  Zi  Z2  C2 

Here  {pq,  zq),  (Pi,zi),  and  (pj,  ^2)  are  the  three  corners  of  the  triangle,  and  Co5  Ci^and  C2  are 
the  simplex  coordinates.  Based  on  the  Cj,  fbe  element  vector  functions  are 


No  =  (2Co-1)Co(H^),  (C.2) 

Ni  =  (2C,-l)Ci(^  +  p),  (C.3) 

N2  =  (2C2-1)C2(^  +  p),  (C.4) 

F3  =  4CoCi  +  p)  ?  (^-5) 

-^4  =  4CiC2(^  +  p)) 

=  4CoC2(^  +  p)>  (^•’^) 


C-1 


Wo  =  pCi^rC2  =  ^{b2P  +  C2z}, 

#1  =  pC2^r(l  =  ^{blP  +  C^z}, 

W2  =  K2V.Co  =  ^{&oP  +  coz}, 

^3  =  pCo^rC2  =  ^{b2P  +  C2z}, 

W4  =  p(o^rCl  =  ^{blP  +  Clz}, 

Ws  =  pCi^rCo=^{boP  +  coz}, 

We  =  pHl{C2^rCo-Co^rC2) 

=  ^^{C2{boP  +  coz)-Co{b2p  +  C2z)}, 

W7  =  p  4^2  (Co'^tCi  ~  Ci^tCo) 

=  ^^{^o{biP  +  ciz)-C^{hop  +  coz)]. 


Here 


d  ^ 

=  5?'’  + 


the  hi  and  c,-  are  given  by  the  inverse  affine  transformation, 


1 

0 

1 _ 

0 

1 

Cl 

1 

”  A 

a\  hi  Cl 

p 

1 

(N 

z 

and  A  is  twice  the  area  of  the  triangle, 


2  2  2 

A  =  Pi  [Zii  -  Zi2)  =  ^  pihi  =  ^  ZiCi 
t=0  «=0  i=0 


(C.8) 

(C.9) 

(C.IO) 

(C.ll) 

(C.12) 

(C.13) 

(C.14) 

(C.15) 


(C.16) 


(C.17) 


(C.18) 


(i,  il,and  i2  are  modulo  3).  These  functions  are  illustrated  in  Figures  C.2  and  C.3.  These 
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functions  have  three  important  properties  that  make  them  ideal  for  electromagnetic  calcu¬ 
lation  in  cylindrical  coordinates. 

First,  when  filling  the  S  and  T  matrices,  functions  based  on  the  same  node  or  edge  of 
neighboring  elements  (not  necessarily  with  the  same  element  function  number  No  —  W7)  are 
assigned  the  same  coefficient  Xi.  As  a  result,  the  node  and  edge  functions  naturally  enforce 
tangential  field  continuity  (as  illustrated  in  Figure  C.4)  without  necessarily  forcing  normal 
field  continuity,  due  to  the  two  face  functions. 

Second,  my  element  functions  are  modified  from  the  more  standard  form  [35]  to  force 
field  regularity  at  p  =  0  (the  axis  condition).  The  lasing  modes  I  seek  are  analogous  to  the 
(m  =  1)  LPoi  fiber  modes  [58],  where  the  azimuthal  field  dependence  is  sin(^)/  cos{<f>).  The 
proper  axis  condition  for  the  m  =  1  modes  is  [74] 

(C.19) 

and 

lim£^2=0.  (C.20) 

p-fO 

This  condition  is  included  in  (C.2)  -  (C.15)  by  weighting  the  “standard”  [35]  edge  and  face 
functions  by  p,  and  including  the  node  functions  as  part  of  the  p  field  expansion.  Note  that 
I  may  find  m  =  0  and  m  >  1  modes  using  a  basis  similar  to  (C.2)  -  (C.15),  altered  to 
accommodate  the  proper  axis  condition  for  these  modes. 

Third  and  least  obvious,  these  elements  properly  model  the  null  space  of  the  curl 
operator,  which  has  been  shown  to  eliminate  the  occurrence  of  spurious  solutions  (modes) 
[40,72]. 


lim  Ep  = 

p-i-O 
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Figure  C.2 


Illustration  of  the  six  node  based  functions  {No  —  Ns)  for  a  typical  element; 
the  functions  are  overlayed  on  a  triangle  outline  of  the  element. 
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Figure  C.3  Illustration  of  the  six  edge  based  {Wq  —  W^)  and  the  two  face  based  \  Wq  —  Wr) 
functions  for  a  typical  element;  the  functions  are  overlayed  on  a  triangle  outline 
of  the  element.  These  illustrations  do  not  include  the  p  weighting  present  in 
(C.8)  -  (C.15). 
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Element  B,  Wq 


Figure  C.4  Illustration  of  the  natural  tangential  continuity  between  mesh  elements.  By 
assigning  the  same  coefficient  Xi  to  Wi  of  element  A  and  Wq  of  element  B,  the 
vector  sum  of  the  two  functions  is  tangential  to  their  common  edge. 
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Appendix  D.  Absorbing  Layer  Design 

In  this  appendix,  I  give  the  details  of  a  new  absorbing  layer  design.  I  use  a  two  step  process 
to  optimize  b{p),  the  absorbing  layer  loss  function:  First,  I  obtain  a  rough  estimate  using 
an  asymptotic  description  of  the  fields  and  reflections.  Second,  I  fine-tune  the  layer  by 
minimizing  the  exact  reflection  values  as  obtained  by  a  rigorous  transfer  matrix  calculation. 
In  both  steps  I  model  the  radial  AL  (1)5)  as  a  set  of  discrete  cylinders,  as  shown  in  Figure 
D.l,  approximating  a(p)  as  =  a((/9„  -|-  p^-i)  /2),  where  AL  cylinder  n  is  defined  by 
[Pn-i )  Pn]  •  To  simplify  further,  I  ignore  the  2:  dependence  of  the  material  parameters  and 
perform  my  analysis  using  a  single  set  of  (sr,  the  rough  mean  material  values  for  the 
VCSEL.  This  should  be  sufficient,  since  I  am  not  concerned  about  interfacial  reflections 
between  the  various  axial  layers,  and  desire  only  to  suppress  nonphysical  radial  reflections.  I 
choose  b{p)  to  minimize  the  reflection  of  cylindrical  waves  incident  on  the  radial  boundaries. 

These  waves  are  constructed  from  the  z  and  ^  field  components  given  in  each  region 
(cylinder)  by 


Ez,n(^Pi  ■2) 

Z) 


E,j),niPj  z) 


(*»/>)  +  e<^-% 

—iuJSoSran  dEg 

kl  dp 

zcxjpfQa^  dHz 

K  dp 


(D.l) 

(D.2) 

(D.3) 


(D.4) 


where  Co  =  1  (60  =  0)  in  Cly  and  I  have  assumed  nonmagnetic  materials  {pj.  =  1).  The 
coefficients  and  are  the  magnitude  for  the  outward  and  inward  propagating  waves, 
and  pol  =  TE  or  TM  labels  the  (uncoupled)  polarizations  of  the  m  =  0  modes.  It  has  been 
shown  [45]  that  absorbing  layers  that  perform  well  for  these  m  =  0  modes  will  also  work  well 
for  the  m  =  1  modes  that  I  am  interested  in. 
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I  can  describe  cylindrical  waves  incident  at  an  angle  0 — measured  from  the  normal  to 
Fs  in  the  p  —  z  plane  (ups) — the  radial  interface  by  writing  the  radial  (A:„)  and  axial  (/S„) 
propagation  constants  as 

kn  —  cos  (^)  , 

=  korjan  sin  {9) . 

This  choice  ensures  that  and  satisfy  the  dispersion  relation 

in  each  region,  where  the  (generally  complex)  index  of  refraction  is 

r]  =  (D.8) 


D.l  Initial  Design:  Theory  of  Small  Reflections 

For  the  initial  AL  design  analysis  I  employ  the  theory  of  small  reflections  [60],  approx¬ 
imating  the  total  Fresnel  reflection  for  a  cylindrical  wave  incident  on  the  AL  at  pQ  as 

1  ^o|  +  I  ^1 1  +  I  1  ^1,2  +  .  .  ■  +  I  r^cyi  I  ■  (1^-9) 

Here  Ncyi  is  the  total  number  of  cylinders  in  the  AL,  r„  is  the  Fresnel  reflection  coefficient  at 
/9„,  and  gives  the  attenuation  due  to  the  propagation  from  pg  to  and  back.  I  use  |  r^l 

rather  than  in  (D.9)  to  minimize  the  interference  effects,  since  I  am  after  a  broadband 
AL  design  optimized  for  all  angles  of  incidence.  I  estimate  r„  using  the  asymptotic  form 
for  cylindrical  waves  in  each  region  n.  Due  to  the  form  of  (4.16)  and  (D.l)  -  (D.4),  the 
magnitude  of  the  TE  and  TM  Fresnel  reflection  coefficients  are  the  same,  hence  I  need  only 


(D.7) 


(D.5) 

(D.6) 
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perform  the  analysis  once.  Taking  the  asymptotic  form  for  the  Hankel  functions, 


lim  {Kp) 

p—*-oo 


,±iknP 


(D.IO) 


(+  applies  to  (1)  and  -  applies  to  (2))  and  enforcing  tangential  field  continuity,  I  find 


dn  ^n+1 

COS  (^„+i)  -  cos  ( 

«»)| 

cos (^n+l)  +  cos(0n) 

(D.ll) 


for  n  =  0  to  Ncyi  -  1.  and  ^n+i  are  the  (complex)  propagation  angles  with  respect  to 
riTs  in  regions  n  and  n  +  1,  respectively.  I  set  |  =  1  to  enforce  the  perfect  conductor 

boundary  condition  at  the  mesh  edge  (F).  Keeping  only  the  first  (lowest  order)  term  in  the 
derivatives  of  (D.IO),  I  approximate 


^1,2 . n  —  ^1^2  •  •■ifi 


(D,12) 


where 


—  ^ — i2knin 


tn  —  ^ 


(D.13) 


and  kn  and  In  =  Pn  —  Pn-i  radial  propagation  constant  and  cylinder  thickness  in 

region  n.  Substituting  (D.ll)  -  (D.13)  into  (D.9)  and  taking  the  magnitude  yields, 


r|  « 


Ncyl  n 

^  kni  n 


n=0  n=0 


(D.14) 


which  I  may  solve  for  a  given  a{p)  to  find  |r  (0j)l,  the  magnitude  of  the  total  Fresnel  reflection 
for  a  general  (e.g.,  includes  axial  propagation  exp{i/3z))  cylindrical  wave  incident  on  the  AL 
at  angle  6i.  Note  that  (D.14)  is  equivalent  to  the  result  obtained  via  a  plane  wave,  planar 
interface  analysis,  and  may  therefore  be  directly  applied  to  the  axial  AL  design  (fly  and 
Db),  in  addition  to  the  radial  design  (fls)- 
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D.2  Final  Design:  Transfer  Matrix  Solution 


To  finalize  my  radial  AL  design  I  use  a  2  x  2  transfer  matrix  solution  for  the  TE  or  TM 
fields.  Enforcing  continuity  of  (D.l)  and  (D.4),  or  (D.2)  and  (D.3)  at  each  radial  interface 
I  have 


L 


n 


(D.16) 


where 


(D.16) 


Rn+l  = 


(D.17) 


and 

I  have  dropped  the  pol=TE/TM  superscript  since  (D.15)  holds  for  both  polarizations.  At 
the  mesh  boundary  (T),  I  enforce  E<j,  =  0  which  relates  the  coefficients  in  the  outermost  AL 
cylinder  by 


^Ncyi  =  -Bn, 


Hi'\kN,„PNj 


(D.18) 


For  convenience,  I  choose  An^,  =  and  Bat..,  =  Beginning 

in  the  outer  AL  cylinder  and  working  inward  via  repeatedly  applying  (D.15),  I  calculate  Aq 
and  Bq.  The  total  Fresnel  reflection  coefficient  for  an  outward  propagating  wave  incident  on 
the  AL  is  given  by 


(D.19) 


D.S  Results  for  Structures  Analyzed 

To  design  my  AL,  I  ran  the  analysis  outlined  in  Sections  D.l  and  D.2  with  VCSEL 
material  parameters  er  =  9  and  /j.^  =  1.  I  assumed  a  polynomial  form  for  the  AL  loss 
function, 

order 

Kx)  =  E 

»=:1 
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Here  order  refers  to  the  polynomial  order,  bi  are  the  polynomial  coefficients,  and 


_  P  ~  Pts 
Pr-  Pts 


(D.21) 


is  the  normalized  radial  coordinate  ranging  from  0  at  F^,  the  VCSEL-AL  interface,  to  1  at 
r,  the  problem/mesh  boundary.  Through  a  numerical  optimization,  I  found  that  a  second 
order  polynomial  with  coefficients. 


h  =  0.1178,  (D.22) 

h  =  0.7433,  (D.23) 

worked  best  for  6(x)*  addition,  I  found  an  AL  thickness  of  1.5A  (at  A  =  870  nm)  was 
superior,  with  little  change  in  the  AL  properties  for  greater  thicknesses.  The  reflection 
coefficient  (calculated  with  the  transfer  matrix  approach)  verses  angle  is  given  in  Figure  D.2. 
The  key  point  to  recognize  from  Figure  D.2  is  the  intensity  reflection,  RdB  =  201ogio(|  r\)  < 
—50dB  for  6i  <  40°,  indicating  a  good  “wide-angle”  absorber  design. 

To  provide  a  more  practical  test  of  my  AL  design,  I  analyzed  the  lA-lTHIN  VCSEL 
(introduced  in  Section  4.4),  tracking  the  lasing  mode  (defined  as  the  lowest  loss  mode, 
e.g.,  the  mode  with  the  smallest  Im(^))  eigenvalue  as  a  function  of  which  is 

the  separation  between  the  VCSEL-AL  boundary  and  the  oxide  aperture.  I  tested  two 
different  oxide  aperture  radii,  0.5  pm  and  1.0  pm,  using  my  1.5A  (=  0.435  pm),  second  order 

polynomial  AL  design.  The  resulting  real  and  imaginary  parts  of  the  eigenvalue  are  shown  in 

Figures  D.3  and  D.4,  respectively.  Clearly,  the  eigenvalues  converge  for  pp^  -  p^^  >  0.9  pm 
(corresponding  to  i=s  3A),  denoting  the  minimum  allowable  separation  between  the  oxide 
aperture  and  Fs.  It  is  interesting  to  note  that  for  the  1.0  pm  oxide  aperture,  the  variation 
in  the  eigenvalue  from  p^^,  =  1.2  pm  to  2.0  pm  corresponds  to  a  variation  in  the  resonant 
wavelength  of  less  than  two  angstroms,  implying  that  the  AL  works  well  even  when  placed 
very  close  to  the  aperture.  No  analogous  test  was  performed  on  the  axial  AL  design;  I  assume 
that,  based  on  the  good  results  of  the  radial  AL  design  and  the  high  reflectance  of  the  oxide 
DBRs,  the  optimized  radial  AL  design  should  work  for  the  axial  AL  as  well.  To  perform 
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Reflectance 


1 


Figure  D.2  Calculated  intensity  reflection  R  =  |rp  of  the  radial  absorbing  layer  (AL) 
using  the  rigorous  transfer  matrix  approach.  The  calculation  is  based  on  a 
1.5A  thick  AL  at  870  nm,  positioned  with  =  2.0  fim.  Thirty  layers  were 
used  to  discretize  the  AL  loss  function.  Inset  is  the  intensity  reflection  in 
decibels  R^b  =  20logio(|»'|). 
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the  lasing  mode  analysis  of  Chapter  IV,  I  used  =  2.0  nm  and  varied  the  oxide  aperture 
radius  (p^^)  from  0.4 //m  to  1.0 /^m. 


D-8 


Figure  D.3  Real  part  of  the  eigenvalue  ^  vs  separation  between  oxide  radius  (Pox) 

radial  AL  radius  (prs)  lA-lTHIN  structure  for  two  different  oxide 

radii,  0.5  pm  and  1.0  pm. 
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Figure  D.4  Imaginary  part  of  the  eigenvalue  ^  vs  separation  between  oxide  radius  {po^) 
and  radial  AL  radius  (pps)  lA-lTHIN  structure  for  two  different  oxide 

radii,  0.5  //m  and  1.0  pm. 


D-10 


Bibliography 


1.  Abramowitz,  M.  and  LA.  Stegun,  editors.  Handbook  of  Mathematical  Functions  with 
Formulas,  Graphs,  and  Mathematical  Tables.  National  Bureau  of  Standards  Applied 
Mathematics  Series  55,  Department  of  Commerce,  1964. 

2.  Agrawal,  G.P.  and  N.K.  Dutta.  Long-Wavelength  Semiconductor  Lasers.  Van  Nostrand 
Reinhold,  1986. 

3.  Arfken,  G.  Mathematical  Methods  for  Physicists  (3rd  Edition).  Academic  Press,  1958. 

4.  Baba,  T.,  et  al.  “Spontaneous  Emission  Factor  of  a  Microcavity  DBR  Surface- Emitting 
Laser,”  IEEE  Journal  of  Quantum  Electronics,  ;27;1347-1358  (1991). 

5.  Baba,  T.,  et  al.  “Spontaneous  Emission  Factor  of  a  Microcavity  DBR  Surface  Emitting 
Laser  (II)  -  Effects  of  Electron  Quantum  Confinements,”  IEEE  Journal  of  Quantum 
Electronics,  ;25(5):1310-1319  (1992). 

6.  Babic,  D.L,  et  al.  “Room-Temperature  Continuous- Wave  Operation  of  1.54-//m  Vertical- 
Cavity  Lasers,”  IEEE  Photonics  Technology  Letters,  7(11):1225-1227  (November  1995). 

7.  Balanis,  C.A.  Advanced  Engineering  Electromagnetics.  John  Wiley  and  Sons,  1989. 

8.  Bjork,  G.,  et  al.  “Spontaneous-Emission  Coupling  Factor  of  Planar  Dielectric  Microcav¬ 
ity  Lasers,”  Physical  Review  A,  J7(5):4451-4463  (1993). 

9.  Bjork,  G.,  et  al.  “Definition  of  a  Laser  Threshold,”  Physical  Review  A,  5<?(2):1675-1679 
(1994). 

10.  Bjork,  G.,  et  al.  “Modification  of  Spontaneous  Emission  Rate  in  Planar  Dielectric  Mi¬ 
crocavity  Structures,”  Physical  Review  A,  .^.^(1):669-681  (1991). 

11.  Bjork,  G.  and  Y.  Yamamoto.  “Analysis  of  Semiconductor  Microcavity  Lasers  using  Rate 
Equations,”  IEEE  Journal  of  Quantum  Electronics,  ;27(ll):2386-2396  (November  1991). 

12.  Burak,  D.  and  R.  Binder.  “Cold-Cavity  Vectorial  Eigenmodes  of  VCSELs,”  IEEE  Jour¬ 
nal  of  Quantum  Electronics,  55(7):1205-1215  (July  1997). 

13.  Cendes,  Z.J.  “Vector  Finite  Elements  for  Electromagnetic  Field  Computation,”  IEEE 
Transactions  on  Magnetics,  ;27(5):3953-3966  (September  1991). 

14.  Choquette,  K.D.,  et  al.  “Scalability  of  Small-Aperture  Selectively  Oxidized  Vertical 
Cavity  Lasers,”  Applied  Physics  Letters,  70(7):823-825  (February  1997). 

15.  Choquette,  K.D.,  et  al.  “Fabrication  and  Performance  of  Selectively  Oxidized  Vertical- 
Cavity  Lasers,”  IEEE  Photonics  Technology  Letters,  7(11);1237-1239  (1995). 

16.  Cohen-Tannoudji,  C.,  et  al.  Quantum  Mechanics,  II.  John  Wiley  and  Sons,  1977. 

17.  Coldren,  L.A.,  et  al.  “Recent  Advances  and  Important  Issues  in  Vertical- Cavity  Lasers.” 
Vertical- Cavity  Surface- Emitting  Lasers.  Number  3003  in  Proceedings  of  the  SPIE.  2-13. 
1997. 


BIB-1 


18.  DeMartini,  F.,  et  al.  “Thresholdless  Microlaser,”  Physical  Review  A,  .^^(7):4220-4233 
(October  1992). 

19.  Deppe,  D.G.,  et  al.  “Oxide- Confined  VCSELs  with  Quantum  Well  and  Quantum  Dot 
Active  Regions.”  LEOS  ’97  Conference  Proceedings.  287-288.  November  1997. 

20.  Deppe,  D.G.,  et  al.  “Low- Threshold  Vertical-Cavity  Surface-Emitting  Lasers  Based  on 
Oxide- Confinement  and  High  Contrast  Distributed  Bragg  Reflectors,”  IEEE  Journal  of 
Selected  Topics  in  Quantum  Electronics,  5(3):893-904  (June  1997). 

21.  Deppe,  D.G.,  et  al.  “Eigenmode  Confinement  in  the  Dielectrically  Apertured  Fabry- 
Perot  Microcavity,”  IEEE  Photonics  Technology  Letters,  5(6):713-715  (June  1997). 

22.  Feld,  S.,  et  al.,  “Microcavity  Surface  Emitting  Lasers.”  Patent  Disclosure,  1997.  AF  Inv 
D00227. 

23.  Gelfand,  I.M.  and  S.V.  Fomin.  Calculus  of  Variations.  Prentice-Hall,  Inc.,  1963. 

24.  Gloge,  D.  “Weakly  Guiding  Fibers,”  Applied  Optics,  10:2252-2258  (1971). 

25.  Graglia,  R.D.,  et  al.  “Higher  Order  Interpolatory  Vecgtor  Bases  for  Computational 
Electromagnetics,”  IEEE  Transactions  on  Antennas  and  Propagation,  45(3):329-3A2 
(March  1997). 

26.  Hadley,  G.R.  “Effective  Index  Model  for  Vertical-Cavity  Surface-Emitting  Lasers,”  Op¬ 
tics  Letters,  ;20(13):1483-1485  (July  1995). 

27.  Hadley,  G.R.,  et  al.  “Understanding  Waveguiding  in  Vertical-Cavity  Surface- Emitting 
Lasers.”  Technical  Digest,  Conference  on  Lasers  and  Electro-optics.  1996  Technical 
Digest  Series.  425.  May  1996. 

28.  Harrington,  R.F.  Time-Harmonic  Electromagnetic  Fields.  New  York:  McGraw-Hill 
Book  Company,  1961. 

29.  Hegblom,  E.R.,  et  al.  “Scattering  Losses  from  Dielectric  Apertures  in  Vertical-Cavity 
Lasers,”  IEEE  Journal  of  Selected  Topics  in  Quantum  Electronics,  5(2):379-389  (April 
1997). 

30.  Hegblom,  E.R.,  et  al.  “Towards  Microcavity  Vertical  Cavity  Lasers:  Aperture  and 
Cavity  Design  for  High  Efficiency  and  Low  Threshold.”  IEEE  Summer  Topical  Meetings: 
Vertical- Cavity  Lasers.  August  1997. 

31.  Huffaker,  D.L.,  et  al.  “Native-Oxide  Defined  Ring  Contact  for  Low  Threshold  Vertical- 
Cavity  Lasers,”  Applied  Physics  Letters,  65(l):97-99  (July  1994). 

32.  Huffaker,  D.L.,  et  al.  “Threshold  Characteristics  of  Planar  and  Index-Guided  Microcav¬ 
ity  Lasers,”  Applied  Physics  Letters,  ^7(l):4-6  (July  1995). 

33.  Huffaker,  D.L.,  et  al.  “Low  Threshold  Half-Wave  Vertical-Cavity  Lasers,”  Electronics 
Letters,  SO (23):  1946-1947  (November  1994). 

34.  Ikegami,  T.  “Reflectivity  of  Mode  at  Facet  and  Oscillation  Mode  in  Double- 
Heterostructure  Injection  Lasers,”  IEEE  Journal  of  Quantum  Electronics,  QE-8 (6):470- 
476  (June  1972). 


BIB-2 


35.  Itoh,  T.,  et  al.,  editors.  Finite  Element  Software  for  Microwave  Engineering.  John  Wiley 
and  Sons,  Inc.,  1996. 

36.  Jackson,  J.D.  Classical  Electrodynamics  (2nd  Edition).  John  Wiley  and  Sons,  1975. 

37.  Jin,  J.  The  Finite  Element  Method  in  Electromagnetics.  John  Wiley  and  Sons,  Inc., 
1993. 

38.  Kardontchik,  J.E.  “Mode  Reflectivity  of  Narrow  Stripe- Geometry  Double  Heterostruc¬ 
ture  Lasers,”  IEEE  Journal  of  Quantum  Electronics,  QE- 1 8 (8)’.1279— 1286  (August 
1982). 

39.  Kendall,  P.C.,  et  al.  “Theory  for  Calculating  Approximate  Values  for  the  Propaga¬ 
tion  Constants  of  an  Optical  Rib  Waveguide  by  Weighting  the  Refractive  Indices,”  lEE 
Proceedings,  ISf  pt  A(8):699-702  (September  1987). 

40.  Lee,  J.-F.,  et  al.  “Tangential  Vector  Finite  Elements  for  Electronagnetic  Field  Compu¬ 
tation,”  IEEE  Transactions  on  Magnetics,  ;27(5):4032-4035  (September  1991). 

41.  Lewin,  L.  “A  Method  for  the  Calculation  of  the  Radiation-Pattern  and  Mode-Conversion 
Properties  of  a  Solid-State  Heterojunction  Laser,”  IEEE  Transactions  on  Microwave 
Theory  and  Techniques,  MTT-23 {7):576-58b  (July  1975). 

42.  Loehr,  J.P.  Physics  of  Strained  Quantum  Well  Lasers.  Kluwer  Academic  Publishers, 
1997. 

43.  Loehr,  J.P.,  et  al.  “Effects  of  Native  Oxides  and  Optical  Confinement  on  Microcav¬ 
ity  VCSEL  Spontaneous  Emission.”  Physics  and  Simulation  of  Optoelectronic  Devices 
VI 3283.  Proceedings  of  the  SPIE.  January  1998. 

44.  Lott,  J.A.  and  M.J.  Noble.  “Model  of  Intra  and  Extracavity  Photodetection  for  Planar 
Resonant  Cavity  Light  Emitting  Diodes.”  24th  International  Symposium  on  Compound 
Semiconductors.  1997. 

45.  Maloney,  J.,  et  al.  “Generalization  of  PML  to  Cylindrical  Geometries.”  13th  Annual 
Review  of  Progress  in  Applied  Computational  Electromagnetics.  900-908.  1997. 

46.  The  MathWorks,  Inc.  Partial  Differential  Equation  Toolbox  2.0:  User’s  Guide,  1997. 
Documentation  for  Beta  test. 

47.  Meystre,  P.  and  M.  Sargent  III.  Elements  of  Quantum  Optics.  Springer- Verlag,  1991. 

48.  Michalzik,  R.  and  K.J.  Ebeling.  “Generalized  BV  Diagrams  for  Higher  Order  Transverse 
Modes  in  Planar  Vertical- Cavity  Laser  Diodes,”  IEEE  Journal  of  Quantum  Electronics, 
57:1371-1379  (1995). 

49.  Nedelec,  J.  C.  “Mixed  Finite  Elements  in  R®,”  Numerische  Mathematik,  55:315-341 
(1980). 

50.  Noble,  M.J.,  et  al.  “Analysis  of  Microcavity  VCSEL  Lasing  Modes  using  a  Full  Vector 
Weighted  Index  Method,”  IEEE  Journal  of  Quantum  Electronics  (1997).  In  Review. 


BIB-3 


51.  Noble,  M.J.,  et  al.  “Calculation  of  Microcavity  VCSEL  Field  Modes  using  a  Doubly 
Iterative  Weighted  Index  Method.”  Physics  and  Simulation  of  Optoelectronic  Devices 
V.  Number  2994  in  Proceedings  of  the  SPIE.  259-266.  February  1997. 

52.  Noble,  M.J.,  et  al.  “Calculation  of  VCSEL  lasing  mode  threshold  gain  with  the  weighted 
index  method.”  LEOS  ’97  Conference  Proceedings 2.  354-355.  November  1997. 

53.  Noble,  M.J.,  et  al.  “Optimization  of  Microcavity  VCSELs.”  Vertical  Cavity  Surface 
Emitting  Lasers  II.  Number  3286  in  Proceedings  of  the  SPIE.  January  1998. 

54.  Noble,  M.J.,  et  al.  “Analysis  of  Three  Dimensionally  Confined  Microcavity  Surface  Emit¬ 
ting  Lasers  using  Vector  Finite  Elements.”  IEEE  Summer  Topical  Meetings:  Vertical- 
Cavity  Lasers.  August  1997. 

55.  Noble,  M.J.,  et  al.  “Design  of  Low- Threshold,  Single- Transverse  Mode  VCSELs  using  the 
Transverse  Confinement  Optical  Mode  Separation.”  Advances  in  Vertical  Cavity  Surface 
Emitting  Lasers  15.  OSA  Trends  in  Optics  and  Photonics,  edited  by  C.  Chang-Hasnain, 
Optical  Society  of  America,  1997. 

56.  Noble,  M.J.,  et  al.  “Rigorous  Optical  Analysis  of  Oxide- Apertured  Microcavity  VCSELs 
using  Vector  Finite  Elements,”  IEEE  Journal  of  Quantum  Electronics  (1998).  In  Review. 

57.  Noble,  M.J.,  et  al.  “Calculation  and  Measurement  of  Resonant-Mode  Blueshifts  in  Oxide- 
Apertured  VCSELs,”  IEEE  Photonics  Technology  Letters.,  70(4)  (April  1998). 

58.  Okoshi,  T.  Optical  Fibers.  New  York:  Academic  Press,  1982. 

59.  Polstyanko,  S.V.  and  J.-F.  Lee.  “Hi(cwr/)  Tangential  Vector  Finite  Element  Method  for 
Modeling  Anisotropic  Optical  Fibers,”  Journal  of  Lightwave  Technology,  75(11):2290- 
2295  (November  1995). 

60.  Pozar,  D.M.  Microwave  Engineering.  Addison- Wesley,  1990. 

61.  Robertson,  M.J.,  et  al.  “The  Weighed  Index  Method:  A  New  Technique  For  Analyzing 
Planar  Optical  Waveguides,”  IEEE  Journal  of  Lightwave  Technology,  7(12):2105-2108 
(December  1989). 

62.  Saad,  Y.  Numerical  Methods  for  Large  Eigenvalue  Problems.  Manchester  University 
Press,  1992. 

63.  Sacks,  Z.S.,  et  al.  “A  Perfectly  Matched  Anisotropic  Absorber  for  Use  as  an  Absorbing 
Boundary  Condition,”  IEEE  Transactions  on  Antennas  and  Propagation,  J5:1460-1463 
(1995). 

64.  Sadiku,  M.N.O.  Numerical  Techniques  in  Electromagnetics.  CRC  Press,  1992. 

65.  Saleh,  B.E.A.  and  M.C.  Teich.  Fundamentals  of  Photonics.  John  Wiley  and  Sons,  Inc., 
1991. 

66.  Shin,  H.-E.,  et  al.  “780  nm  Oxidised  Vertical-Cavity  Surface- Emitting  Lasers  with 
Alo.11Gao.89As  Quantum  Wells,”  Electronics  Letters,  S:2(14):1287-1288  (July  1996). 

67.  Shin,  J.-H.,  et  al.  “Very  Small  Oxide-Confined  Vertical  Microcavity  Lasers  with  High- 
Contrast  AlGaAs/Alj;Oy  Mirrors,”  IEEE  Photonics  Technology  Letters  (June  1998). 


BIB-4 


68.  Silvester,  P.P.  and  R.L.  Ferrari.  Finite  Elements  for  Electrical  Engineers  (3rd  Edition). 
Cambridge  University  Press,  1996. 

69.  Snitzer,  E.  “Cylindrical  Dielectric  Waveguide  Modes,”  Journal  of  the  Optical  Society  of 
America,  (5):491-498  (May  1961). 

70.  Snyder,  A.W.  and  J.D.  Love.  Optical  Waveguide  Theory.  Chapman  and  Hall,  1983. 

71.  Soda,  H.,  et  al.  “GalnAsP/InP  Surface  Emitting  Injection  Lasers,”  Japanese  Journal  of 
Applied  Physics,  75(12):2329-2330  (December  1979). 

72.  Sun,  D.,  et  al.  “Spurious  Modes  in  Finite-Element  Methods,”  IEEE  Antennas  and 
Propagation  Magazine,  57(5):12-24  (October  1995). 

73.  Vurgaftman,  1.  and  J.  Singh.  “Spatial  and  Spectral  Characteristics  of  Spontaneous  Emis¬ 
sion  from  Semiconductor  Quantum  Wells  in  Microscopic  Cylindrical  Cavities,”  Applied 
Physics  Letters,  ^7(26):3865-3867  (December  1995). 

74.  Wong,  M.F.,  et  al.  “Axisymmetric  Edge-Based  Finite  Element  Formulation  For  Bodies 
Of  Revolution:  Application  To  Dielectric  Resonators.”  IEEE  Microwave  Theory  and 
Techniques  (MTT-S)  Digest.  285-288.  1995. 

75.  Yamamoto,  Y.,  et  al.  “Microcavity  Semiconductor  Laser  with  Enhanced  Spontaneous 
Emission,”  Physical  Review  A,  (l):657-667  (July  1991). 

76.  Yang,  G.M.,  et  al.  “Ultralow  Threshold  Current  Vertical-Cavity  Surface- Emitting  Lasers 
Obtained  with  Selective  Oxidation,”  Electronics  Letters,  5i(ll):886-888  (1995). 

77.  Yokoyama,  H.  and  K.  Ujihara,  editors.  Spontaneous  Emission  and  Laser  Oscillation  in 
Microcavities.  CRC  Press,  1995. 


BIB-5 


Vita 


Captain  Noble  was  born  on  18  Aug  1967  in  Ft.  Meade,  MD.  He  attended  Rensselaer 
Polytechnic  Institute  where  he  graduated  in  1989  with  a  Bachelor  of  Science  Degree  in 
Electrical  Engineering.  Upon  graduation,  he  accepted  a  commission  in  the  United  States 
Air  Force.  In  July  1989,  Captain  Noble  reported  to  his  first  duty  station  at  Norton  AFB,  CA. 
During  his  four  year  tour  at  Norton  AFB  Captain  Noble  performed  acquisition  management 
and  engineering  duties  at  the  Ballistic  Missile  Organization.  In  May  1993,  Captain  Noble 
was  assigned  to  the  Air  Force  Institute  of  Technology,  where  he  has  earned  a  Masters  of 
Science  degree  in  Electrical  Engineering  in  1994,  and  continued  on  to  pursue  the  Doctor  of 
Philosophy  degree. 


Permanent  address:  5117  Pumphery  Drive 

Fairfax,  Virginia  22032 


VITA-1 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704^0188 


1 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing 
the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information 
Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

18  Mar  98  Dissertation 

4.  TITLE  AND  SUBTITLE 

Optical  Physics  of  Microcavity  Surface  Emitting  Lasers 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

M.J.  Noble 

7.  PERFORMING  ORGANIZATION  NAMEIS)  AND  ADDRESS(ES) 

AFIT/ENG,  2950  P  Street,  Bldg  640,  WPAFB,  OH  45433 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

AFIT/DS/ENG/98-03 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Dr.  John  P.  Loehr 

AFRL/AADP,  Bldg  620,  WPAFB,  OH  45433 
(937)  255-2227  x3388 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

IZa.  DISTRIBUTION  AVAILABILITY  STATEMENT 

UNLIMITED 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

This  dissertation  consists  of  an  in-depth  theoretical  analysis  of  the  optical  impact  of  transverse  index  confinement  on  the 
lasing  modes  of  microcavity  surface  emitting  lasers.  Two  different  variational  techniques  were  employed  to  calculate  the 
relevant  optical  laser  parameters;  modal  resonance,  field  profile,  confinement  factor,  and  threshold  gain. 

Through  a  semianalytic  field  analysis,  two  new  confinement  related  effects  were  discovered:  a  polarization  dependent  change 
in  mirror  reflectance  with  aperture  radius,  and  a  mode  dependent  blueshift.  Through  a  numerical  field  analysis,  the  physical 
mechanisms  governing  diffraction  were  ascertained.  These  are:  the  transverse  confinement  of  the  optical  mode,  the  degree 
of  tilt  of  the  mode  propagation  vector  away  from  normal  to  the  mirror  surfaces,  and  the  density  of  parasitic  modes  in  the 
spectral  vicinity  of  the  lasing  mode.  Based  on  these  results,  a  new  semianalytic  technique  was  derived  for  estimating 
diffiraction  loss  via  rigorous  coupling  to  the  parasitic  modes.  Finally,  this  model  was  combined  with  a  gain  model  and  a  new 
mode  density  model  to  compute  spontaneous  emission  rate,  light  verses  current,  threshold  current,  and  spontaneous  emission 
factors  for  a  low  threshold  oxide-apertured  device. 

14.  SUBJECT  TERMS 

Microelectronics,  Semiconductor  Lasers 

15.  NUMBER  OF  PAGES 

188 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION 

OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED  UNCLASSIFIED 

20.  LIMITATION  OF 

ABSTRACT 

UL 

Standard  Form  298  (Rev.  2*89)  (EG) 

Prescribed  by  ANSI  Std.  239.18 

Designed  using  Perform  Pro,  WHS/DIOR,  Oct  94 


